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2.0 INTRODUCTION 

Measuring trends in ozone, and most other geophysical variables, requires that a small systematic 
change with time be determined from signals that have large periodic and aperiodic variations. Their time 
scales range from the day-to-day changes due to atmospheric motions through seasonal and annual varia- 
tions to 1 1 -year cycles resulting from changes in the sun ultraviolet output. Aperiodic variations include 
the irregular quasi-biennial oscillation, with a period of roughly 26-28 months, approximately 4-year 
variations, and other sources of interannual differences. 

Because the magnitude of all of these variations is not well known and highly variable, it is necessary 
to measure over more than one period of the variations to remove their effects. This means at least 2 01 - 
more times the 1 l-year sunspot cycle. Thus, the first requirement is for a long-term data record. The second 
related requirement is that the record be consistent; a small effect is being sought, and changes in instru- 
mentation or data analysis method will obscure changes in the atmosphere. A third requirement is tor 
reasonable global sampling, to ensure that the effects are representative of the entire Earth. Therefore, the 
various observational methods relevant to trend detection are reviewed in Section 2. 1 to characterize tneir 
quality and time and space coverage. Available data are then examined for long-term trends or recent 
changes in ozone total content (Section 2.2) and vertical distribution (2.3), as well as in related parameters 
such as stratospheric temperature (2.4), source gases and tropospheric ozone (2.5), and aerosols (2.6). The 
relation between trends in total column ozone and variations observed in the solar ultraviolet radiation at 
the ground are discussed in Section 2.7, and outstanding issues in trends detection are emphasized in 

Section 2.8. 


2.1 OBSERVATIONAL METHODS RELEVANT TO TREND DETECTION 

2.1.1 Introduction 

This section briefly describes the measurement techniques that have been used in the past or are 
expected to be used in the future to quantitatively measure trends in total ozone or the ozone vertical 
distribution. More detailed descriptions of these techniques have been presented elsewhere, in the Ozone 
Trends Panel Report (WMO, 1989, hereafter denoted OTP) and in previous reports (WMO, 1986; NASA, 
1988), earlier ozone assessments, and in scientific journals and reports. The purpose of this section is to 
describe these observations briefly and to characterize their quality, as well as their time and space coverage . 

For the individual measurements, the most critical quantity is stability, or the absence of time- 
dependent systematic errors. A systematic error, i.e., an error that is consistently present, is the amount 
by which the mean of a large number of individual observations of the same ozone value could differ from 
the “true” value. This is also referred to as the accuracy. The relevance of the accuracy to trend deter- 
mination is further discussed in .Section 2.1.6. 

In the following sections, the term precision is used to refer to the random variations or spread of 
values that an instrument would report when observing a constant ozone value. It is sometimes referred 
to as repeatability, and depends on the random errors of the measuring system. This is important for many 
studies of atmospheric processes, and could be important for trend studies when only a few observations 
of a particular kind are available. However, most quantities compared in trend studies involve averaging 
a large amount of data, reducing the spread, so precision is often not of primary importance. 
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The section is organized as follows: Section 2.1.2 describes the measurements that have resulted in 
the data that have been used to derive trends, or closely related observations. The techniques that are 
expected to add to the trend data for future analyses are presented in Section 2.1.3. Both of these sections 
are divided into ground-based observations (including balloons or rockets launched from the ground) and 
satellite observations. Intercomparisons of different types of data, and what they indicate about the 
capabilities of the different techniques, are described in 2.1.4, which is divided into sections on present 
and future trend data. This portion of the chapter concludes with brief discussions of the sampling 
requirements in time and space (2. 1 .5) and an overview of trend-measuring capabilities (2. 1 .6). 

2.1.2 Measurements Used in Current Trend Analyses 

Measurements of ozone in the atmosphere have been made for many years. However, most of the 
data that have the length of record, accuracy, and stability to allow a search for trends only go back to 
about the time of the IGY (International Geophysical Year) in 1957. Figure 2.1-1 gives a time-line for the 
types of measurement that have been used for trend analyses, and the periods involved. These measure- 
ments are outlined below. 

2.1. 2.1 Total Ozone Measurements 
Ground-based Measurements 


Absorption spectroscopy affords a sensitive means for monitoring the column abundance of ozone 
from the ground. The attenuation of monochromatic radiation is related to the number of absorbing 


OZONE MEASUREMENTS USED IN TREND ANALYSES 
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Figure 2.1-1. Time lines of measurements available for use in ozone trend analyses at the present time. 
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molecules in the optical path as they undergo state transitions that absorb incoming radiation. Absorption 
features in the ultraviolet and visible regions of the spectrum mainly involve electronic transitions and 
hence are generally relatively insensitive to temperature and pressure, which is a feature that considerably 
simplifies the reduction of the data to slant-column abundances. 

The source of incoming light may be a celestial body, usually the sun, but also the moon, or a star, 
or it may be the light scattered from the sunlit zenith sky. Direct-light absorption measurements can be 
carried out only during clear periods from the ground, since an unobstructed view of the light source is 
required. Scattered light measurements can be obtained in cloudy weather by some absorption systems. 


Dobson Spectrophotometer 

The standard instrument in the Global Ozone Observing System is the Dobson spectrophotometer 
(Dobson, 1957). The instrument is a quartz double monochromator which measures the relative intensities 
of pair wavelengths in the Huggins ozone absorption band (300-350 nm) from which total ozone in a vertical 
column of the atmosphere can be deduced. Most precise observations are made on direct sunlight and 
standard double-pair wavelengths designated AD (305.5/325.4 and 317.6/339.8 nm). Less precise observa- 
tions are made on clear or cloudy zenith skylight or, infrequently, moonlight. (For a more detailed 
description of the Dobson spectrophotometer and its operation, see OTP.) Direct sun AD observations 
generally provide the most precise measurements if the secant of the zenith angle, p, is less t an . 

Long-term ozone measurement precision for the Dobson spectrophotometer, for annual means, is 
estimated to be 1% (at the 2 a level), based on the standard deviation <r from analyses ot mean data from 
individual stations (WMO, 1980, 1981, 1982). Attainment of this precision requires that observations at all 
times be made on correct wavelengths; that changes in spectral characteristics ot the instrument be 
accounted for using data derived from periodic standard lamp tests; that recalibrat.on of the wedge be 
performed, if evidence suggests it is needed; that periodic recahbrations of the instrument be performed 
by the Langley method (Dobson and Normand, 1962) or through intercalibration with a primary standard 
Dobson instrument (Komhyr et al„ 1989); that observational errors be minimized; that the temperature 
dependence of the ozone absorption coefficients be taken into account; and that relations derived empiri- 
cally between direct sun and clear or cloudy zenith sky observations be adequately quantified. For use ot 
the total ozone data in global trend analyses, it is necessary, furthermore, that observations be made on a 
sufficient number of days each month in order to obtain representative monthly mean data (see Section 
2.1.5), and that the observations not be unduly influenced by local interfering absorbing species, such as 
sulfur dioxide and nitrogen dioxide or ozone produced photochemically in locally polluted air. 

Until 1968, Dobson instrument calibrations at different stations were generally conducted randomly 
and independently. From 1974 onwards, increasing numbers of instruments have been modernized, re ur- 
bished, and calibrated by direct intercomparison with the WMO designated World Primary Standard 
Dobson Spectrophotometer No. 83, maintained at the World Dobson Spectrophotometer Centra a o- 
ratory in Boulder, Colorado. In subsequent intercomparisons, these instruments have shown typical 
calibration changes of 0-2% for direct sun observations on AD wavelengths. The long-term (196. 1987) 
ozone measurement stability of Primary Standard Dobson Instrument No. 83 is reported to have been 
maintained at ±0.5% (Komhyr et al„ 1989). Since the mid-1970s, virtually all (90) Dobson instruments of 
the global Dobson instrument station network have been calibrated several times either directly with 
instrument No. 83 or indirectly through intercalibrations with Secondary Standard Dobson spectrop oto- 
meters calibrated in Boulder in 1977. 
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The uneven geographical distribution of the existing Dobson spectrophotometer network (Figure 2.1- 
2) gives rise to a spatial sampling error when attempts arc made to determine global ozone content and 
trends. At most, the Dobson instruments can provide trend information for specific regions of the globe. 
An important utility of the Dobson instruments is their ability to provide correlative data for satellite 
instruments that measure ozone on a global scale (Fleig et al., 1986; 1989a; 1989b; Bojkov et al„ 1988; 
McPeters and Komhyr, 1989). but that are subject to calibration drifts and are not highly sensitive to 
tropospheric ozone. 

Filter Ozonometers M-83 and M-124 

Since 1957, routine ground-based total column ozone measurements have been made at more than 40 
stations in the USSR using a filter ozonometer instrument designated as type M-83 (not to be confused 
with the Standard Dobson ozone spectrophotometer, which was the 83rd instrument manufactured of the 
Dobson type). The filter-type instrument is based upon the same principle as the Dobson spectrophotometer 
in using differential absorption of ultraviolet radiation in the 300-350 nm Huggins band of ozone. The 
M-83 instrument, however, uses two broadband filters and measures the relative attenuation of the solar 
ultraviolet radiances either directly from the sun or indirectly from the zenith sky (Gustin, 1963). 

Direct intercomparisons between M-83 filter instruments and Dobson spectrophotometers prior to 
1971 (Bojkov, 1969) revealed that the M-83 records show 6% less ozone when the observations are restricted 
to an air mass p < 1 .5, and 20% to 30% more ozone when data are taken for p > 2.0. A strong dependence 
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on turbidity was also detected, with 9% to 14% higher ozone readings when the surface visibility was less 
than 5 km. These strong deviations for » > 2.0 make many of the high-lat.tudc measurements in the L SSR 
very uncertain, especially for measurements prior to 1972, as described below. 

Improved filters were introduced into the M-83 instrument starting in 1972-1973 (Gustin, 1978). The 
new filters have maximum transmittance at 301 nm and 326 nm, and their band passes are less than those 
in the earlier version: 22 nm (291-312 nm) and 15 nm (319-334 nm). Comparison of N.mbus-4 BU V satellite 
overpasses over M-83 stations in the USSR demonstrated a standard deviation of about 50 Dobson Units 
before 1973 and about 25 DU afterward (WMO, 1980. 1983). The Nimbus-4 BUV overpasses of Dobson 
stations maintained a standard deviation of about 17 DU during the 1970-1977 lifetime ot the satellite. 

A much newer, reportedly improved instrument designated as M-124 has been installed in many 
stations since 1986 (Gustin and Sokolenko. 1985), but no ozone data have been reported yet tor this 
instrument. No trend data with the M-124 can be expected for about a decade unless the data can be 
satisfactorily cross-calibrated with the M-83 data from the same location. 


Satellite observations 


Total Ozone Monitoring Spectrometer (TOMS) 


The TOMS, launched on the Nimbus-7 spacecraft in 1978, is an instrument (Heath et al., 1975) whose 
primary measurement goal is to obtain contiguous mapping of the total column ozone amount over the 
globe (Bowman and Krueger. 1985; Schoeberl et al., 1986). To achieve this, TOMS step scans across the 
sub-orbital track, sampling radiation backscattered from the underlying surface and atmosphere. Ozone 
column amounts are inferred by utilizing the wavelength dependence of the Earth's ultraviolet albedo in 
the Huggins band of the ozone absorption spectrum. The TOMS raw data are measurements of the direc 
and backscattered solar U V radiation at six fixed wavelength channels (312.5, 317.5, 331.2 339.8 360 and 
380 nm). Data from the first four channels are used in pairs to provide three estimates ot the total column 
ozone amount by the differential absorption method. The remaining two channels, which are free ot ozone 
absorption, are used to determine the effective background albedo. 


It was recognized from the outset that this technique was intrinsically capable ot very high accuracy 
and stability, since the requirement was for a relative measurement of the ratio of Earth's backscattered 
UV radiance to the solar UV irradiance at the same wavelength. Because both measurements could in 
principle be made with the same instrument, the determination of albedo as a function of wavelength should 
not depend on either the absolute calibration of the instrument nor on long-term variations in the sensitivity 
of the instrument. However, a serious uncertainty is introduced by the use of the ditfuser plate, which is 
not common to both measurements but only used to transform the solar irradiance into a radiance that is 
comparable in magnitude to the backscattered Earth radiance, and can be measured in the same manner. 
Observations of the solar irradiance indicated that the diffuser plate was degrading and becoming less 
reflective with time. Cebula et al. (1988) developed a model of the diffuser and instrument degradation 
which was adopted by the Ozone Processing Team (OPT). Application of this model in the reduction of 
the TOMS data led to a downward drift of the TOMS results compared to those of the ground-based 
Dobson network (see above) reported by Fleig et al. (1986, 1989a), Bojkov et al. (1988). 


Extensive analysis reported by the Trends Panel indicated that the OPT model had large uncertainties 
and had almost certainly underestimated the true diffuser plate degradation. The errors were large enoug 
to allow agreement with the Dobson data. Subsequently, some results for trends in total ozone in the 
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Trends Report were based on normalizing the TOMS data to the Dobson values. This allows extension to 

global coverage, but is at best only an approximation and provides no independent information from TOMS 
on the trends. 

Subsequently, a major advance in the analysis of the TOMS data has been made. Bhartia et al. (1989) 
pointed out that, for the related Solar Backscatter Ultraviolet (SBUV) experiment, there was a wavelength 
at 305.8 nm whose data could be used with those from the SBUV 312.5 nm channel to form a new “D 
pair. Results from this pair are less sensitive to diffuser degradation for two reasons. First, the degradation 
varies with wavelength, but the D pair wavelengths are only 6.7 nm apart, compared to the 18.7 for the A 
pair. In addition, because the difference in ozone absorption coefficients is larger for the D pair than for 
the other pairs, results are estimated to be only 0.22 times as sensitive to diffuser drift than the archived 
values, which are based on a latitude-dependent weighted sum of the A, B, and C pairs. Because of the 
large absorption coefficients. D pair values can only be obtained in the tropics, where the total column 
amount is small, but a comparison showed that they follow the Dobson values much more closely than the 
archived values do. This provides confirmation that the archived data are in error, as well as strong 
corroboration that the error mechanism has been identified (Bhartia et al., 1989; Hudson et al., 1989 
private communication). 

Generalizing from these insights, McPeters et al. (1989, private communication) have developed the 
pairjustification method. The basis is that different pairs have different sensitivity to diffuser degradation 
errors, resulting in a drift between pairs than can be measured accurately. By requiring that all pairs obtain 
the same total ozone result, and assuming that the uncorrected diffuser degradation is approximately 
linearly dependent on wavelength, a unique determination of the total ozone is obtained that is independent 
ot diffuser degradation. The implementation makes use of pairs denoted A, B, and C (312.5/331.2; 317 5/ 
339.8; 331.2/339.8 nm). There appears to be a fundamental limit of 0.5-1. 0% per decade to the stability 
that can be obtained. Test results are described in Section 2. 1 . 4 . 

Solar Backscatter Ultraviolet Spectrometer 


The Solar Backscatter Ultraviolet (SBUV) spectrometer (described below in Section 2. 1.2.2) also 
measures the solar ultraviolet radiation that is backscattered by the Earth and atmosphere. Data from four 
of its wavelength channels, which are the same as four of the TOMS channels, are used in the same way 
to determine total ozone amounts, but only in the nadir directly below the spacecraft. SBUV ozone values 
arc systematically slightly lower than those from TOMS, but they vary with time in the same way. Because 

SBUV also relies on the same diffuser plate, its total ozone values have also decreased with time due to 
diffuser degradation. 

2.1. 2.2 Profile Measurements 
Ground-based Measurements 

Standard Dobson Umkehr 

The Dobson instrument can also be used to obtain information on the vertical distribution of ozone. 
Measurements of the downward scattered radiation at two wavelengths are made for solar zenith angles Z 
fiom 60-86.5 degrees. From the ratio of these radiances as a function of Z, which reverses for Z near 90 
degrees, the vertical profile can be interred. (The method takes its name from the German word for a 
reversal.) The solution is usually given in terms of the layer mean ozone partial pressure (nb) of ozone in 
layers about 5-km thick, but with layer 1, the troposphere, about 10-km thick. 
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Although the retrieval provides data for Umkehr layers 1 through 9 (0-48 km), the averaging kernels 
for the present standard algorithm indicate that only data for layers 4 to 8 (19-43 km) should be used for 
trend analysis, while other layers should be used with caution. The vertical reso'ution of the retnevad for 
these layers is from 1 1 to 14 km. The systematic uncertainty or bias is estimated to be 5-20 ze in the (J 1 r , 
and 5-12% by DeLuisi et al. (1989a). Random errors, due to the dependence of the ozone absorption 
coefficients on the varying temperature, tropospheric and stratospheric aerosols, thin clouds invisible from 
the surface, and instrumental effects, are of the order of 5-10%, some of which could be reduced by 
monitoring the clarity of the zenith sky. Neither the systematic nor random errors should seriously attect 
estimation of long-term trends. Trend estimation will be affected by instrument drifts, by changes in 
calibration, by stratospheric aerosols, and by real atmospheric temperature trends. Calculations indicate 
that realistic temperature trends (< 2°C/decade) will induce errors < 0.2%/decade in ozone trends for layers 
4 through 8 (Figure 3-12, Chapter 3, OTP). Serious errors from stratospheric aerosols are easily recognized 
in the data and such contaminated data may be edited out before trend estimation. Moderate and smaller 
errors may be adequately corrected by the methods of DeLuisi et al. ( 1989a) or by statistical methods using 
proxy data (Reinsel et al., 1989). Finally, instrumental factors may be handled by statistical methods (e.g., 
sudden calibration changes as by Reinsel et al., 1989) or by re-evaluation of the observations using new 

calibration data. 


To reduce the time required for an Umkehr measurement, the Short method was developed and tested 
by Mateer and DeLuisi (1984). This method uses the A-C-D Dobson wavelength pairs in place of the C 
pair that is used for the Standard Umkehr method. By using the triple pair, ozone profile information is 
obtained when measurements are made during a solar zenith angle change of 80 to 89 degrees. In contrast, 
the C pair Standard Umkehr method that has provided most of the past data requires a solar zenith angle 
change of 60 to 90 degrees. An exploratory ozone profile retrieval algorithm for the Short Umkehr has 
been developed for testing. An operational algorithm will be developed after the present Standard Umkehr 
algorithm (Mateer and Duetsch, 1964) has been updated. 


Although Umkehr observations date back to the 1930s, the current Umkehr profile archive at the 
World Ozone Data Centre (viz., retrievals using the present standard C-Umkehr algorithm) begins about 
the time of the IGY (ca. July 1957). The archive comprises over 35,000 ozone profile retrievals for some 
68 stations, uniformly processed using the standard retrieval algorithm (Mateer and duetsch 1964). For 
26 of these stations, the entire station record consists offewer than 100 profiles. Reinsel et al. (1984), found 
only 13 stations with a sufficient number of observations and length of record for use in their trend analyses, 
while Reinsel et al. (1989) only used 10. One of the stations, Mont Louis, no longer makes > observations^ 
The record of Aspendale (the only Southern Hemisphere station in the trend set) ends in 198_. From 1983 
onward the record is being continued at Melbourne, 25 km away. The remaining trend stations cover the 
latitude range from 24 to 53° North. 


It appears that additional stations, especially those in the Automated Dobson Network (Komhyr et 
al., 1985), will have a sufficient number of observations and length of record to be used tor trend estimation 

within the next 5 years. There are presently seven automated Dobson instru ! k th! 
Standard and Short Umkehr measurements. These instruments were developed by NOAA/GMCC (Komhyr 
et al 1985) and are located at Perth (Australia); Lauder (New Zealand); Huancayo (Peru); Mauna Loa 
(Hawaii); Boulder, Colorado (United States); Haute Provence (France); and Fairbanks, Alaska (United 
States). The frequency of Umkehr measurements is significantly increased by the automated Dobson 
because observers are not required, and a shorter observing time (for the Short Umkehr) decreases the 
chance of cloud interference. Ancillary measurements of sky conditions are also made at the automated 
Dobson sites. These measurements consist of zenith-sky clouds and turbidity. In addition, lidar measure- 
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ments of stratospheric aerosols for correcting errors to the Umkehr profiles (DeLuisi et al., 1989a) are 
routinely made at Mauna Loa, Boulder, and Haute Provence. 

Balloon Ozonesondes 


Balloon ozonesondes are compact, lightweight, balloon-borne instruments flown with standard mete- 
orological radiosondes for measurement of ozone, air pressure, and air temperature up to altitudes of about 
30 km. Regener (1960, 1964) developed a very fast response chemiluminescent ozonesonde, but these 
instruments often exhibited considerable variations in sensitivity to ozone (which was later corrected). 
Thus, while an early measurement program provided a great deal of useful information on the ozone 
climatology over the North American continent, the variable response characteristics associated with the 
measurements and the relative shortness of the record render the data unsuitable for ozone trend studies 
in the troposphere and stratosphere. 

Brewer and Milford (1960) described an electrochemical ozone detector employing the well-known 
oxidation of potassium iodide (KI) by ozone as the basic reaction. Subsequently, Griggs (1961) investigated 
the physical and chemical aspects ot a similar balloon-borne instrument commonly referred to as the Brewer 
“bubbler” ozonesonde. Versions of the instrument are manufactured in the United States, G.D.R., and 
India. Because the air pumps of the Brewer ozonesondes are lubricated with a thin film of oil that may 
destroy ozone, they must be conditioned with ozone prior to flight time to minimize the ozone destruction. 
With proper conditioning, ozone losses can be kept to a few percent, as has been the case at Hohenpeis- 
senberg Observatory (F.R.G.). A few soundings from routine ozone measurement programs conducted in 
past years, however, have exhibited ozone losses, in extreme cases of up to 50%. Improvement in data 
quality is achieved through normalization of all soundings to Dobson spectrophotometer total ozone. The 
normalization factor for any sounding is a constant by which ozone values at all altitudes are multiplied. 
Use of a constant multiplication factor for normalization may, however, not be justified for soundings that 
exhibit large ozone losses within the instrument (Hilsenrath et al., 1986). Brewer ozonesonde data having 
normalization factors that range from 0.95-1.25 are generally acceptable; in an analysis later in this chapter 
the range 0.9-1. 2 is used. The average correction factor for more than 1,000 ozonesondes flown over 20 
years at Hohenpeissenberg is 1.07. Other uncertainties associated with ozone measurements with Brewer 
ozonesondes stem from variations in pump air flow rates above 10 mb, particularly for lightly lubricated 
pumps that become “dry” after several hours of operation . Build-up of Agl on the platinum anode during 
operation may also affect sensor performance. The significance of errors due to these effects has not been 
adequately assessed. 


A considerable amount of useful atmospheric ozone vertical distribution data has nevertheless been 
obtained with the Brewer ozonesondes. The data are suitable for ozone trend analyses, particularly from 
stations that have maintained unchanged instrument preflight conditioning and test procedures throughout 
the entire measurement program. 

Komhyr (1969) developed an electrochemical concentration cell (ECC) ozonesonde utilizing the 
reaction of ozone with KI, but with platinum cathode and anode electrodes contained in separate sensor 
chambers connected by an ion bridge. This sensor evolved from a carbon-iodine (Cl) ozonesonde (Komhyr, 
1964) that is no longer in use in the United States, but continues to be used in Japan. Because the chemical 
composition of the electrodes of the ECC sensor remain unchanged during operation, the sensor can be 
used indefinitely without deterioration of performance. The ECC sonde incorporates a miniature air pump 
fabricated from Teflon reinforced with glass fibers. The pump is not lubricated, so that minimal conditioning 
with ozone is required to prevent ozone loss within it. Early versions (type 3A) of the ECC ozonesonde 
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employed a pump of rectangular cross section. Because of non-uniformity of manufacture, some of the 
pumps exhibited excessive leakage at pressure altitudes above 15 mb, which in extreme cases have een 
underestimating the ozone by more than 15%. A newer instrument version (type 4A) employs a Teflon 
pump of circular cross section with external O-ring seals that is claimed to render the pumps and instruments 
suitable for use at higher altitudes. 


ECC ozonesondes using a 1% KI cathode electrolyte yield measured total ozone amounts that agree 
closely with Dobson spectrophotometer total ozone. For example, for 525 soundings made by NOAA since 
1984 over a wide range of operating conditions, from the tropics to the polar regions, the mean ECC sonde- 
Dobson spectrophotometer total ozone normalization factor was 1.01 ±0.05 (lo) (Komhyr et al., 1989) 
ECC ozonesondes are. therefore, suitable for use at stations where independent measurements ot total 
ozone are not available, e.g., in polar regions during the polar night. 

Ozone measurement uncertainty for ECC sondes is estimated to be ± 10% in the troposphere, ±5% 
in the stratosphere up to 10 mb, with the uncertainty increasing to 20±5%at 3 mb (Hi ■ senrath et al., 1986F 
Two ECC ozonesondes flown with an ozone UV-photometer during the MAP-GLOBUS campaign of 1983 
(Aimedieu et al. 1987) gave ozone values that differed from the photometer values by 2.1 ± : 1 .1% at 8.1 ± 1.1 
mb, and by -0.6±3.0% at 3.9±0.4 mb. The low ozone values measured above 10 mb, if real, may stem 
from application during data processing of inadequate pump efficiency corrections. Variability in the data 
at these altitudes is attributable at least in part to variations in manufacture of the pump components. It 
must be pointed out that small, continuing improvements over the years have led to more accurate 
measurements, but they may also result in spurious indications of trends. 


Rocket ozonesondes 

Over the last 2 to 3 decades, several groups in various countries have developed and used rocket- 
sondes In the USSR, rocket optical ozonesondes (Brezgin et al. 1977) and chemiluminescent son es 
(Konkov and Perov, 1976; Perov and Khrgian, 1980; Perov and Tishin, 1985) have been '";™ duce ^ A 
solar photometer (Subbaraya and Lai, 1981) and another optical sonde (Somayajulu et al., 1981) have been 
developed in India. However, the longest data record has been collected in the United States by the ROC 
and ROCOZ-A systems, which are described here. 

The ROCOZ-A ozonesonde is a four-filter, sequential-sampling, ultraviolet radiometer. 1 he instru- 
ment was originally developed for stratospheric soundings aboard an ARCAS rocket (Kreuger and McBride, 
1968a,b). This instrument was subsequently modified for launch aboard a smaller diameter Super-Lok. 
launch vehicle. In 1982, an instrument improvement program was initiated at NASA s Goddard Space 
Flight Center/Wallops Flight Facility (WFF). Changes were made to the center waydengths and spectral 
shapes of the ultraviolet filters. An integrated calibration facility was established (Holland et al., 1985), 
and new data reduction procedures were designed (Barnes et al., 1986). A description ot the present design 
of the radiometer has been published by Barnes and Simeth (1986). Because of the short data record of the 
present instrument, these data are not directly applicable to trend studies at this time. 

The ROCOZ-A ozonesonde is propelled aloft by a Super-Loki booster rocket. At rocket burnout, the 
instrument and its carrier coast to a nominal apogee of 70 km, where the payload is ejected for deployment 
on a parachute. The radiometer measures the solar ultraviolet irradiance over its filter wavelengths as it 
descends through the atmosphere. The amount of ozone in the path between the radiometer and the sun is 
then calculated from the attenuation of solar flux as the instrument falls. In addition, radar from the launch 
site measures the height of the payload throughout its descent. Finally, knowledge of the solar zenith angle 
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allows calculation of the fundamental ozone value measured by the radiometer; that is, ozone column 
amount as a function of geometric altitude. Ozone number density is the derivative of ozone column amount 
versus altitude. Combined with auxiliary atmospheric soundings for pressure and temperature, ROCOZ-A 
can duplicate the fundamental ozone values of backscatter ultraviolet (Barnes, 1988) and solar occultation 
measurements (Cunnold et al., 1989) from satellites. ROCOZ-A comparisons have also been made with 
the ultraviolet spectrometer on the SME satellite (Barnes et al., 1987b). 

Details of the measurements of the precision of ROCOZ-A ozone column amounts and ozone number 
densities have been published (Holland et al., 1985; Barnes et al., 1986). In addition, there are published 
results from an equatorial ozone measurement campaign (Barnes et al., 1987a) that found very low 
variability in stratospheric ozone, temperature, and pressure. From the results of this campaign, it is 
possible to estimate the precision of the remaining measurements in the ROCOZ-A data set. The full set 
of precision estimates for ROCOZ-A is given in Barnes et al. (1989). For ozone number densities and 
mixing ratios, the profile-to-profile repeatability of ROCOZ-A measurements is estimated at 3 to 4%. 

The estimates of the accuracy of the ROCOZ-A profiles are also given in Barnes et al. (1989). For 
ROCOZ-A ozone amount, the accuracy estimates come from an internal, unpublished error analysis. The 
analysis is based on errors of the effective ozone absorption coefficients used to convert the radiometer 
readings into ozone profiles, plus the differences between the ozone values at altitudes where two ROCOZ- 
A channels give simultaneous readings (Barnes et al., 1986). A laboratory flight simulator, based on long 
pathlength photometry (DeMore and Patapoff. 1976; Torres and Bandy, 1978), has been constructed to 
measure the accuracy of ROCOZ-A ozone readings. Publication of a detailed error analysis will follow the 
conclusion ot experiments with the simulator. It will complete the primary characterization of the ROCOZ- 
A ozonesonde. ROCOZ-A ozone number densities and ozone mixing ratios are estimated to be accurate 
to 5 to 1% and 6 to 8%, respectively. 

The vertical resolution of ROCOZ-A ozone profiles is 4 km (Barnes et al., 1986). This resolution 
comes from the data reduction algorithm for the profiles, since measurements from the instrument are less 
than 100 meters apart during flight. ROCOZ-A ozone column amounts are smoothed before differentiation 
for ozone density. The vertical resolution of the smoothed profiles has been adjusted to roughly match the 
8-km resolution quoted for the two limb scanning instruments on the Solar Mesospheric Explorer (Rusch 
et al., 1984), and the I- to 5-km vertical resolution for SAGE II (Mauldin et al., 1985b). 


Satellite Observations 

The Stratospheric Aerosol and Gas Experiment (SAGE I and II) 

SAGE I and SAGE II are both satellite-borne multi-wavelength radiometers employing solar extinc- 
tion (occultation) techniques to measure stratospheric aerosols and gases. Ozone profiles are determined 
from measurements of attenuation ot solar radiation by ozone in the most intensely absorbing portion of 
the Chappuis band, at 600 nm. SAGE I was launched aboard the dedicated Applications Explorer Mission- 
2 in February 1979 and operated for 34 months until November 1981 , when the spacecraft electrical system 
failed. SAGE II was launched on the Earth Radiation Budget Satellite in October 1984, and has operated 
continuously since then. Both SAGE 1 and SAGE II are in approximately 600 km circular orbits with 
inclination angles of 56 and 57 degrees, respectively, such that the latitudinal coverage is almost identical. 
Detailed descriptions of the instruments are given by McCormick et al. ( 1979), and Mauldin et al. ( 1985a, b). 
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The capabilities of the two instruments are not identical because of small differences in the instrument 
configurations and data processing algorithms. 

In the solar occultation technique, measurements are made of the solar radiation transmitted through 
the atmosphere as the sun sets (or rises) behind it relative to the spacecraft. The transmission for a given 
ray path is the ratio of the signal strength for that ray path to the signal when the sun is above the 
atmosphere. The vertical distribution of ozone is determined from the changes in the transmission. er 
wavelengths allow the concentrations of other gases and aerosols to be determined, the aerosol effects on 
the ozone transmittance may then be corrected. The retrieval algorithms are described in the OTP. Ch. 3, 
or in Chu and McCormick (1979), Mauldin and Chu (1982), and by Chu (1986). Because of the high s.gnal- 
to-noise ratio, a vertical resolution of 1 km is achieved in the stratosphere. 

It is important to note that the measurements performed by the SAGE instruments are self-calibrating, 
in that only relative radiance measurements are required to determine the transmission and, therefore, the 
distribution of atmospheric species such as ozone. Consequently, no absolute radiance calibration is 
necessary. The only requirement is that the instrument retain constant respons.vity for the duration of 
each spacecraft sunrise or sunset, usually about 100 seconds. However, the position of the line of sight 
must be accurately known, which requires very accurate data on the position of the spacecraft during 
observing events and a reasonably stable spacecraft. The OTP analysis indicates that absolute accuracy of 
the ozone values determined by SAGE 1 and SAGE II is about 6 to 9%, depending on aUdude^However 
the stability, or ability to detect changes, was 2 to 7% for SAGE 1 and 1.3 to 4 % for SAG . 
uncertainty in the difference between SAGE 1 and SAGE II ozone values allow changes of 2/, over the 
several years between their observing periods to be detected from 25 to 45 km. 

A limitation of the occultation technique is the relatively small amount of data obtained, two vertical 
profiles per orbit, and the changing location of the observed latitudes, rendering comparisons between 
occultation instruments difficult. The ozone trends detection uncertainty estimated above for the SAGE 
and II ozone data do not include errors caused by the undersampling of the natural variability of the ozone 
distribution in both the spatial and temporal domains. 

Ozone data from SAGE II between November 1984 to November 1988 have been archived at the 
National Space Science Data Center (NSSDC). The SAGE I ozone data have also been archived at NSSDC. 
Reprocessing of the SAGE I ozone data using the updated temperature correction information from Nt - 
NMC has been performed for the ozone trends study, and the data are currently bemg rearchived at 
NSSDC. The new temperature data introduce only small differences at high altitudes in the low-latitude 

regions. 


Solar Backscatter Ultraviolet Spectrometer (SBUV) 


The SBUV is a downward viewing double monochromator that was launched on the Nimbus-7 
spacecraft in 1978 to measure the UV albedo of the atmosphere and surface for the purpose of determining 
the vertical profile of ozone, in addition to total ozone. Singer and Wentworth (1957) suggested that 
observations from above the atmosphere, in which the fraction of sunlight reflected back to space (the 
planetary albedo) is measured as a function of wavelength, could be used to deduce the concentration o 
ozone as a function of pressure (i.e„ altitude). Other experiments utilizing the same principle have flown 
on Kosmos-65 , OGO-4, Nimbus-4, Atmospheric Explorer D, and most recently, NO A A-9 and the Japanese 

OHZORA satellite. 
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The SBU V made measurements in 12 wavelength channels, located at (in nm) 255.5, 273.5, 283.0, 
287.6, 292.2, 297.5, 301.9, 305.8, 312.5, 317.5, 331.2, and 339.8. Measurements in the channel at 255.5 nm 
were not used because they were contaminated by NO fluorescence; the next seven were used for extracting 
profile information, while the last tour (which are common with TOMS) are used to determine the total 
ozone. The solutions are given in terms of ozone amounts within the Umkehr layers. The analysis presented 
in the OTP shows that only data in layers 6 to 9 (approximately 28 to 48 km) are independent and reliable 
for trend studies. The vertical resolution of the solutions is 8 to 10 km. 


Again, in principle, this technique requires only the measurement of the ratio of the backscattered 
Earth radiance to the solar irradiance, with the attendant insensitivity to calibration accuracy and stability, 
but the use of a diffuser plate involves another optical element in the measurement of the solar irradiance 
which eliminates this advantage. On Nimbus-7, SBUV and TOMS shared the same diffuser plate. SBUV 
observations of the solar irradiance also showed that the diffuser plate was degrading with time, most 
rapidly at the shortest wavelengths. The diffuser model of Cebula et al. (1988) was developed from the 
SBUV observations and used by the OPT for the SBUV data reduction also. The archived results showed 
very large (25%) downward ozone changes in the upper stratosphere over 8 years, which was the reason 
for the original appointment of the Ozone Trends Panel. In this case also, it was found that the uncertainties 
in the diffuser model were much larger than expected, and that the changes in the stratosphere were 
probably smaller than indicated by the archived data and could even be slightly positive. In this case, the 
D pair again supports the assertion that the degradation is larger than predicted by the model, but there 
has been no method developed to date for finding an independent determination of the diffuser degradation 
at the shorter wavelengths. Thus, at the moment the SBUV data can provide no independent information 
on long-term trends in the ozone profile (see OTP). 

Limb Infrared Monitor of the Stratosphere 

The Limb Infrared Monitor of the Stratosphere (LIMS) is a six-channel limb scanning infrared 
radiometer that also was launched on Nimbus-7 in 1978 (Gille and Russell, 1984; OTP). This type of 
instrument measures the infrared radiation emitted by atmospheric molecules as the instrument scans 
across the limb. Because of the geometry, this technique has an inherently high vertical resolution and the 
ability to sound to high altitudes. Because it measures infrared emission, it can obtain measurements at all 
local times, resulting in very dense coverage. LIMS obtained ozone data from 15 to 64 km altitude, with a 
vertical resolution of 2.5 km. The absolute accuracy was about 10%, and the precision a few percent. 

Because of the small signals involved, it is necessary to cool the detectors. LIMS life was limited by 
the technology of that period, which dictated the use of a solid cryogen cooler. This resulted in a 7-month 
lifetime. Otherwise, it should be useful for trend measurements, especially since it provides a 1978/79 
determination. LIMS data were used in the OTP to help evaluate other data. 

2.1 .3 Measurements for Future Trend Analyses 


Many of the methods described previously will continue to be used to obtain data for trend studies. 
Where these methods are reasonably accurate and stable with time, this is essential, because it extends 
the length of the data record. However, as technology advances, improvements to old techniques are 
continually being introduced and new measuring systems are being developed. A number of these are 
shown in the timelines on Figure 2.1-3, and are described below. 
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Figure 2 


.1-3. Time lines for ozone measuring systems that are projected to be operating over the next decade. 


2.1 .3.1 Total Ozone Measurements 
Ground-based Measurements 

The Brewer Spectrophotometer 

The Atmospheric Environment Service of Canada developed the fully automated Brewer ozone 
spectrophotometer in the early 1980s. Development of the Brewer instrument was based, in part, on earlier 
work by Wardle et al. (1963) and Brewer ( 1973) and was carried out with the goal to replace or supplement 
the Dobson instrument in the world network. There are presently about 30 instruments operatmg m 15 
countries, including the seven Canadian stations (Kerr et al., 1988). At present, five stations are reporting 
data to the WODC. The record of routine measurements of total ozone made with the Brewer spectropho- 
tometer started in 1 982 . These records are presently not of sufficient length to carry out a proper independent 

trend analysis. 

The spectrophotometer is a modified Ebert type with a 1,800-lines/mm holographic grating used in 
the second order. It simultaneously measures the intensity of light at 5 wavelengths in the ultraviolet 
absorption spectrum of ozone with a resolution of 0.6 nm. The operational wavelengths are 306.3 nm, 3 10.0 

nm 313 5 nm, 316.8 nm,and 320.0 nm. Measurements at these wavelengths allow correction for the effects 

of sulfur dioxide, a potential interferent for Dobson ozone measurements in polluted air. A more complete 
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description of the automated instrument is given by Kerr et al. (1985). The Brewer instrument has been 
calibrated on an absolute scale using the ozone absorption coefficients of Bass and Paur (1985). Intercom- 
parison with the Dobson instrument has shown that the Dobson AD direct sun total ozone measurement 
is about 3% larger than the Brewer measurement. When this bias is removed, the long-term agreement 
between field Brewer and Dobson instruments at Toronto and Edmonton has been within 1% (Kerr et al 
1988; 1989a). 

Although the method to measure total ozone with the Brewer instrument is similar to that for the 
Dobson instrument, there are differences in sampling and reporting daily data. These differences arise 
primarily because the Brewer instrument is fully automatic and thereby capable of sampling continuously 
throughout the day. Measurements are screened and only those of good quality are used to determine the 
daily mean total ozone value that is reported. In general, good quality direct sun measurements can be 
made on about 75% of the days at a typical mid-latitude station. 

High Resolution Visible/Ultraviolet Absorption Spectroscopy 

Absorption spectroscopy affords a sensitive means for monitoring the column abundance of several 
stratospheric species (e.g., Noxon et al, 1979; McKenzie and Johnston, 1982; Pommereau et al., 1988a, 
1988b). A useful approach is to measure the absorption of light scattered by the sunlit zenith sky. An 
important feature of such measurements is that they can be carried out on cloudy as well as clear days. 
The optical paths (or air mass factors) relevant to a species in a stratospheric layer become very large at 
solar zenith angles greater than about 88 degrees, maximizing the slant-column abundance of absorber and 
hence the observed percent absorption. 

Modern optics, electronics, computers, and detectors allow atmospheric absorption spectroscopy to 
be done with very high precision, typically 0.05% absorption. Multichannel detectors eliminate the uncer- 
tainties that scanning optics introduce. Least-squares fitting of the observed spectra to those recorded in 
the laboratory with the same instrument for several atmospheric gases allows the abundance of more than 
one atmospheric species to be determined simultaneously. Under favorable spectral conditions, the sen- 
sitivity stated above can be enhanced a factor of 50 by such fitting. The stratospheric species that are most 
amenable to observing with ultraviolet/visible absorption spectroscopy have proven to be O,, NO,, NO,, 
and OCIO. The characteristics of one instrument that has been applied to these species is used here as an 
example of the type that is achievable today. 

The example (Mount et al, 1987) is a crossed Czerny-Turner spectrometer, with a grating chosen so 
that the appropriate color filters provide second-order light in the red region of the spectrum, approximately 
605—685 nm (for NO, observations) and third-order light in the blue region, approximately 400 to 450 nm 
(forO„ NO : , and OCIO observations). The instrument is used at approximately 0.5-nm spectral resolution. 
The detection system is a Reticon diode array cooled by a refrigerator to about -70°C; the array contains 
1 ,024 independent silicon diodes that simultaneously measure the spectrum over the indicated wavelengths. 
Eleven diodes cover the full width at half maximum of the instrumental spectral profile. Thus, spectral 
lines are highly oversampled. 

An important advantage of simultaneous measurement of the desired spectral interval (as opposed to 
scanning) is the elimination of time-dependent changes during the course of a measurement (particularly 
atmospheric scintillation effects). These can be quite important when the measured species absorb only a 
few tenths of a percent of the incoming light, as is the case for ozone and OCIO in general and for NO; at 
small zenith angles. After least squares fitting the observed sky data to the reference absorption spectra 
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measured in the laboratory, the one-sigma standard errors of the fitted column abundances are typically 
2-5% for large solar zenith angles for NO : . Ozone columns can be determined to 10-15% accuracy using 
the Chappuis bands. Future development is expected to improve this considerably. While most efforts 
with this technique have been on the measurement of other gases, the technology presumably could be 
applied to ozone measurements as well. 

The largest source of systematic error in the measurements is the evaluation of the air mass factors 
for scattered light that are used to convert measured slant column abundances to vertical column abun- 
dances. These errors largely affect the accuracy but not the precision of the data. Estimated error limits in 
air mass factors are no more than 10% for solar zenith angles near 90 degrees, and they are substantially 
less for smaller angles. 


Satellite Observations 

Future Flights of TOMS Instruments 

TOMS data have proven of great value in observing and studying the polar ozone depletions, as well 
as in providing information on the global decrease of total ozone. With the development of independent 
techniques for determining the diffuser reflectivity, TOMS data should be able to provide a determination 
of total ozone decrease over all areas of the globe that is independent of the ground-based Dobson network, 
although it should agree with the latter in areas where both obtain measurements. For these reasons, it is 
extremely important to continue TOMS-type measurements. The TOMS instrument on Nimbus-7 is now 
over 10 years old, and efforts have been made for future flights of TOMS-type instruments on several 

spacecraft. 

First flight of a new, technically upgraded instrument is scheduled for the fall of 1991 on board the 
Soviet Union’s Meteor-3 spacecraft. Unlike the sun-synchronous orbit of Nimbus-7, the Meteor-3 oibit 
will precess with a period of 225 days. Thus, the coverage will not be as uniform as that obtained with 
previous data, and adjustments will have to be made for diurnal variations. TOMS has also been selected 
to fly on a small U.S. Explorer-class satellite in 1993, and on the Japanese ADEOS satellite, scheduled for 
launch in February 1995. In these cases the orbit will again be sun-synchronous. 

An upgraded TOMS will also be part of the Global Ozone Monitoring Radiometer (GOMR), an 
operational instrument to make long-term ozone measurements. The GOMR, discussed further below, will 
be included in the complement of operational instruments on the NOAA afternoon free-flier. beginning in 
about 1997. 

2.1 .3.2 Profile Measurements 

Ground-based Measurements 

Brewer Umkehr 

Although Umkehr observations have been made with the Brewer Ozone Spectrophotometer for 
several years (Mateer et al., 1985, Kerr et al„ 1989b), this has been a developmental period during which 
the observational technique and the algorithm have changed. The present body of data is insufficient for 
trend analysis. 
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A Brewer Umkehr observation consists of zenith sky intensity ratio measurements for three wave- 
length pairs over the solar zenith angle range from 60 to 90 degrees. The profile retrieval data cover the 
altitude range of layers I to 10 (0-53 km). The averaging kernels for the present algorithm show that the 
retrievals for layers 4 to 8 ( 19-43 km) and marginally layer 9 (43-48 km) are suitable for trend estimation. 
Until a complete error analysis has been carried out for the Brewer Umkehr system, those for the standard 
Dobson Umkehr may be considered as only approximately applicable. Results of an intercomparison of 
Brewer Umkehr, Dobson Umkehr, and ozonesonde data have indicated that the quality of the Brewer 
Umkehr technique is comparable to that of the Dobson "short” Umkehr method (McElroy et al., 1989). 

Lidar 

The lidar measurements of the ozone vertical distribution are based on the Differential Absorption 
Laser technique (DIAL), which requires the simultaneous emission of two laser wavelengths characterized 
by a different ozone absorption cross section. Part of the laser radiation is scattered back to the surface, 
where it is collected by a telescope, detected by a photomultiplier and time-sampled to retrieve the altitude 
resolution of the measurement. The derivative of the logarithm of each signal is computed and the ozone 
number density is obtained from the difference of the derivatives, divided by the differential ozone cross 
section. The lidar technique has been described by Megie et al. (1977), Uchino et al. (1978), Pelon and 
Megie (1982), Werner et al. (1983), Pelon et al. (1986), Ancellet et al. (1988), Godin et al. (1989), and 
McDermid and Godin (1989). 

In principle the technique is self-calibrating, and thus particularly suited for the detection of trends. 
Two types of error affect the ozone lidar measurement. The random or statistical error is related to the 
signal-to-noise ratio. The systematic error is related on one hand to errors in the value of the ozone 
absorption cross sections used and their dependence on temperature, and on the other to uncertainties in 
the different physical processes, such as the Rayleigh and Mie scattering or the absorption by other species 
(N0 2 , S0 2 ). The systematic error can be reduced by additional measurements of temperature and aerosol 
vertical profiles or use of empirical models. The random error which results from this correction is less 
than 3/c in the whole altitude range. The choice of the two wavelengths used in the experiment depends 
directly on the altitude range monitored. In the troposphere, the limiting factor is the systematic error, 
which becomes smaller as the difference in wavelength is reduced. A wavelength separation of less than 
10 nm is used, and the time necessary to obtain a tropospheric ozone profile is about 10 min for an altitude 
resolution superior to 1 km and an overall uncertainty less than 5%. In the stratosphere, the separation of 
the wavelengths needs to be larger (—50 nm) in order to limit the statistical error. Furthermore, the range 
solution has to be degraded in the upper range of the measurement (above 30 km) to account for the rapid 
decrease of the signal-to-noise ratio with altitude. With the present systems, obtaining an ozone profile in 
the 15-50 km altitude range requires an integration time of 2 to 3 hours for an altitude resolution of 0.5-8 
km and a corresponding uncertainty of 2% to 10%. 

Microwave Radiometry 

Microwave radiometers (MR) are being used increasingly to investigate the middle atmosphere. An 
MR measures the thermal emission from rotational transitions of the molecule under investigation. The 
frequency employed is typically below 300 GHz (1 mm wavelength) to reduce the effects of tropospheric 
opacity (e.g., clouds or aerosols) and local thermodynamic equilibrium (LTE) conditions apply to —80 km 
altitude. The sensor operates on the superheterodyne principle that provides a very high spectral resolution, 
allowing the determination of the exact shape of the emission line. Since the shape of the emission lines 
are dominated by pressure broadening to 80 km altitude, profiles can be retrieved from an exact measure- 
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ment of the emission line shape. At low altitudes (~ 15 km) the emission line widths become very large, 
and this, together with instrumental problems (base line and non-linearities) makes low-altitude retrievals 
very difficult. The altitude range over which high-quality ozone profiles can be determined from MR 
measurements is typically 15 to 70 km (Zommerfeld and Kunzi, 1989) with an altitude resolution of -10 
km. The accuracy claimed for this sensor is 1 part per million by volume, and the statistical error is of the 
order of 15%. The measuring time required for a complete profile is approximately 1 hour. Another 
instrument, developed by De la Noe et al. (1988) has obtained measurements of the ozone distribution with 
reported random errors from 1 to 5% at 40 km. Connor et al. (1987) present an error analysis for data 
obtained in Antarctica during 1986; they show measurement accuracies of 15-19% (depending on altitude), 
supported by intercomparisons with other observations. A detailed study performed recently by Bevilacqua 
and Olivero (1988) showed that the best vertical resolution attainable is 6 or 7 km. 


It should be noted that the accuracy can be improved by using observing sites at high altitudes, thus 
reducing the attenuation by tropospheric water vapor. For a mid-latitude site at 500 m above sea level 
(Bern, Switzerland), the observing statistics for favorable conditions are >80% of the time during tall, 
winter, and spring, and -50% of the time for summer. Model calculations (Kunzi and Rubin, 1988) indicate 
that for >90% of available observing time the sensor has to be at —3,000 m altitude for the tropics, whereas 
in polar regions the site can be at any altitude. 

MR instruments seem particularly well suited to perform long-term trend measurements of total 
stratospheric and mesospheric ozone and its vertical distribution, since an MR is easily calibrated with 
black body radiators and thereby has a response that can be made stable over long periods. Furthermore, 
the operation and data retrieval of an MR can be fully automated and require only little maintenance by 
skilled operators. However, at present, only little experience is available with MRs in ozone research and 
more intercomparisons and experience with sites in different climatic regions is needed. 

Satellite Observations 

Solar Backscatter Ultraviolet Radiometer— Version 2 (SBUV/2) 

NOAA is currently flying the Solar Backscatter Ultraviolet Radiometer— Version 2 (SBUV/2) as an 
operational instrument for measurement of both total ozone and vertical profiles. The first SBUV/2 was 
launched on NOAA-9 in 1984 and continues to operate in overlap with the instrument launched on NOAA- 
1 1 in 1988. This series of instruments is planned to be in orbit through the lifetime of the present afternoon 
series of NOAA polar orbiting satellites (ca 1996). These satellites sound the atmosphere. Unfortunately, 
an onboard check of the diffuser plate incorporated in the NOAA-9 SBUV/2 failed to function properly. 
At the same time, the processed ozone data from this instrument indicate a time-dependent bias with 
respect to the ground-based information as well as the Nimbus-7 SBUV data that is opposite in sign to our 
general understanding of diffuser plate degradation. No data from either instrument are now available. 
This time-dependent bias is still under examination; if the cause for this disparity is discovered, the data 
will be reprocessed accordingly. An attempt has been made to correct the onboard diffuser plate check on 
the SBUV/2 follow-on units. 

The failure to obtain data from the NOAA-9 SBUV/2 is especially damaging, since the first 2 years 
overlapped with the last 2 years of the SBUV. The SBUV/2 data could have been used to correct the SBUV 
diffuser model, allowing a reliable profile trend to be obtained over the 8-year life of SBUV . Alternatively, 
the differences between the first 2 years of SBUV and SBUV/2 could have been used to get a direct 
difference, as was done for data from SAGE I and SAGE II. 
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Upper Atmosphere Research Satellite 

The Upper Atmosphere Research Satellite (UARS) is scheduled for launch by the U.S. National 
Aeronautics and Space Administration (NASA) in late 1991. The instruments have been designed for 
research purposes and the study of dynamical, chemical, and physical processes in the stratosphere, 
mesosphere, and lower thermosphere. Although it is not designed to measure ozone trends, a few of the 
instruments may provide data useful for trend studies. 


Ideally, each of the ozone measurements will be sufficiently accurate that it can be compared with 
the existing body of previous data. However, because of the small trends and the often larger systematic 
biases between different types of instruments, the clearest indications are likely to occur when similar 
instruments are used. Some possibilities are: comparisons of ozone measurements among infrared limb- 
scanners such as UARS 1R limb scanners (CLAES and ISAMS) as compared to LIMS and to future EOS 
instruments (see following); comparisons between measurements made by the UARS Microwave Limb 
Sounder (MLS) and the future EOS MLS; comparisons between measurements made by the infrared 
occultation instruments (UARS HALOE, and the later 1LAS [see following]); and perhaps also comparisons 
among the visible occultation observations performed by the successive SAGE instruments. 

Infrared Limb Atmospheric Spectrometer (HAS) 

The 1LAS is under definition study by the Japanese Environmental Agency for possible launch on 
the ADEOS satellite, now scheduled for February 1995. It is a solar occultation spectrometer, baselined 
to use three monochromators in the middle infrared. In addition to the measurement of ozone (9-10 

micrometer bands), several other trace species and aerosols will be measured over the 10-60 km altitude 
range. 

Earth Observing System (EOS) 

The U.S. National Aeronautics and Space Administration (NASA), the European Space Agency 
(ESA), and the Japanese Space Agency (NASDA) are developing a co-ordinated international system of 
polar platforms for remote sensing of the atmosphere, oceans, land surface, biosphere, and solid Earth. 
The goals are to begin or continue a set ol baseline observations of key variables that indicate the state of 
the entire Earth system, and continue them for at least 15 years. At this time, instrument selection is still 
underway. 

NASA has selected instruments for the definition phase for NASA’s polar orbiting platforms (NPOP- 
I and -2). It includes, for NPOP-I, an infrared limb sounder and SAGE III occultation instruments that 
should provide observations of ozone with high vertical resolution, high accuracy, and precision from cloud 
tops to 80 km. The limb sounder will also be capable of high horizontal resolution. Launch is scheduled 
for late 1996. In addition, a SAGE III instrument is planned as an attached payload for the Space Station 
to provide profiles at low to mid-latitudes to complement the NPOP-I high latitude coverage. The Space 
Station is scheduled for launch in 1995. For NPOP-2, definition studies are being carried out for limb 
sounders operating in the infrared, sub-millimeter, and microwave portions of the spectrum, which will 
make additional measurements of ozone profiles. 

The Global Ozone Monitoring Radiometer (GOMR) is planned to be the next in the series of NOAA 
operational satellite ozone monitoring systems. It is to fly on the satellite series subsequent to the current 
Advanced TIROS-N with launch about the mid-1990s. The GOMR will supplant the SBUV/2. The GOMR is 
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currently envisioned as two subsystems: an ultraviolet nadir and side-scan sounder based on Nimbus- 
TOMS technology, as noted above, and an infrared limb sounder based on EOS limb sounding instruments. 
The nadir sounder will provide total ozone in the sunlit portion of the Earth while the limb sounder will 
provide stratospheric profiles of ozone and other minor constituents important to ozone photochemistry 
and climate The final configuration of the satellite ozone monitoring system is undergoing evaluation. 


Three instruments providing ozone measurements in the stratosphere have been selected by the 
European Space Agency for flight on the first ESA Polar Platform (EPOP-1). GOMOS (Global Ozone 
Monitoring by Occultation of Stars) is an instrument intended to monitor ozone by stellar occupation, 
observing the full UV/visible/near-infrared spectrum of stars using a double grating spectrometer. SL1A- 
MACHY (Scanning Imaging Absorption Spectrometer) is a combination of two spectrometers operating 
also in the UV/visible/near-infrared part of the spectrum to observe transmitted, reflected, and scattered 
light. It is designed to be used for both nadir and limb sounding plus solar occultation and will provide 
measurements of tropospheric ozone. MIPAS (Michelson Interferometer for Passive Atmospheric Sound- 
ing) is a limb sounder operating in the mid-infrared region of the spectrum. All these instruments are 
presently in Phase A studies. 


2.1.4 Comparisons of Errors/Capabilities 


2.1 .4.1 Intercomparisons Relevant to Present Data 


Total Ozone— Dobson-TOMS-SBUV 

TOMS values for total ozone are in good agreement with those from SBUV, being systematically 
about 1% higher for reasons that are not completely understood. Fleig et al. (1989b), Reinsel et al. (1988), 
Bojkov et al. (1988), and OTP examined the drift of SBUV relative to an ensemble of 41 Dobson stations 
and found that there was a roughly linear drift between SBUV and Dobson of -0.38% per year over 6 
years. The intercomparison of TOMS data with the network of Dobson stations was also discussed as a 
function of time in OTP. This shows that TOMS total ozone declined slowly relative to the 4 1 -station 
ensemble between 1979 and mid- 1982— about -0.25% per year— and declined rapidly thereafter, at a rate 
of about -0.53% per year through October 1986, the end of the period studied. The decline for SBUV was 
similar. These results are confirmed by a recent study performed by McPeters et al. (private communication, 
1989) using a network of 39 stations as shown in Figure 2. 1 -4a. 

While total ozone measured by SBUV/TOMS clearly declined relative to that measured by the Dobson 
network further tests were necessary to determine whether the problem is due to degradation of the 
satellite instruments or to changes in individual Dobson instruments. One such test involved satellite 
overpass of Mauna Loa, Hawaii, on those days during which the World Primary Standard Dobson instru- 
ment, instrument No. 83, was undergoing its regular recalibration at that location. Komhyr et al. (1989) 
report that the stability of this instrument has been maintained to 0.5-1% since 1962. The trend of TOMS 
relative to this instrument, also shown in Figure 2. l-4a, is very similar to that obtained relative to the 39- 
station ensemble— little change from launch through mid- 1982 followed by a rapid decline thereafter 
amounting to approximately -0.7% per year through 1988. 

As described above, a purely internal check for instrument degradation is possible for SBUV. The 
“D" wavelength pair, consisting of wavelengths 305.8-312.5 nm is sensitive to ozone because of its large 
ozone cross section and is very insensitive to wavelength-dependent degradation errors because of the 
small wavelength separation. The drift of normal SBUV ozone relative to D pair ozone exhibits the same 
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Figure 2.1 -4a. Comparison between archived TOMS total ozone data and Dobson data. These TOMS 
measurements have been corrected for diffuser plate degradation by the method described by Cebula Park 
and Heath ( 1 988) and the Ozone Processing Team ( 1 989) . T riangles show the difference between TOMS and 
the World Primary Dobson Instrument No. 83. (There are reasons to think the 1980 Dobson data are of lower 
quality than those for other years. ) Squares show differences between TOMS and a network of 39 Dobson stations. 
Both comparisons show a sharp drop of TOMS relative to Dobson ozone, especially since 1982. 


pattern as that observed for SBUV/TOMS relative to Dobson and for TOMS relative to instrument No 
83— little change before 1982 followed by rapid change thereafter. 

The conclusion that there was about 3.4% drift between SBUV/TOMS and the Dobson network 
between 1979 and 1987, with much of the drift occurring after 1982, seems inescapable. The fact that the 
same pattern of relative stability in the TOMS-Dobson comparisons before 1982, followed by a sharp 
decline between 1982 and 1986, was repeated for so many of the ensemble of 41 independent Dobson 
stations was strong evidence for SBUV/TOMS instrument degradation, and insufficient correction. The 
complete confirmation by comparison with the World Standard Dobson instrument and by internal evidence 
of relative pair drift furnish compelling additional evidence that this trend relative to Dobson is a robust 
result. As pointed out previously, this drift of SBUV/TOMS is believed to result from degradation of the 
diffuser plate used to measure solar flux, with incomplete correction. 

TOMS data for the ground station overpasses have been reprocessed using the “pair justification” 
calibration. The comparison of them to the data obtained by Dobson 83 during its calibration at Mauna 
Loa are shown by the triangles in Figure 2. 1 -4b. from which it can be seen that there is no large, systematic 
drift between them greater than a few tenths of a percent over 9 years. Because the techniques are so 
different, this agreement between the Dobson and TOMS support the conclusion that both are correct 
Further confirmation is provided by the squares in Figure 2. 1 -4b. which presents comparison data between 
TOMS and a network of 39 stations. They suggest a drift of 0.5% over 8 years, close to the level of 
calibration stability achievable for this technique (it could also be explained by about a 1 %/year increase 
in tropospheric ozone). 
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Fiaure 2 1 -4b Comparison between TOMS total ozone data, where diffuser plate degradation was corrected 

memod, and Dobson data. Triangles she J JSS.T"' 

Dobson Instrument No. 83 calibration measurements. (There are reasons to think the 1 980 Dobson data are of 
lower aualitv than those for other years.) The squares show the differences with the network of 39 stations used 
or Figure 2 1 4a. These results show that once corrected, TOMS has not drifted relative to instrument 83 between 
1979 and 1988, and has only changed by about 0.5% relative to the 39 station network between 1979-1987. 

This may be a real effect (see text). 


Vertical Profiles 


SAGE 1/ SAGE II 

The accuracy of estimating ozone trends using SAGE I and SAGE II ozone data is discussed in OTF\ 
Chapter 2. It is limited by the systematic errors between the two instrurnents. These errors range from 
5 8% below 25 km to 1 5% at 30 km. 1.6% at 40 km, and 3.0% at 50 km he.ght (see F.gure 2.1-3c). Tins 
assetsmenT doe s riot incibde any spatial or ienrpora, satnpiing error assorted with ihose two data sets. 

SAGE 1 1 SAGE II and SBUV 

An intercomparison of SAGE 1/11 and SBUV ozone measurements was performed for the periods 
February 1979 through November 1981 and October 1984 through December 1986. It is described in de ai 
in OTP. For this comparison, the SAGE I/II profiles were integrated, from concentration as a fu " ct ‘°" ° 
geometric altitude, the primary SAGE product, into Umkehr-layer amounts. In Umkehr layers 7 8, and 
9 SAGE I ozone layer amounts were consistently lower than SBUV in 1979-1981 by or A u 
U ozone layer amounts were consistently higher than SBUV in 1984-1986 by 10 ± % The results , 

Umkehr layer 6 consistently showed SAGE 1 about 4% higher than SBUV "J increased 

about 10% higher than SBUV in 1984-1986. The 1980-1985 offset between SAGE and SBUV increased 
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monotomcally with altitude between Umkehr layers 6 and 9. Coupled with the other information on 
degradation of the SBUV diffuser plate and the self-calibrating nature of the SAGE instruments, this 
comparison supports the conclusion that the large ozone changes reported by SBUV at high altitudes are 
due to larger diffuser degradation at the shorter wavelengths that allowed in the SBUV data reduction. 

SBUV, SAGE /, and LIMS Ozone Intercomparisons 


In OTP (Chapter 5) the ozone observations during March and April 1979 from three simultaneous 
satellite instruments, namely the SBUV, SAGE I, and LIMS, were examined. (For a discussion of LIMS 
see OTP or Gille and Russell, 1984 and Remsberg et a!., 1984.) The analysis is an intercomparison of the 
ozone data from these instruments in terms of ozone column abundance in four Umkehr layers (layers 6 
to 9). In the case of SAGE 1 and LIMS, it was necessary to convert their primary measured ozone quantities 
to ozone Umkehr layer amounts for comparison purposes. Since the sampling locations of the SAGE 
instrument are distributed in a narrow latitudinal belt for a 24-hour period, as opposed to the nearly global 
observations of the SBUV and LIMS for the same period, the comparison analysis is carried out by 
following the SAGE I latitudinal sampling progression. Because the three satellite instruments never 
sampled at the same locations and times, the comparison was made in terms of the zonally averaged ozone 
layer amount in these Umkehr layers. 

For the case of April 1979 and SAGE I sunset, the standard deviation of the SAGE I instrument about 
its zonal means is less than 10% except at high latitudes. The standard deviation for LIMS about its zonal 
mean is about 4%, and that for SBUV is about 2%. All three satellite systems show the same pattern with 
latitude. The results show that SBUV data are about 3% systematically lower than the mean of the three 
instruments in layer 6 but not systematically different from the other instruments in layers 7, 8, and 9. For 
each of the three satellites, for each of the four Umkehr layers, and for each of the four times, the root- 
mean-square deviation from the average of the three instruments is given in Table 2.1-1. For 15 out of 48 
entries in Table 2. 1-1, the deviations are less than 3%; for 23 out of 48 entries the deviations are between 
3 and 5%; and for 10 out of 48 entries the deviations are greater than 5% but less than 10%. This set of data 


Table 2.1-1. Estimated overall percentage differences of the calculated zonal mean ozone layer 

amount of SBUV, SAGE, and LIMS with respect to the average of these three instru- 
ments 


1 Sunrise (March 79) 

2 Sunrise (March 79) 

3 Sunrise (April 79) 

4 Sunrise (April 79) 


1 Sunrise (March 79) 

2 Sunrise (March 79) 

3 Sunrise (April 79) 

4 Sunrise (April 79) 


Layer 6 
SBUV SAGE 

4.5 6.3 

2.8 3.1 

5.1 3.7 

3.8 3.6 

Layer 8 
SBUV SAGE 

1.6 3.6 

3.0 4.1 

1.6 3.2 


UMKEHR LAYER 


LIMS 

SBUV 

2.6 

2.4 

1.6 

3.8 

1.3 

2.9 

2.1 

2.8 


LIMS 

SBUV 

4.4 

4.2 

4.2 

3.3 

3.9 

3.0 

3.3 

6.2 


Layer 7 

SAGE 

LIMS 

4.2 

4.6 

2.8 

4.0 

1.5 

3.7 

1.4 

3.4 

Layer 9 

SAGE 

LIMS 

6.1 

8.9 

7.7 

8.6 

7.3 

9.5 

1.8 

6.9 
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shows that three totally different satellite systems agree with each other to about 4% in measuring zonal 
mean ozone amounts in the four Umkehr layers in the middle to upper stratosphere in April, 1979. The 
conclusion is that this is the level of agreement that has been demonstrated between different measuring 
techniques as of this time. The differences are probably due in part to differences in systematic errors, 
such as the absorption coefficients, as well as to effects of the individual techniques. A further inference 
is that when one is looking for small changes, at this time it is probably necessary to use instruments of 
the same or closely related types, with the same systematic errors. 

SAGE I and SAGE II Comparisons with Ozonesondes 

As a further step in understanding data quality, electrochemical concentration cell (ECC) and Mast 
ozonesonde profiles from the World Ozone Data Center (WODC) data base have been compared with 
SAGE I and SAGE II ozone profiles, and discussed in OTP (Chapter 5). SAGE I and SAGE II data sets, 
composed of six profile pairs each, were chosen for comparison at the following six stations: Edmonton 
(53NJ 13W), Churchill (58N,94W), Hohenpeissenberg (47N,1 IE), Wallops Island (38N,75W), Lindenberg 
(52N,14E) and Praha (50N,14E). 

SAGE I and SAGE II measurements over these stations occurred at regular intervals as either of 
these instruments sampled the latitude band containing any particular station. The six profile pairs repre- 
senting the SAGE I period are shown in the top six panels of Figure 2.1-5. The WODC profiles were 
converted to ozone number density, the fundamental unit of ozone measured by both SAGE I and SAGE 
II. The ozone number density profiles are plotted as a function of pressure from approximately 250 to 5 
mb. The SAGE I profiles are plotted as the thin curve bounded by two dashed curves, the dashed curves 
representing la uncertainties on the SAGE I measurement. The dates of each measurement, the SAGE I 
event type (sunrise or sunset measurement), the spatial and temporal separation of the station-to-SAGE 
measurement tangent point, and the correction factor tor the ozonesonde are provided on each graph. 
Typically, SAGE averages over a cylinder of diameter 1 km and length 200 km. The ozonesonde correction 
factor represents the ratio between the total ozone as measured by the Dobson instrument at the station 
and the total ozone derived by integration of the ozonesonde profile. Correction factors deviating greatly 
from unity are indicative, generally, of ozonesonde profiles of questionable reliability. In the analysis done 
here, ozonesonde profiles were not selectively screened out by rejecting those with correction factors 
outside of the intervals mentioned above. 

The WODC and SAGE I profile pairs in Figure 2.1-5 had a maximum spatial separation of 481 km at 
Praha and a minimum separation of 71 km at Lindenberg. The maximum and minimum temporal separation 
between the profile pairs occurred at Hohenpeissenberg and at Lindenberg, respectively. Thus the profile 
pair at Lindenberg is the closest coincidence between the in situ ozone measurement and the remote sensing 
ozone measurement. This ozonesonde profile is almost entirely contained within the uncertainties of the 
SAGE I measurement. The precision of an ozonesonde measurement, which is not provided by WODC, 
has been estimated at 5-10% from ground to tropopause, and at 5% from tropopause to 25 km (Barnes et 
al., 1985; Hilsenrath et al., 1986). Similar results are shown on Figure 2. 1-6 for the SAGE II instrument. 

Estimates of the differences between the ozonesonde and satellite ozone measurements are given in 
the lower panels of Figures 2.1-5 and 2.1.6 over the altitude range from 15.5 to 30.5 km. Each difference 
profile was calculated by first interpolating the six ozonesonde profiles to the altitudes at which SAGE I 
and SAGE II measure ozone. For each altitude, the six ozonesonde values were averaged as were the six 
SAGE I or SAGE II values. The mean values were used to compute the percentage difference between 
the ozonesonde and the satellite measurements relative to the mean of the ozonesondes. The mean 
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Figure 2.1-5. Samples of profile pairs measured by ozonesondes and the SAGE I instrument appear in the 
upper six panels. The lower panel is the mean percentage difference between the ozonesondes and the SAGE 
I ozone. Ozonesonde profiles were interpolated to the altitudes 15.5, 16.5, ..., 30.5 km. The dashed curves 
represent one standard error. 
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Figure 2.1-6. Samples of profile pairs measured by ozonesondes and the SAGE II instrument appear in the 
upper six panels. The lower panel is the mean percentage difference between the ozonesondes and the SAGE 
II ozone. Ozonesonde profiles were interpolated to the altitudes 15.5, 16.5, ..., 30.5 km. The dashed curves 
represent one standard error. 
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percentage difference is plotted as the solid curve in the lower panels of Figures 2.1-5 and 2.1-6 with the 
dashed curves representing the standard error of the mean percentage difference. Using these relatively 
small samples, the differences between the satellite and the ozonesonde measurements are on the order of 
5-10%. 

A similar analysis to that described above was performed on a larger sample of ozonesonde and 
satellite pairs from the same set of stations listed above. An ozonesonde profile and a satellite profile were 
considered a pair if the satellite measurement location was within 1,000 km of the sonde launch station and 
within 5 days of the launch time. This criterion yielded 2,306 SAGE I-WODC pairs and 1,474 SAGE II- 
WODC profile pairs. These sets of pairs had a mean temporal separation of 2.4 ± 1.4 days and a mean 
spatial separation of 700 ± 200 km. The results of averaging, to form the percentage differences relative 
to the ozonesondes, are shown in Figure 2.1-7, with the horizontal bars representing the 95% confidence 
interval of the mean percentage difference. Between 20.5 and 27.5 km, the differences are mostly negative 
and centered at approximately -6%. Above 27.5 km, the differences computed using both SAGE I and 
SAGE II ozone are consistently tending toward higher values, a possible indication of ozonesonde pump 
inefficiency at high altitudes. Below 20.5 km, the SAGE II comparisons with the ozonesondes show an 
essentially constant difference with magnitude less than 5%, but with increasing uncertainties owing to the 
spatial and temporal ozone gradients below the ozone concentration peak. 

The conclusion is that comparisons of SAGE I and SAGE II profiles, colocated in time and space 
with ozonesonde profiles archived in the WODC data base, reveal that both the SAGE I and SAGE II 
satellite measurements consistently agree to within ±7% above approximately 20.5 km. 

SAGE-Umkehr 

The results of the trends analysis performed on the SAGE I/II data and Umkehr are discussed in 
Chapter 5 of OTP. At present, a direct comparison between the two data sets, as far as ozone trend is 



Figure2.1-7. Mean differences between SAGE I and ozonesondes (open squares) and SAGE II (solid squares) 
and ozonesondes. Error bars are the 95% confidence interval on the mean (±2<x). 
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concerned, has not been performed. The reprocessed SAGE I data that were used for the ozone trends 
study are being archived at NSSDC. Comparisons between Umkehr data and the previously archived 
version of the SAGE I ozone data showed some biases between the two data sets within the time period 
from 1979 to 1981 (Newchurch et al., 1987). An updated comparison with the revised SAGE I ozone data 
has not been performed. No comparison between the SAGE II ozone data and the Umkehr data has been 
performed. The changes derived from SAGE I-SAGE II differences and the Umkehr data are discussed in 

Section 2.3. 


2.1 .4.2 Intercomparisons Relevant to Future Trend Analyses 

Total Ozone— Dobson Versus Brewer and TOMS 

The Brewer data have been used in a number of intercomparisons which are of importance to present 
trend analyses and will be relevant to future trend studies. A number of sites, including five stations in 
Canada, have been making ground-based total ozone measurements on a routine basis with both the Brewer 
and Dobson instruments. These data, plus the TOMS measurements, provide three independent sets of 
routine field measurements of the same atmospheric constituent over a time period of a few years. 
Intercomparison between these data sets has been useful in evaluating monitoring performance and capa- 
bilities. Results of comparisons between Brewer and Dobson data at Toronto and Edmonton have shown 
an agreement within 1% and a statistically insignificant long-term drift (Kerr et al., 1989a; 1989b). Com- 
parison between TOMS data and Brewer data has shown that the TOMS total ozone has been decreasing 
relative to that of the Brewer at a rate of 0.6% per year at Toronto between 1983 and 1987, and 0.61% per 
year at Edmonton between 1984 and 1987 (Evans and Kerr, 1989). These values are similar to the 0.64% 
per year decrease of TOMS data relative to the ensemble of Dobson data for the same time period (Figure 
2.1-4). 


Vertical Profiles 

Micro wave-Umkehr 

Only one long-term (greater than 1 year) intercomparison of ozone profiles obtained from Umkehr 
observations and microwave measurements has been published (Lobsiger, 1987). From April 1984 to April 
1985, 334 daytime ozone profiles were obtained with a microwave sensor operated at the University of 
Bern, Switzerland. For Umkehr layers 6 through 9, simultaneous data were available for 211 days from 
the Arosa Umkehr station at a distance of approximately 200 km. The seasonal variations observed by 
both instruments are very similar and the corresponding correlation coefficients for layers 6-8 are 0.87, 
0.87, and 0.74, respectively. 

SBUV — Microwave 

A comparison was carried out between ozone profiles derived from microwave data obtained at the 
University of Bern, Switzerland and SBUV data. For the latter, archived data (i.e., with no allowance for 
SBUV instrument drift) were mapped using the Kalman filter technique and evaluated at Bern. Daily data 
were calculated at pressure levels between 30 and 0.7 mb; when plotted as time series they showed regular 
variations with season. The Bern data were available for four periods between December 1981 and June 
1987. All data were smoothed with a gaussian filter 10-days wide, to reduce short-term fluctuations. 


189 


GLOBAL TRENDS 


In all periods, there were appreciable biases. Averaged over altitude, the mean differences ranged 
from 30% to - 15%. The level showing the least bias was at 7 mb, where differences ranged from - 3% to 
- 13%. The standard deviations around these biases between 10 and I mb are 7-10% in the last three 
periods, with correlations of about 0.5. Because there were changes in the instrument, observational 
method, and data reduction, it is likely that the changes between the first and second periods, and between 
the third and fourth, are due to instrumental effects. (During these observational periods, differences 
between SBUV and the Arosa Umkehr observations showed differences of 0 to -6%, with standard 
deviations of about 7% and correlations of about 0.5.) The results are not appreciably changed if the 
comparisons are carried out for ozone amounts integrated over Umkehr layers. The shapes of the vertical 
profiles were also compared. There were occasions when the agreement was reasonably good, but other 
situations when there were large differences. There was no obvious correlation of these differences with 
season. 

In conclusion, these microwave observations, while somewhat noisy, and the SBUV observations 
appear to follow ozone variations and obtain roughly similar results for the day-to-day changes. However, 
there are large biases which are not entirely understood, and which appear to be rather sensitive to the 
details of the observations. In order to look for trends, the observing and data reduction will need to be 
standardized and held constant. As absorption by tropospheric water vapor is a major interference, 
observing from high, dry sites should reduce some of the observational problems, as noted previously. 

De la Noe et al. (1988) reported an intercomparison between microwave ozone profiles and profiles 
obtained by the BUV instrument on the OHZORA satellite during September 1985. Over the altitude range 
37-55 km, they found an average rms difference of less than 10% for distances between satellite and 
microwave soundings below 1,500 km, and in general less than 20% for larger distances. They also report 
comparisons with SAGE II, and find average rms differences of 18% during 1985. Connor et al. (1987) 
report approximately 10% agreement in microwave-ozonosonde comparisons. 

LIDAR — Umkehr — Satellite 


Recently, the capability to make lidar measurements of the ozone number density has been demon- 
strated in France, the Federal Republic of Germany, Japan, and the United States. Two systems are 
presently operational for systematic stratospheric measurements. The first one is located at the Observatoire 
de Haute-Provence (OHP), France and has been operating since October 1986. The second one, located 
at Table Mountain Observatory (TMO), California, USA began systematic measurements in February 
1988. Comparisons have been performed with other instruments such as the Umkehr, SBUV, and SAGE 
II. Concerning the short Umkehr method, a comparison based on 43 profiles at OHP indicated that lidar 
was about 3% higher than Umkehr in layers 5 and 6, but with the Umkehr values larger than the lidar 
values by about 7% in layer 7. Concerning SBUV, the comparisons have involved 12 measurements 
performed in July 1985 and 1986. The lidar values are generally higher than SBUV, with an average 
difference of about 10%. The comparison with SAGE is more direct in that it does not require conversion 
to mixing ratios and pressure coordinates, but there are fewer matchups. In general, the difference is less 
than 15% for each profile. Finally, a comparison of the ozone seasonal variation at different altitudes as 
determined by the lidar, SBUV, and Umkehr instruments has also been performed at OHP and shows that 
the shapes are very similar, but that the lidar data show a significantly larger amplitude (Godin et al., 1989). 

An intercomparison campaign involving the TMO lidar, the mobile lidar station recently developed 
by the Goddard Space Flight Center, SAGE II. rocket- and balloon-borne instruments (i. e., ROCOZ-A 
and ECC sondes) was organized at Table Mountain Observatory in October-November 1988. The corn- 
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parison between the two lidars, which could not be operated at the same time because of interferences 
between the two systems (located 600 meters apart), shows a day-to-day agreement within 10% in the 20 
km-40 km altitude range. This is somewhat greater than the error bars of the two systems, but it can be 
explained by the fact that the instruments didn’t observe exactly the same air masses at a period when the 
day-to-day variability of ozone was quite high. The comparison with the other instruments involved in the 
campaign is based on 3 days of observation. The results show an agreement within the instruments’ error 
bars in most cases. However, it is clear that the agreement depends strongly on local meteorological 
conditions, because the different measurements are not simultaneous. In particular, the agreement was 
less satisfactory during the passage of an airmass ot tropical origin characterized by a lower ozone content. 

2.1 .5 Temporal Sampling Requirements for Total Ozone Trend Detection 

Total ozone measurements are made on a routine basis at many stations around the world. The data 
are then reported to the WMO World Ozone Data Center in Toronto, Canada and printed regularly in a 
series of publications entitled Ozone Data for the World. Ideally a station takes several readings each day, 
with the times at which these readings are taken centered around local noon. These measurements are 
combined to form daily means which are then reported to Ozone Data for the World for publication. The 
daily means are averaged to form monthly means, and it is the monthly means that have been used in all 
statistical trend analyses. In practice, however, problems such as local weather, funding, and operator 
availability ensure that total ozone readings are not made every day. Understanding the effects of such 
missing data on the calculated trends is important in assessing the results of trend analyses and in designing 
future monitoring networks. A recent study by Tiao et al. (1989) examined the consequences of a reduction 
in the daily sampling frequency on trend estimates using data from Ozone Data for the World. The stations 
with the most complete daily records were chosen and year-round trends were estimated from data sets 
where the monthly total ozone means were calculated by using first all the data (1/1), then every other daily 
value (1/2), next one in four daily values (1/4), and so on down to one in sixteen daily values (1/16): in the 
last case just two daily values were thus used to calculate each monthly mean. The results are shown m 
Table 2.1.2. The trend estimates and their standard errors do not vary very much as the sampling rate is 
reduced to the 1/4 rate, although some differences are found at Mauna Loa and Hohenpeissenberg. These 
results seem to indicate that for the period 1970-1985 the precision of the long-term trend estimates will 
not be appreciably different at stations that sampled every other day than they would have been if that 
station had sampled every day. The reason for this insensitivity is the high day-to-day autocorrelation 
found in total ozone readings: the larger the day-to-day autocorrelation, the less sensitive the monthly 
mean (and hence the trend estimate) is to missing daily values. It should be noted that Tiao et al. eliminated 
the daily values systematically and spaced the remaining daily observations evenly through the month. 
These results are in agreement with the observation that for any particular calendar month, at a particular 
station, the ratio of the standard deviation of the daily total ozone measurements within the month to the 


Table 2.1-2 Summary of Ozone Trend Estimates (in % per year) in total ozone and associated 
standard error 


Sampling rate 

1/1 

1/2 

1/4 

1/8 

1/16 


Tateno 

- .002 ± .082 
- .022 ± .082 
-.010 ± .080 

- .046 ± .073 

- .009 ± .082 


Hohenpeissenberg 

- .078 ± .061 

- .099 ± .066 
-.136 ± .080 
-.103 ± .083 
-.187 ± .084 


Mauna Loa 

-.382 ± .090 
-.370 ± .091 

- .300 ± .095 

- .348 ± .092 
-.322 ± .093 
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standard deviation of the monthly means from year to year is typically found to be about 2, indicating that 
there are about 4 truly independent observations each month. 

Tiao et al. (1989) addressed a related issue which concerns estimating relationship between variables 
in two neighboring locations. Employing total ozone data from two European stations (Arosa and Hoh- 
enpeissenberg) and those from three Japanese stations (Sapporo, Tateno, and Kagoshima), they have found 
that the correlation estimates based on monthly averages of daily measurements are insensitive to temporal 
sampling rate when the measurements are taken on the same days across the locations. On the other hand, 
the estimates deteriorate as the time lag between measurements taken at the two locations is increased. 
The main reason for this phenomenon is likely to be that relations on a monthly scale are primarily due to 
short-term contemporaneous correlations on a daily basis. Although their empirical findings are based on 
the same variable (total ozone) at two different locations, similar phenomenon is likely to occur between 
different atmospheric parameters at the same location. Thus, in designing a measurement program when 
relations between variables are of interest, it is important to exercise care in arranging the timings of 
different measurements both across and within stations. 

Tiao et al. ( 1989) also investigated how different lengths of record and differing degrees of month-to- 
month autocorrelation within the data record affect the trend estimate and its standard error, and concluded 
that the precision of the trend estimates depends critically on the degree of the month-to-month autocor- 
relation in the measurements. If the month-to-month autocorrelation is large, it is indicative of unexplained 
medium to long-term variations in the data, and the standard error of the trend estimate is consequently 
large. It will be most difficult to provide early warning of a trend in total ozone in an area where the month- 
to-month autocorrelation is naturally high. In practice this means the tropics, and unfortunately there are 
not many long-term series of measurements in this region. 


2.1.6 Concluding Remarks 

In the previous sections, numerous past and emerging techniques for the measurement of ozone have 
been described. Where possible their accuracies and precisions have been presented. A few further remarks 
on the application of these techniques to the detection of trends are in order. 

Following the definitions given in the introduction, if the systematic errors are truly consistent (as, 
for example, if they are dependent on a poorly known absorption coefficient), the changes or trends that 
can be detected may be smaller than the absolute error of a system’s measurements. This is often the case; 
applicability then depends on how stable the system or instrument remains over time. 

Ground-based instruments can be checked and recalibrated as necessary. These are generally regarded 
as being capable of long-term stable operation. This is technically true; actual results may depend on 
funding, the importance attached to the measurements, the skill of the observer, and other non-technical 
factors. An advantage is that the same instrument is used repeatedly for the observations, so instrument- 
to-instrument variations do not come into play at a single location. 

Balloon and rocket measurements may be calibrated immediately before flight, but only rarely 
afterward. In addition, the conditions under which observations are made are generally different from those 
in which the calibrations are carried out. Perhaps most seriously, small instrument-to-instrument variations 
add an additional source of noise to attempts to look for trends. Because there is usually a relatively limited 
number of such observations, they have been of limited use for trend studies to date. 
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The features of satellite instruments are in some ways the converse of ground-based instruments. 
They cannot be checked in orbit, although they usually have a capability for in-flight calibration of at least 
some of the critical features of their designs. However, again it is the same instrument used at all locations 
as well as all times, so gradients and changes may be detected, depending on instrument stability. There 
is probably not enough experience to say much about instrument-to-instrument differences; the experience 
with SAGE I and II is encouraging that these effects can be small, but the SBUV-SBUV/2 problems are a 
cause for concern. Comparisons between different techniques are often illuminating. In principle the 
differences could be large, but the differences less than ~4% found among SAGE, SBUV, and LIMS in 
1979 is an indication of the present level of agreement between different techniques that can be achieved. 
As noted above, these differences are probably due in part to differences in systematic errors, implying 
that at the present time, when one is looking for small changes, it is necessary to use very similar if not 
identical techniques, including instruments and data reduction algorithms, to keep the systematic errors 
constant. In future it may be possible, with a great deal of care, to reduce the systematic errors to the point 
that data from different instruments can be combined in trend analyses. This will require further laboratory 
measurements of the ozone absorption coefficients and their temperature and pressure dependences, 
comparisons of the retrieval algorithms, and field intercomparisons under as many conditions as possible. 


A large amount of data has been collected over many years by instruments whose characteristics are 
known reasonably well. Initially, few of these measurements were made for the purpose of determining 
long-term changes. Fortunately, some of these data have the requisite accuracy and stability to provide a 
reasonable data base for efforts to detect and quantitatively measure changes in the ozone overburden and 
its vertical distribution, such as those reported by the OTP and later in this chapter. 


Looking to the future, the measurement situation for total ozone is moderately satisfactory. The 
present network of Dobson stations is continuing for the most part, and provides a reasonable network for 
detecting regional and probably zonal trends in the middle latitudes of the Northern Hemisphere. New 
ground-based instrumentation appears to provide compatible data, allowing a continuation of the measure- 

ment series. 


For global observations, the TOMS instrument continues to operate, having already provided almost 
one solar cycle (11 years) of data. A new technique for internal calibration appears to allow the TOMS to 
provide an entirely separate data record for total ozone that can be compared and checked against the 
ground-based network. Additional TOMS instruments are planned for flight on a variety of spacecraft 
during the 1990s, eventually reaching operational status late in the century. In the interim, if the problems 
with the instrument and data reduction are solved, the operational SBUV/2 could potentially provide 
another very valuable data record from the mid-1980s until the operational TOMS is launched. 

The situation for the measurement of trends in vertical profiles is somewhat less satisfactory. In the 
past the ground-based Umkehr techniques were subject to interference by aerosols. It is believed that these 
may now be corrected, using new techniques and independent measurements of aerosol distributions, but 
further testing and validation is necessary. The new lidar and microwave methods also require further work 
to fully understand and characterize them. In addition there are sampling problems; the number of stations 
taking regular Umkehr observations is limited, and most of them are in the mid-latitudes of the Northern 
Hemisphere. At present, there are very few stations employing the new technology. 

For the present, only data from the SAGE instruments appear to have the long-term stability to allow 
trend detection. SAGE II (and its spacecraft) are already over 5 years old, and there is not a planned 
successor until 1996. A failure could result in a lengthy gap. In addition, these data are subject to the 
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limitations of the sampling of the occultation technique and the attendant uncertainties in comparing 
different years (see OTP, Chapter 5). Again, if the problems with the SBUV/2 can be successfully resolved, 
they could potentially make an extremely valuable contribution. Other profile measurements will be 
obtained by research instruments on UARS, but it is unclear how well they will be able to extend the 
previous data records. An infrared limb scanner is scheduled to become available in 1996. 

Present planning for ozone trend measurements suggests a systems approach, based on complemen- 
tary ground-based and satellite observations. In this, satellite measurements can be checked for unexpected 
drifts through comparison with a network of carefully calibrated, maintained, and operated ground-based 
instruments. In turn, the satellite instruments act as transfer standards, indicating whether any of the 
ground instruments are drifting relative to the others. This will be especially challenging with respect to 
the measurement of vertical profiles. The stations in the Network for the Detection of Stratospheric Change 
(NDSC), which will include new instruments (lidar, UV, 1R, and microwave spectrometers), are being 
planned with this function in mind. 


2.2 TRENDS IN TOTAL OZONE 
2.2.1 Introduction 

This section summarizes the current understanding of trends in total column ozone. The objective of 
the analyses reported here are: (1) to re-examine the data set used by the Ozone Trends Panel (OTP) using 
several different statistical models which make varying assumptions in their treatment of the measurements, 
(2) to update trends derived from the Dobson network using the most recent data available, and (3) to 
examine the regional behavior of changes in total ozone using both ground-based and satellite-based 
measurements. 

It is essential to recognize the merits and limitations of both the Dobson network and orbiting remote 
sensors. The Dobson network provides a long-term data record whose duration is suitable for detailed 
statistical analysis. The deficiency is a limited geographic coverage, with a bias toward middle latitudes of 
the Northern Hemisphere. The data record obtained by satellites has global coverage, although at present 
the length of record is inadequate for derivation of long-term trends. In addition, the trends that exist in 
data sets obtained by satellites contain drifts of instrumental origin. A correction is therefore required, and 
at present this must be done in an approximate way by comparison with the Dobson network. However, 
new approaches to data processing may minimize this problem in the future. Nonetheless, satellite mea- 
surements are useful in identifying geographic variations in ozone change. 

A major effort in the Ozone Trends Panel Report centered on developing a provisionally revised 
Dobson data set. This made use of ozone measurements from nearby stations, temperatures measured at 
50 mb, and detailed studies of the records from each instrument to identify problems in the Dobson data 
set. In addition, comparison of co-located ozone measurements made by the Nimbus-7 TOMS instrument 
and the Dobson network allowed identification and correction of inconsistencies in the latter data base. 
Further information appears in Chapter 4 of the Ozone Trends Panel Report. The updated trend analyses 
reported here adopt the provisionally revised data base for each Dobson station through December 1986, 
the latest information available to the Ozone Trends Panel in early 1988. As of this writing, available data 
after 1986 have received a screening similar to that performed for the provisionally revised Dobson data 
set, and a number of adjustments have been made to data from 1987 and 1988 to smooth erratic variations. 
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2.2.2 Trend Analysis of Dobson Total Ozone Data 

2.2.2.1 Result of the 1988 Ozone Trends Panel and Comparison with Other Studies 

In 1988, the Ozone Trends Panel completed a statistical analysis of Dobson data from a network of 
stations, most of which were located between 30 and 64 degrees North latitude. The Panel focused on data 
over the period 1965-1986 for three reasons: It contained two complete solar cycles; effects of the 
atmospheric nuclear bomb tests should have subsided by 1965; and many more Dobson stations were 
operating in 1965 than in 1957-1958, thereby providing a more substantial data base than was available in 
earlier years. The percent changes in total ozone over the 17-year period 1970-1986 were reported by 
season and latitude zone. These results appear in the first column of Table 2.2-1, labeled as OTP. This 
work built on statistical studies of ozone performed over the last decade (e.g., Bloomfield et al., 1983; Hill 
et al., 1986; St John et al., 1981; Reinsel et al., 1987b). Over the Northern Hemispheric latitude band JO- 
64 degrees North, there was a wintertime change of about -4% for this period (ranging from - 2 % to -6% 


Table 2.2-1. Comparison of trend results by season and latitude (17 years % change between 1970 
and 1986 with one standard error in parentheses) 


Year Round 


Zone/ 

Latitude 

OTP 

AS1 

UW&C1 

UW&C2 

AS2 

53-64°N/55°N 

40-52°N/45°N 

30-39°N/35°N 

-2.3(0. 7) 
-3. 0(0. 8) 
-1. 7(0.7) 

- 1.6(0. 8) 
-2. 6(0. 7) 
-0.8(0. 8) 

-20(0.7) 

- 1. 9(0.6) 

- 1. 0(0.7) 

-2.0(0. 7) 
- 1.6(0. 6) 
- 1. 3(0.6) 

- 1.8(1. 0) 

- 1.9(0. 8) 

- 1. 9(0.9) 

Average 

-2. 3(0.5) 

- 1.7(0. 5) 

- 1. 7(0.5) 

- 1.6(0. 6) 

- 1.9(0. 8) 



Winter 




Zone/ 

Latitude 

OTP 

AS1 

UW&C1 

UW&C2 

AS2 

53-64°N/55°N 

40-52°N/45°N 

30-39°N/35°N 

-6. 2(1. 5) 
-4.7(1 .5) 
-2.3(1. 3) 

-4.2(1.6) 

-5.10.4) 

-1.70.3) 

— 4. 7(0. 8) 

— 3. 4(0. 7) 
-1. 5(0.8) 

-50(0.8) 
-3. 6(0.6) 
-20(0.7) 

-50(1.0) 
-4. 0(0. 7) 
-2. 8(0.8) 

Average 

— 4.4(1 .0) 

-3.70.0) 

— 3. 2(0. 6) 

-3.6(0. 6) 

-4. 0(0. 7) 



Summer 




Zone/ 

Latitude 

OTP 

AS1 

UW&C1 

UW&C2 

AS2 

53-64°N/55°N 

40-52°N/45°N 

30-39°N/35°N 

-0.2(0. 8) 

— 1 .9(0.7) 

- 1. 9(0.8) 

- 10(0.8) 
- 1. 2(0.7) 
-0. 6(0.8) 

— 0.7(0. 7) 

- 1. 5(0.7) 
-1. 2(0.8) 

- 0.9(0. 7) 

- 10(0.5) 

- 1.2(0. 6) 

-0. 6(1.0) 
- 1. 2(0.8) 
- 1. 2(0.9) 

Average 

-1. 3(0.5) 

- 1. 0(0.5) 

- 10(0.5) 

-10(0.5) 

-1.2(0. 8) 


Notes: 

OTP 

AS1 

UW&C1 

UW&C2 

AS2 


Ozone Trends Panel result based on zonal series. 

Allied-Signal results using zonal series. 

Universities of Wisconsin and Chicago results using regional/latitudinal analysis where the same trend is 

assumed to apply to all stations in a single region. 

Universities of Wisconsin and Chicago results using regional/latitudinal analysis with a linear trend in latitude. 
Allied-Signal results using regional/latitudinal analysis of individual stations including a nuclear weapons 
correction. The text describes additional differences between the AS2 and UW&C1 and UW&C2 models. 
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with increasing latitude), a summertime change of about - 1%, and a year-round change of about -2%. 
The statistical model included terms to account for changes in ozone related to the quasi-biennial oscillation 
(QBO) cycle in winds and the 1 1-year solar cycle. The long-term trends that remained, as reported in Table 
2.2-1, could not be explained by these known natural variations. 

Recent analyses by scientists at Allied-Signal, Inc. and the Universities of Wisconsin and Chicago 
have subsequently verified the findings of the Ozone Trends Panel, including a new analysis of the sensitivity 
of derived trends to (I) the length of data records, (2) regional differences, and (3) the statistical techniques 
applied. Of particular interest were the major perturbations in ozone before 1965 and after 1982 and their 
effects on derived trends, as well as the differences observed among different geographical regions (e.g.. 
North America, Europe, and Japan). The nuclear weapons tests in the early 1960s and natural phenomena, 
such as the El Chichon volcanic eruption and El Nino events in the 1980s, may be related to these ozone 
fluctuations. Figure 2.2-1 is a time history of residual zonal mean total ozone in Dobson Units (DU) with 
most known natural and man-made variations subtracted out. This shows some large perturbations in the 
post- 1982 period following, but not necessarily related to, the El Chichon volcanic eruption. With respect 
to regional differences. Figure 2.2-2 (a)-(c) shows larger decreases over North America and Europe than 
over Japan after allowing for known variability. The nonrandom fluctuations that remain in Figure 2.2-2 
arise from sources of variation not included in the statistical model, for example the El Nino oscillation, 
and an incomplete subtraction of known sources such as the QBO. However, these nonrandom fluctuations 

are removed by the autoregressive filter. The European wintertime dip in 1982-1983 is particularly notice- 
able. 


Table 2.2-1 presents a comparison of five separate sets of trend results, including the Ozone Trends 
Panel findings, over the 17-year period 1970 through 1986. In these and all subsequent Dobson analyses, 
the statistical models account for variations in ozone related to the solar cycle, the QBO, and, where 
applicable, nuclear weapons tests during the 1960s. The derived long-term trends cannot be explained by 
any of these mechanisms. The time dependence of ozone changes related to nuclear weapons tests is based 
on calculations by the Lawrence Livermore Laboratory chemical-transport model (Wuebbles and Kinnison, 
1989). However, the magnitude of the ozone change is determined from the statistical analysis and is 
allowed to vary with latitude. 

The choice ot the period 1970-1986 for the trend assessment corresponds to the period during which 
theory predicts the onset of effects on the ozone layer related to chlorofluorocarbon (CFC) release. The 
trend is assumed to start in January 1970 and be zero prior to this time. Sensitivity to the starting date is 
discussed later. The form of the trend model used in all of the studies is described in Appendix 2. A. 

The results in the second column of Table 2.2-1, labeled AS1, were developed by Allied-Signal 
researchers with the same zonal mean series of total ozone data used by the Ozone Trends Panel for latitude 
belts 30-39 N, 40-52 N and 53-64°N. Results in the third and fourth columns were obtained by researchers 
from the Universities of Wisconsin and Chicago (designated UW&C). Unlike the first two sets of results, 
which analyzed the data after it was pooled into the zonal mean series shown in Figure 2.2-3, the UW&C 
studies derived a separate trend of each of the 25 stations within the 26-60°N latitude band (see Table 2.2- 
2, excluding Reykjavik and Uccle). The trends were then fitted (or regressed) against latitude using two 
different methods. In the first method, labeled UW&C1 , all stations within each of three different latitude 
zones, 30-39°N , 39-50°N and 50-64°N, were assumed to have a common average trend. In the second 
method, labeled UW&C2, the trends were regressed as a linear function of latitude (Bojkov et al., 1989). 
Both UW&C methods allow for apparent regional differences over North America, Europe, Japan, and 
the Middle East. The first method more closely corresponds to the OTP and AS I analyses based on the 
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Figure 2.2-2a. Regional plots of Dobson total ozone after removal of seasonal, solar, QBO, and nuclear 

t fl 6 !! eC S f ° r N ° rth America ( t0 P)> Eur ope (center), and Japan (bottom), all months. The deviations are 
p otted in Dobson units where the upper panels give monthly values and the lower panels are smoothed. 
Nonrandom fluctuations are removed from the trend analysis by an autoregressive filter. Vertical lines at 
the end of 1986 denote the end of the data set available to the Ozone Trends Panel. 


zonal mean series, while the second method provides a capability for estimating trends as a function of 
latitude. Overall the second method uses more information contained in the data than does analysis of a 
zonal series, but it also makes a stronger assumption concerning the relationship of trends to latitude, 
namely that they are linear. In Table 2.2-1 the trend estimates are given for three latitude zones based on 
UW&C method I , and the predicted changes are given for latitudes 55°N, 45°N, and 35°N based on UW&C 
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Figure 2.2-2b. Same as 2.2-2a for winter months (December-March). 


method 2. Separate trends are derived for year-round, winter (December-March), and summer (May- 
August), respectively. 

The last column in Table 2.2-1, labeled AS2, presents results obtained by Allied-Signal researchers 
using a technique similar to method 2 in the UW&C study but which includes adjustments for nuclear 
weapons tests in the 1960s. Results are based on only the 25 stations located at latitudes between 26 N and 
64°N (Table 2.2-2). The AS2 results do not use the two stations from Egypt and Pakistan to form a separate 
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SUMMER 
NORTH AMERICA 



SUMMER 



SUMMER 



Figure 2.2-2c. Same as 2.2-2a for summer months (May-August). 


Mid-East region because of the disproportionate influence these two stations can have on the summary 
results. 1 


Overall, the five sets of results show a consistent pattern. They indicate a wintertime (December- 
March) change ranging from -2% to -6% between 1970 and 1986 which correlates with latitude, that is, 
the wintertime trends become more negative with increasing latitude. The value derived by the Ozone 
Trends Panel for winter near 55°N latitude tends to be 1-2% more negative than found in the other studies. 
Also, standard errors associated with the wintertime trends are smaller in those studies that accounted for 
the additional regional and latitudinal information. 
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ZONAL AVERAGE SERIES 



30° - 39° N 



Figure 2.2-3. Plots of Dobson total ozone monthly averages combined into time series for each latitudinal 
zone 53-64°N, 40-52°N, and 30-39°N. The vertical scale is in Dobson units. 


Although for a given latitude, the five sets of summertime trend estimates may differ by as much as 
about 1% the change in summertime averaged over all latitudes was about \/c in each study, 
significant latitudinaf pattern appears in the summertime trends. Similarly, in the year-round trends, no 
significant latitudinal pattern appears. The concept of “significance” here refers to the stat.s ‘ ica > ^ se 
exceeding two standard errors. For each study, the average year-round change is approximately 2% o 

17 years. 

The difference between the two studies which used zonal series (labeled OTP and AS I in Table 2.2- 
1) arise mostly from the differing lengths of record. For example, if the pre-1965 data and nuc ear test 
adjustment are excluded from the AS1 analysis, as was done in the Ozone Trends Panel study, then the 
revised AS1 results would be close to the OTP trends. Other differences among the groups ot results m 
Table 2 2-1 can be ascribed to the inclusion of a nuclear weapons adjustment, which was used in the 
columns marked AS1 and AS2, the use of four versus three regions and slightly different groups of stations 
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Table 2.2-2. List of Dobson total ozone stations at latitudes between 26°N and 64°N (data from these 
stations were revised by the Ozone Trends Panel through 1986) 


Region 

Station 

Latitude 

Start Date 

End Date 

Now available 

N. America 

Churchill 

59°N 

1/65 

12/86 

9/88 


Edmonton 

54°N 

1/58 

11/86 

3/88 


Goose 

53°N 

1/62 

11/86 

9/88 


Caribou 

47°N 

5/58 

12/86 

10/88 


Bismarck 

47°N 

1/58 

12/86 

10/87 


Toronto 

44°N 

1/60 

11/86 

9/88 


Boulder 

40°N 

1/64 

12/86 

10/88 


Wallops Island 

38°N 

1/70 

12/86 

10/88 


Nashville 

36°N 

1/63 

12/86 

10/88 


Tallahassee 

30°N 

5/64 

11/86 

4/88 

Europe 

Reykjavik 

64°N 

1/58 

10/86 

10/88 


Lerwick 

60°N 

1/58 

11/86 

10/88 


Leningrad 

60°N 

8/68 

12/85 

12/85 


Belsk 

52°N 

1/63 

12/86 

8/88 


Bracknell 

51°N 

1/69 

12/86 

10/88 


Hradec Kralove 

50°N 

1/62 

6/86 

10/88 


Uccle 

50°N 

7/71 

12/86 

8/88 


Hohenpeissenberg 

48°N 

1/67 

12/86 

1 1/88 


Arosa 

47°N 

1/58 

12/86 

10/88 


Vigna DiValle 

42°N 

1/58 

12/86 

1 1/87 


Cagliari 

39°N 

1/58 

12/86 

6/88 

Japan 

Sapporo 

43°N 

1/58 

12/86 

10/88 


Tateno 

36°N 

1/58 

12/86 

10/88 


Kagoshima 

32°N 

1/58 

12/86 

10/88 


Naha 

26°N 

4/74 

12/86 

10/88 

Mid East 

Cairo 

30°N 

11/74 

10/86 

10/88 


Quetta 

30°N 

1/58 

12/86 

10/88 


' . &C '^"p VerSUS ^ S2) ’ the analys,s of averaged zonal data versus the analysis of individual 

nR ‘° ns (OTP and versus UW&C2 an ^ AS2), and small effects due to the choice of solar cycle and 

QBO adjustments. The essent.al information in Table 2.2-1 is that all of the analysis techniques produce 
similar latitudinal and seasonal patterns. 

Trends derived from the AS2 model are shown on a station-by-station basis in Figures 2 2-4(a)-(c) 
Throughout the remainder of this section, trends will be expressed in percent per decade. This allows for 
easy comparison of trends derived from data records of differing lengths. For the year-round case all the 
Morth American (N) trends are negative, all the European (E) trends are negative except Reykjavik at 64°N 
and the four Japanese (J) trends are balanced on either side of zero change.The most northerly Japanese 
station, Sapporo, lies at 43°N. For the wintertime period, a strong latitudinal pattern emerges with the most 
negative trends appearing at the most northern latitudes. Only the Japanese stations show little or no 

wintertime change. In the summertime, the trends, although tending toward the negative side, show no 
latitudinal gradient. 
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Figure 2.2-4. Individual trend estimates for the period 1 970-1 986 using the AS2 model of Allied-Signal, 
Inc. Trends appear versus latitude and are expressed in % per decade for (a) year-round, (b) winter 

(December-March), and (c) summer (May-August) periods. The plotting symbol is E = Europe, N = North 

America, and J = Japan. 


Fisletov (1989) performed an independent trend analysis of the provisionally revised data from 30 
Dobson stations and, in addition, examined the effect of including measurements from the 22 most reliable 
M-83 instruments in the USSR network. Trends based only on the Dobson record through 1986 lie in the 
range defined in Table 2.2-1 and show the same latitudinal and seasonal pattern. 
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2.2.2.2 Sensitivity Analysis 

To quantify the dependence of the derived trends on factors such as data length, nuclear weapons 
adjustment, type of solar adjustment, inclusion of data from Soviet instruments, and nonseasonal versus 
seasonal trend models, a sensitivity analysis was done. The effects on the trends are summarized in Table 
2.2-3 using the Allied-Signal models as baselines. The sensitivities are expressed as the change in trend, in 
% per decade, due to the factors considered. The results are similar to those obtained by the Ozone Trends 
Panel. 


The exclusion of post-1982 data has a large effect as shown by item 1 in Table 2.2-3 and Figure 2.2-5 
(a)-(c). The trends are noticeably less negative or more positive when only data through October 1982 are 
used. This date corresponds approximately to the arrival of the El Chichon aerosol at middle latitude. 
Interestingly, the effect on the trends from eliminating post-October 1982 data is larger in summer than in 
the winter. This is particularly evident when the data are analyzed using a regional/latitudinal model where 
the effect of dropping post-October 1982 data is to produce a summertime trend 0.8% per decade less 
negative and a wintertime trend 0.3% per decade less negative. The effects in summer and winter are both 
about 1% per decade using the zonal analysis. After 1982, summertime ozone in all three regions (North 
America, Europe, and Japan) went down and stayed down until 1987, whereas wintertime ozone showed 
some degree of recovery in all three regions (see Figure 2.2-2). Comparison of Figure 2.2-5(a)-(c) with 
Figure 2.2-4 (a)-(c) shows how the trend patterns change when the post-October 1982 data are dropped. 
This behavior may be a chemical response to volcanic aerosol from the El Chichon eruption or to a true 
change in the nature of the long-term trend in ozone. Since derived trends are highly influenced by extreme 
values or sustained dips at or near the end of the record, the addition of extra data from 1987 and 1988 is 
important to see if this nonlinear pattern is continuing (see Sections 2. 2. 2. 3 and 2. 2. 3. 2.). 


Table 2.2-3. Summary o f sensitivity analysis (% per decade effect on trends) 



Zonal Analysis' 


Regional/Latitudinal Analysis 6 

Sensitivity Scenario 

Year-Round 

Winter 

Summer 

Year-Round 

Winter 

Summer 

1. Exclude Post 10/82 Data 

+ 0.8 

+ 0.9 

+ 1.0 

+ 0.3 

+ 0.3 

+ 0 8 

2. Uniform Model 

+ 0.6 



— 

+ 0.4 



3. Exclude Pre-1965 Data 

-0.6 

-0.6 

-0.5 

-0.3 

-0.1 

-0 4 

4. Exclude Nuclear Term 

+ 0.2 

+ 0.2 

+ 0.1 

+ 0.1 

0 0 

0 0 

5. Add Russian Stations 

-0.1 

+ 0.4 

-0.6 

NA 

NA 

NA 

6. Smooth Solar Term 

0.0 

0.0 

0.0 

0.0 

0.0 

0 0 

7. Remove Solar Term 

-0.1 

-0.1 

-0.1 

-0.1 

-0.2 

-0 1 

8. Remove QBO 

+ 0.1 

0.0 

+ 0.1 

0.0 

-0.1 

+ 0 I 

9. Start Trend in 1/65 vs. 1/70 

+ 0.4 

+ 0.7 

+ 0.3 

+ 0.2 

+ 0.5 

+ 0 2 

10. Start Trend in 1/75 vs. 1/70 

-0.5 

-0.9 

-0.4 

-0.4 

-0.6 

-0.4 

11. Add 1987 and 1988 Data 

NA 

NA 

NA 

+ 0.2 

+ 0.1 

+ 0.3 

Sensitivity effect averaged across three zones 30-39°N, 40 

-52°N, 53-64°N usina AS1 model. 



^Sensitivity effect predicted at 45°N using AS2 model. 






Note: Rounding is to nearest 0. 1 % decade. 






The table entry indicates the amount (in % decade) to add to the trend from the standard 

scenario to obtain the trend for the 

sensitivity scenario. 







Example: if data after 10/82 are omitted from the AS1 standard calculation, 

the zonal year-round trend increases in a nnsitivp 

direction by +0 8% decade; i.e., from - 

1.0 (standard) to - 

0.2 (sensitivity calculation). 




Winter = December-March; Summer = May-Auqust 
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Figure 2.2-5. Individual trend estimates for the period 1 970-1 986 using the AS2 model of Allied-Signal, 

Inc Same as Figure 2.2-4 except trends are fitted only through October 1 982. Trends appear versus latitude 
and are expressed in % per decade for (a) year-round, (b) winter (December-March), and (c) summer (May- 
August) periods. The plotting symbol is E = Europe, N = North America, and J - Japan. 


If the statistical model assumes that trends are the same for all months (i.e., a uniform trend), then 
the trends are closer to zero by about 0.5% per decade. Analyses suggest that the uniform trend model is 
not appropriate for use on data poleward of 40°N, since there are statistically significant winter versus 
summer differences in the trends. 
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As mentioned earlier, exclusion of pre-1965 data because of possible nuclear effects and the paucity 
of data leads to more negative trends. In a regional latitudinal analysis, the effect on the winter trends is 
small. However, in a zonal analysis like that done by the Ozone Trends Panel, the wintertime change is 
made more negative by about 0.6% per decade when the pre-1965 data are excluded. So, depending on the 
formulation of the model and the way data are handled, the wintertime trend can be cjuite sensitive to the 
re- 1 965 data. 


An analysis by Allied-Signal researchers shows that inclusion of data from Soviet instruments in the 
zones poleward of 40°N makes the winter trend less negative or more positive by about 0.4% per decade 
when averaged across the latitude belt 30-64°N, but enhances the summer reduction, with little effect on 
the year-round trend. 


Exclusion of the nuclear adjustment makes the trends less negative by about 0.2% per decade. The 
treatment of the solar cycle and the QBO in winds has little or no effect on the trends. 

The sensitivity of the derived trend to the start date has interest because of issues related to (I) the 
trends being larger after October 1982 than in the earlier period, and (2) the trend start point being at or 
near extremes of the solar cycle. As pointed out above, the trend estimates are relatively insensitive to the 
solar cycle. However, as seen in Table 2.2-3, the trend estimates are sensitive to the start date. For example, 
in the zonal analysis, the derived winter trend is less negative by 0.7% per decade when the trend is started 
5 years earlier (i.e. , January 1965) and more negative by 0.9% per decade when the trend is started 5 years 
later (i.e., January 1975). This could be due in large part to influence of the post-October 1982 data. That 
is, the influence of the post-October 1982 data becomes larger when the trend period is taken to be shorter 
relative to the baseline period. 

In the sensitivity analysis shown in Table 2.2-3, it should be noted that, in most cases, the regional/ 
latitudinal analysis is less sensitive to the factors considered than is the zonal series analysis. Hence, the 
former approach is more robust and is currently the analysis of choice for the Allied-Signal and Wisconsin- 
Chicago research groups. Another reason for its preference is that variability within and across regions 
can be taken into account along with the latitudinal relationship. 

In a separate regional sensitivity analysis (Table 2.2-4) using the Wisconsin-Chicago model, a com- 
parison of trends was made using both the ‘ provisionally revised” Dobson data set and the Dobson data 
as originally published. This showed that over Europe the “provisionally revised” data gave trends that 
were about 1% per decade more negative in all seasons than did the original data. Much smaller differences 
between trends in the two data sets resulted over North America and Japan. The differences in trends 
caused by the data revisions and especially by the seasonal versus uniform trend model explain why earlier 
published studies (Hill et al., 1986, Reinsel and Tiao, 1987a) show no significant trends, since the previous 


Table 2.2-4. Regional sensitivity analysis comparing trends from the “Provisionally Revised’’ Dob- 
son data set with the original published data 



Trend Difference 

= (Revised - Original) in % per Decade 


Japan 

North America 

Europe 


(26°N-43°N) 

(30°-59°N) 

(39°N-60°N) 

Year-Round 

0.2 

0.0 

- K0 

Winter 

0.2 

0.4 

- 1.2 

Summer 

0.2 

-0.2 

-0.8 
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studies used a uniform trend model on the original Dobson data. Additional discussion of this issue appears 
in Chapter 4 of the Ozone Trends Panel Report. 

In summary, the choice of model (seasonal versus uniform, and zonal versus regional/latitudinal), 
timing of the trend starting date, and the record length have the largest effects on the derived trends. 
However, the general pattern of significant wintertime changes at the middle and higher latitudes remains 
intact under different scenarios. 

2.2.2.3 Update of Trends into 1988 

As of this writing, data were available into 1987-1988 from all but one of the Dobson stations in Table 
2.2-2 (Leningrad). In most cases, data were available through August 1988 or later (see last column of 
Table 2.2-2). When these data were included in the trend analysis, most of the year-round, winter, and 
summer trends became slightly less negative. The effects of the additional data are shown in Table 2.2-3 
based on the AS2 model, and the individual station trends appear in Figure 2.2-6 (a)-(c) for year-round, 
winter, and summer, respectively. At 45°N, the winter trends are less negative by 0.1% per decade. At 
55°N (not listed in Table 2.2-3), the winter trend is less negative by 0.2%. The biggest change was seen in 
the summer, where the trend at 35°N is calculated to be 0.5% per decade less negative. 

Trends in total ozone tend to be less negative in all seasons when the 1987 and 1988 data are added. 
Table 2.2-5 summarizes the trend estimates updated into 1988 by latitude alongside the estimates through 
December 1986, using the AS2 model. Figures 2.2-7(a) through (c) show the differences between trends 
estimated into 1988 versus those through December 1986 at each station. Notice that at most sites the trend 
differences are positive with the addition of the new data. Overall, the year-round trend estimates predicted 
at 35°N and 45°N are no longer statistically significant. The winter trend estimates updated into 1988 at the 
latitudes 45°N and 55°N are significant but closer to zero than those through December 1986, while the 
trend at 35°N is slightly more negative. Table 2.2-5 shows that the addition of data into 1988 does not alter 
the general conclusions of the Ozone Trends Panel, although differences in detail exist. 

Table 2.2-6 summarizes the regional effects of including the 1987 and 1988 data by tabulating the AS2 
model trends using data through 1986, together with the results using all data available into 1988. In both 
situations, the Japanese region shows no significant trend year-round, winter or summer. In both cases the 
year-round and winter trends are significant over North America and Europe, and the summer change is 
significant over North America. In most cases the changes are slightly smaller when the data set is updated 
into 1988 as compared to changes derived through 1986. 

Although including or excluding the solar cycle term in the statistical model has little effect on the 
estimated linear trend in ozone (see Table 2.2-3), there appears to be a measurable solar effect on ozone 
itself. This effect has no apparent relationship to latitude as demonstrated in Figure 2.2-8. The average 
solar coefficient over the 25 Northern Hemisphere stations is 0.03 Dobson units per unit of 10.7 cm flux 
(standard error 0.005). This corresponds to a change in ozone of approximately 1.2% from a typical solar 
maximum to solar minimum, a range taken of 140 units as measured by the 10.7 cm flux. 


2. 2. 2.4 Summary of Dobson Analyses and Comparison with Theory 

Recent studies have both verified the findings of the Ozone Trends Panel regarding changes in Northern 
Hemisphere total ozone over the period 1970 through 1986 and extended the analyses using data through 
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Figure 2.2-6. Individual trend estimates for the period 1970-1988 using the AS2 model of Allied-Signal, 
Inc. Same as Figure 2.2-4 except trends are fitted using updated data into 1988. Trends appear versus 
latitude and are expressed in % per decade for (a) year-round, (b) winter (December-March), and (c) summer 
(May-August) periods. The plotting symbol is E = Europe, N = North America, and J = Japan. 


late 1988. Over the latitude band 26-64°N, the new analyses (see Table 2.2-6) found a wintertime trend 
averaging about —2.2% per decade for the period 1970—1986, with summertime and year-round trends 
averaging about — 0.7% and — 1.1% per decade, respectively. These results are slightly less negative than 
those of the Ozone Trends Panel. The winter reductions become larger with increasing latitude, varying 
from -1.7% per decade near 35°N to -3.0% per decade near 55°N. 
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Table 2.2-5. Trends in total ozone derived from Dobson data updated into 1988 compared with 

trends derived through 1986 (trends are in % per decade with one standard error in 

parentheses). Winter - December— March ; Summer = May— August 


Year Round 


Latitude 

55°N 

45°N 

35°N 

1986 

-1.08L56) 

- 1. 1 0( .49) 

— 1.1 1 ( .53) 

Winter 

1988 

- .99(.58) 

- ,94(.50) 
-,88(.54) 

Latitude 

55°N 

45°N 

35°N 

1986 

-3.02t.59) 
-2.34t.41) 
- I.66(.48) 

Summer 

1988 

-2.70L58) 

-2.26(.47) 

-1.8K.53) 

Latitude 

55°N 

45°N 

35°N 

1986 

— .35(.60) 

- .68(.45) 
- I.0K.5I) 


1988 

- .21 ( .59) 

- ,34( .48) 

- .47(.54) 


Diff. 

.1 

.2 

'y 


Diff. 

.3 

.1 


Diff. 

.1 

.3 

.5 


A variety of sensitivity studies were performed to determine which features of the ozone data and 
the analysis techniques most affected the derived trends. The choice of model (seasonal versus uniform 
and zonal versus regional/latitudinal), timing of the trend starting date, and the record length have the 
largest influences on the derived trends. However, regardless of the model or the particular treatment of 
the data, significant negative trends appear in winter at middle to high latitudes. 

At the time of this writing, total ozone data are available through late 1988 at many stations. Year- 
round and winter trends computed using the most recent data (middle to late 1988) are slightly less negative 
than the trends derived through 1986. while summer trends were noticeably less negative at latitude 35°N. 
Winter trends based on the updated data continue to show an increase with latitude, from about - 1.0% 
per decade near 35°N to -2.7% per decade near 55°N. 

The executive summary of the OTP report presented predicted change in total column ozone based 
on two-dimensional models for summer (May-August) and winter (December-March). Results appeared as 
the difference between computed column ozone averaged over the 11-year period 1976-86 and that for 
1965-75. For summer the computed difference imply a trend of - 0.4 to - 0.6% per decade over the latitude 
range 35°N to 55°N. These values lie within one standard error of the derived summer trends in Table 2.2- 
5. Predicted trends during winter are more negative than those in summer, being approximately -0.5, 
-0.8, and - 1.1% per decade at 35°N, 45°N, and 55°N, respectively. The corresponding observed trends 
in Table 2.2-5 are - 1.8 ± 1.0, -2.3 ±1.0, and -2.7 ±1.2% per decade, where the error bars are 2 a, 
95% confidence limits. Theory and observation both yield trends which become more negative as one 
moves poleward. However, in all cases the models produce trends that are less negative than the obser- 
vations and lie outside of the 2 o error bars. 
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Figure 2.2-7. Differences between estimated trends into 1988 versus trends through 1986. The difference is 
calculated as the trend including data into 1988 minus the trend through 1986, expressed in % per decade at 
each Dobson station between 26°N and 64°N (a) year-round, (b) winter (December-March), and (c) summer 
(May- August). 


2.2.3 Analysis of Satellite-Based Data from SBUV and TOMS Beginning in November 1978 
2.2.3.1 Data Base and Analysis Procedure 

The SBUV and TOMS sensors on the Nimbus-7 satellite began operation in November 1978. Given 
the short length of these data bases and other problems discussed below, it is not feasible to derive true 
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Table 2.2-6. Trends in total ozone derived by region based on data through 1986 and data updated 
through October 1988 (trends are in % per decade with one standard error in paren- 


Region 

Number 
of Stations 

Year Round 
1986 1988 

Winter 

1986 

1988 

Summer 
1986 1988 

N. America 

10 

- 1.7(.3) 

-1.6(.3) 

-2.9(A) 

2.91.3) 

-1.41.4) 

-1.21.4) 

Europe 

11 

- 1 -5(.3) 

- 1.4(.3) 

-3.01.4) 

2.91,4) 

-0.61.4) 

-0.31.3) 

Japan 

4 

+0.21.5) 

+ 0.51.5) 

-0.51.6) 

0.61.5) 

+ 0.31.7) 

+ 0.81.6) 

26-64°N 

25 

-1.K.5) 

— 0.9( .5) 

-2.21.6) - 

2.21.6) 

-0.71.4) 

-0.31.5) 


0.08 
0.07 
0.06 H 
0.05 
0.04 
0.03 
0.02 
0.01 
0.00 
- 0.01 
- 0.02 
* 0.03 
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Figure 2.2-8. Solar coefficient (Dobson units per unit of 1 0.7 cm flux) computed for 25 Northern Hemisphere 
stations from the AS2 model of Allied-Signal, Inc. The plotting symbol is E — Europe, N = North America, 
and J = Japan. The average solar coefficient is 0.030, with a standard error of 0.005. 


long-term trends from the satellite record. However, satellite data cover the entire sunlit portion of the 
globe and, subject to the constraint of a short duration, provide valuable information on the geographic 
patterns of ozone change. Here the satellite data complement the Dobson record, which indicates regional 
variability in ozone trends but has insufficient geographic coverage to allow detailed analysis of this 
structure. Figure 2.2-9 (a)-(b) illustrates the geographic variation in ozone changes based on TOMS results 
normalized to the Dobson network to remove instrument drifts. The plots give the percent differences 
between annual mean ozone over the period 1987-88 and over 1979-80. Decreases appear in almost all 
regions. On the Northern Hemisphere (Figure 2.2-9a) the largest ozone decreases, -4.5 to -7.5%, appear 
at mid-latitudes, especially over the Pacific Ocean. The most prominent feature in the Southern Hemisphere 
is the ozone decrease which becomes more pronounced as one moves toward the polar region. 

The remainder of this section examines the spectral variability of changes in ozone over an 8-year 
period beginning in 1978 based on SBUV, TOMS, and the Dobson network. The results should be used 
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YEARLY AVERAGE TOMS 1987-1988 MINUS 1979- 1980 


PERCENT 




Figure 2.2-9a. Geographic structure in total ozone changes based on TOMS. Results are percent 
difference in annual mean ozone for 1987-1988 and 1979-1980 for Northern Hemisphere. 


only as a measure of geographic structure in ozone changes, not as trends in the same sense derived from 
the entire long-term Dobson record. 

Instrument degradation is a major issue in the analysis of satellite data sets over multi-year periods. 
Flcig et al. (1 986), Me Peters and Komhy r ( 1 989b), and Bojkov et al. ( 1 988) have performed a similar analysis 
for TOMS. Although drifts in the satellite data sets have been removed by comparison to the Dobson 
network, the correction involves a degree of uncertainty. For this analysis, the SBU V data were normalized 
to the Dobson record by fitting a single linear trend to all coincident measurements over the 8-year observing 
period. The TOMS data were normalized to Dobson by comparing coincident measurements during single 
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Figure 2.2-9b. Same as Figure 2.2-9a for Southern Hemisphere. 

months and deriving a linear trend for each month of the year (R.S. Stolarski, private communication, 
1989). 

The information of interest is the geographic pattern of changes in total ozone over 8 or 10 years, 
corrected for solar cycle dependence. A problem arises here in that the data records are less than one solar 
cycle in length. In particular, the 10.7 cm flux (F,«. 7 ) decreases markedly during the period 1978 through 
1986. As a result, initial test calculations showed unacceptably high correlations between total ozone and 
F 107 . To circumvent this situation, the long-term Dobson record was used to determine the historical 
relationship between F,« 7 and total ozone (see Figure 2.2-8). This correlation provided a basis tor removing 
the solar cycle effect from the SBUV and TOMS data sets. The changes in ozone over 8 years, corrected 
in this way, were compared to those in which the F l07 correlation was removed from the satellite time 
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series directly. The geographic pattern of ozone change was the same in both cases. This arises from the 
absence of a clear latitudinal dependence in the solar cycle effect, as shown in Figure 2.2-8. The absolute 
values of the changes in total ozone, however, depend on the treatment of the solar cycle correlation. This 
result again demonstrates that the available satellite data are most useful for examining geographic patterns 
of ozone change. As the data base grows to span more than one solar cycle, studies of true long-term 
trends will become feasible. 

For analysis of the SBUV and TOMS data, the statistical model described in Appendix 2. A was 
revised so that total ozone, adjusted for the solar cycle, was the dependent variable. In addition, the 
analysis omitted autoregressive structure in the noise and the term related to the QBO. To account for 
different variances in a given month of different years, the model was fitted by a weighted estimation 
procedure with weights for each month inversely proportional to the monthly variances. Finally, the 
monthly trends were combined into seasonal values by taking the simple averages of December through 
March as winter and May through August as summer. 


2.2.3.2 Geographic Patterns of Ozone Change and Comparison with the Dobson Network 

As part of the NOAA/National Weather Service (NWS) continuing Stratospheric Monitoring Program, 
the Climate Analysis Center has produced monthly global analyses of the SBUV total ozone data (Nagatani 
et al., 1988). These use the NWS standard rectangular grid array on a polar stereographic projection where 
the grid size is approximately 3 degrees latitude by 3 degrees longitude at middle latitudes. For each 
hemisphere the grid is composed of 4,225 points total. 

Analyses of total ozone from SBUV were done utilizing the statistical model outlined above for each 
grid point in each hemisphere. The monthly ozone changes, corrected for solar cycle dependence, were 
then combined into four winter/summer categories. Results for the Northern Flemisphere appear in Figure 
2-2-10 (a)-(b). All contours have the units of percent change per decade, although the data base encompasses 
only an 8-year period. Note that the analysis ceases at 60°N. This is because the SBUV sensor can acquire 
data only in daylight. Since the solar terminator moves in latitude with the seasons, 60°N is generally the 
highest latitude at which data are available during all months. 

Figure 2.2-10 (a) and (b) show that changes in ozone deduced in winter tend to be more pronounced 
and have more regional variability than those in summer, although substantial variations in space exist in 
both seasons. During winter in the Northern Hemisphere, the changes vary from about -4% per decade 
over northwest U.S.S.R. to + 4 % per decade over northeastern North America, with positive changes 
over the northwest Pacific Ocean and eastern Asia. Note also, the substantial negative changes in the mid- 
latitude ocean areas where no ground-based observations exist. Summertime values also show a variance 
around the hemisphere, but with a reduced range. A problem that arises with a short data set is that the 
derived patterns can be influenced significantly by movement of wave structure in ozone. Derived changes 
may represent a geographic redistribution of ozone, rather than a net increase or decrease. 

Results for the Southern Hemisphere appear in Figure 2.2-1 1 (a)-(b). While spatial variability around 
latitude circles is reduced compared to the Northern Hemisphere, many of the statements for the north 
apply to the Southern Hemisphere as well. One item of particular note is that the largest changes in ozone 
are over the ocean areas. This demonstrates the value of the extensive coverage obtained by a satellite- 
borne instrument. 
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The geographic structure in Figures 2.2-10 and 2.2-11 implies that zonal mean changes in ozone can 
repiesent small differences between larger pos.tive and negative quantities. The zonal mean results based 
on X years ot SBUV measurements indicate decreases in summertime ozone over the Northern Hemisphere 
w ic exceed those in winter. This result, which appears contrary to the trends in Table 2 2-5 arises from 
the short duration ot the data set. Reference to Figure 2.2-2 (c) reveals a dip in ozone beginning in 1982 
which fasts into 986 over North America and Japan. This period coincides with the final 4 years of the 
^BU V data set. A recovery in summertime ozone occurs in all regions after the end of the SBUV record. 

An analysis similar to that presented above is possible using data from the TOMS instrument. 
However, for the purpose ot this section it is sufficient to show that both SBUV and TOMS detect consistent 
changes m ozone at selected points on the globe. The TOMS data used in this analysis were provided by 
R.S. Stolarski and R.D. Hudson (1989, private communication) and have been adjusted by reference to 
the Dobson network as mentioned previously. These data were placed on a two degree latitude by five 
degree longitude grid. Comparisons between ozone changes deduced from SBUV and TOMS were done 
at I, 4 geographic locations eight in the Northern Hemisphere and six in the Southern Hemisphere, as listed 
in I able -w. For each ot the 14 grid points the regression model was applied as for SBUV. 


Figure — -I- Presents the results derived from both SBUV and TOMS at each of the 14 locations 
over the 8-year period November 1978 through October 1986. The two sets of results track each other very 

r" ,r , h , th V largeSt dlfferencc being aboul 3% P er decade at point number 9 (latitude 60°S, longitude 
150 E). This discrepancy may arise from the better horizontal resolution and coverage of TOMS in cases 
where tight gradients exist. I he similarity between curves in Figure 2.2-12 supports the reality of the global 
patterns in column ozone change derived from SBUV. The TOMS data set extends later in time than that 

from SBUV. An analysis of the I()-yearTOMS record, through October 1988, produced a pattern which is 
generally consistent with that in Figure 2.2-12. 


Although both , he SBUV and TOMS da, a am have been adjusted ,o conform in a general sense with 
Ihc Dobson network, this procedure does not force agreement by individual region. It is therefore of interest 


Table 2.2-7. 


Latitudes and longitudes of the 14 grid points selected for comparison of 
ozone derived from SBUV and TOMS 


changes in 


Point 

Number 

I 

7 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 


Latitude 

54°N 

34°N 

34°N 

34°N 

I4°N 

40°N 

44°N 

40°N 

60°S 

44°S 

I0°S 

44°S 

24°S 

44°S 


Longitude 

50°W 

55°W 

I30°W 

80°E 

25°W 

I20°E 

I00°W 

20°E 

150°E 

I00°W 

40°W 

I20°W 

50°W 

20°W 


Season 

winter 

winter 

winter 

winter 

summer 

summer 

summer 

summer 

summer 

summer 

summer 

winter 

winter 

winter 
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Figure 2.2-12. Changes in total ozone based on eight-year data sets from the SBUV and TOMS 
instruments. Values refer to the 14 sites listed in Table 2.2-7 and are expressed in percent per decade, 
although the data sets cover only the eight-year period November 1978 through October 1986. 


to examine the compatibility of the ground-based and satellite results over the short time period encom- 
passed by the latter data set. The investigation uses Dobson ozone measurements only over an 8-year 
period (November 1978-October 1986). As such, the concept of a long-term trend should not be applied. 
The focus here is on the consistency of geographic patterns in ozone change derived from independent 
measurement systems. This analysis applies a statistical regression model described in Appendix 2. A to 
the 25 Dobson stations listed in Table 2.2-2 which are located in Japan (J), Europe (E), and North America 
(N). 


Figure 2.2-13 (a)-(b) presents results for winter and summer over the 8-year period, where the change 
in percent per decade for each station appears at the appropriate latitude. Comparison of the winter map 
of Figure 2.2-10 (a) with its Dobson counterpart shows several interesting features. The Japanese stations 
indicate a change which becomes less negative with increasing latitude. This corresponds well with the 
latitudinal variation derived from SBUV over this region, with an increase in ozone appearing northwest 
of Japan. European Dobson data generally indicate the reverse, with a strong negative change with 
increasing latitude. Changes derived from the Dobson instruments are somewhat larger in magnitude, but 
reveal a latitudinal pattern in agreement with SBUV. Finally, over North America ozone changes derived 
from Dobson instruments show no obvious pattern. Based on the map from SBUV, this may arise from 
the fact that this region is divided into a positive northeast sector and a negative western region. 

For the summer season. Figure 2.2-13 (b) shows ozone changes over Japan which are virtually 
independent of latitude. The result from SBUV in Figure 2.2-9 (b) shows the -4% per decade contour 
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LATITUDE 
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LATITUDE 


Figure 2.2-13. Changes in total ozone derived from Dobson data for the eight-year period over which 
SBUV measurements are available, November 1978-October 1986 for (a) winter and (b) summer. The 
plotting symbol is E = Europe, N = North America, and J = Japan. Results are expressed in percent per 
decade. 


running parallel to the country. Europe also shows a much flatter gradient in summer than in winter, and 
this is substantiated by SBUV which gives a uniform change of about — 4 % per decade over this region. 
Over North America SBUV indicates a very weak gradient, and Figure 2.2-13 (b) reveals ozone changes 
which depend on the specific location with no obvious latitudinal pattern. 

In view of the behavior shown in Figure 2.2-2 and in the SBUV data set, it is of interest to compute 
the changes in ozone after November 1978 contained in the Dobson record. Table 2.2-8 presents results 
from the AS2 model, expressed in percent per decade, based only on 8- and 10-year periods beginning in 
November 1978. Values appear for winter and summer at latitudes 55°N, 45°N, and 35°N and should be 
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Table 2.2.8. Changes in total ozone derived from the Dobson Network for periods coincident with 
the SBUV and TOMS data sets (changes in the 8 and 10 years data sets are expressed 
in equivalent % per decade to facilitate comparisons; one standard error appears in 


parentheses) 

Latitude 

Winter” 

8 Years 10 Years 

Summer” 

8 Years 

10 Years 

55 N 
45 N 

IS N 

-1.8(1. 5) -2.2(1. 1) 

-2.0 (1.0) — 2.7 (0.7) 

-2.2 (1.5) ~3.3 (1.0) 

-4.0 (0.8) 
-4.0 (0.8) 
-3.9 (0.8) 

-2.4 (0.7) 
- 1.6 (0.5) 
-0.8 (0.7) 

-The 8 and 10 year records cover November 1978-October 1986 and November 1978-October 1988, respect, vely. Due to the 
sets, the derived changes should no. be interpreted as long-term trends in the sense presented ,n 

previous tables. Winter = December-March, Summer - May-August. 


compared with results in Table 2.2-5. As expected, standard errors derived for the shorter data set are 
larger than those in Table 2.2-5, but the trends still remain significant at the one sigma level. With an 
adopted uncertainty of one standard error, the wintertime changes derived from the short and long-term 
Dobson data sets agree. However, a very different picture emerges in summer. Here the ozone changes 
expressed in percent per decade derived over a 10-year period are very different from those based on 8 
years of data. Furthermore, both sets of summer results in Table 2.2-8 imply larger changes m ozone per 
unit time than does Table 2.2-5. The difference between the 8 and 10-year results arises from the upturn in 
summetime ozone after 1986 as shown in Figure 2.2-2 (c). The results in Table 2.2-8 should 
as long-term trends, but they are nonetheless important. In particular, they show that after statistical 
correction for known sources of variability, atmospheric ozone still displays a temporal behavior whi i 
more complex than a long-term linear trend. This underscores the difficulty in interpreting derived trends 

in terms of cause and effect. 


2.2.3.3 Conclusion of Satellite Data Analyses 

The analyses described above yield the following conclusions. (1) For an 8-year record beginning in 
November 1978 the patterns of ozone change derived from SBUV and TOMS at the 14 selected sites are 
generally consistent. This suggests that SBUV and TOMS are compatible for detecting 8**^ paUerns 
in ozone change. In specific regions of very tight gradients, however, the high horizontal resolution of 
TOMS may provide a more realistic result. (2) Maps of changes in total ozone based on SBUV show 
significant variability around latitude circles, and a zonal average ozone change can 
difference between highly variable quantities. (3) Seasonal changes in total ozone observed by SBUV f 
November 1978 through October 1986 differ from expectations based solely on the long-term Dobson 
record. A separate analysis of Dobson measurements beginning in November 1978 supports this conclusion. 
However, the rate of change of summertime ozone derived from the short Dobson data set is very sensitive 
to the addition of measurements made after 1986. 


2.3 TRENDS IN VERTICAL OZONE DISTRIBUTION 


2.3.1 Introduction 

This section summarizes present understanding of trends in the vertical ozone distribution. Emphasis 
is given to the altitude range in the upper stratosphere where the percentage change in ozone concentration 
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at mid-latitude due to anthropogenic chlorine perturbation is expected to be the largest. The major problem 

in such an analysis is the scarcity in space and time of the data available from both satellite and ground- 
based measurements. 

As the Ozone Trends Panel (OTP) evaluation has shown that measurements performed by the SBU V 
instrument on the Nimbus-7 satellite cannot be used for trend analysis, the satellite section (2.3.2) is based 
only on the measurements of the SAGE I (February 1979-November 1981) and SAGE II (October 1984- 
December 1988) instruments. As tar as ground-based measurements are concerned, only 10 stations in the 
Northern Hemisphere have a long and reliable enough record of Dobson spectrophotometer Umkehr 
measurements to be used for trend analysis. In addition, such data are affected by the very high aerosol 
load in the stratosphere from 1982 to 1984 following the major eruption of the El Chichon volcano. Finally, 
when considering the available stations performing regular balloon-borne ozonesonde measurements, the 
situation is even worse with regard to data reliability, especially at the uppermost level. 

Therefore, the same merits and limitations of both orbiting remote sensors and the ground-based 
network, as emphasized in Section 2.2.1 for total column ozone trends analysis, are also relevant in this 
section. The limitations are even more stringent here, because the biases due to the limited geographic 
coverage and length of record may be larger due to the restricted data base which can be used. The 
objectives ot the present analysis are thus limited. They are: (I) to update the SAGE II data set through 
1988 as compared to the 1986 limitation of the Ozone Trends Panel Report, (2) to perform a detailed 
statistical analysis of both Umkehr and balloonsonde measurements to evaluate the usefulness of the 
available ground-based data for the detection of trends in the vertical ozone distribution throughout the 
stratosphere. As for the total column ozone study performed in the Ozone Trends Panel Report, a trend 
analysis has been made for each station in order to derive on a station-by-station basis the various natural 
(solar cycle, aerosols) and man-made variations in the ozone concentration at various altitude levels. It 
thus constitutes a completely new data evaluation. 


2.3.2 Comparison of SAGE I and SAGE II Stratospheric Ozone Measurements 

SAGE I and II use the solar occupation technique to measure aerosol extinction and trace gas 
concentration in the atmosphere (McCormick et al„ 1979; Mauldin et al„ 1985 a,b, see Section 2. 1. 1.2). 
The retrieved ozone profiles have a vertical resolution of 1 km, and their range extends from cloud top up 
to an altitude limited by the signal-to-noise ratio; for SAGE I the upper limit is about 55 km and for SAGE 
II it is approximately 65 km. Profiles are smoothed over a 5-km layer at heights above approximately 48 
km to increase s.gnal-to-noise ratio at these altitudes. SAGE I and II instruments utilize essentially the 
same design with similar optical components. Differences between the two instruments include the number 
o spectral channels (7 for SAGE II and 4 for SAGE I), the use of narrower spectral bandwidths on SAGE 
, and the use oi a rectangular field-ot-view on SAGE II versus a circular one on SAGE I These differences 

do not introduce any significant systematic errors between SAGE I and II ozone measurements as will be 
discussed below. 

The principal sources of error in the measurement of an individual ozone profile are radiometric 
imprecision, digitizer truncation, and scan mirror pointing errors. These random errors are approximately 
uncorrelated vertically, and their combined effect on an individual retrieved ozone profile is estimated at 
each point on the profile from the variance of the measurements from approximately 4-5 scans of the 
SAGE mirror across the viewing altitude. Each SAGE profile, however, also possesses an uncertainty in 
reference altitude ot approximately 0.25 km which contributes to uncertainties in profile repeatability 
(Cunnold et al., 1989). The principal systematic errors between SAGE 1 and SAGE II arise from instrument 
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scan mirror calibration and the uncertainty in our knowledge of the absorption cross-section of ozone at 
600 nm convolved with instrument bandpass and solar spectrum. Aerosols can also produce biases in the 
ozone retrievals, but only at those altitudes where aerosol concentrations are large (i.e., mostly below 25 
km altitude). The relative systematic error between the SAGE instruments is estimated to be approximately 
2% (Cunnold et al., 1989). The random error components for both instruments are less than 10% for a 
vertical resolution of 1 km. 

2.3.2.1 Zonal Means 

SAGE I and II scan the atmosphere over the Earth’s limb at the terminator twice on each of 15 daily 
orbits. The orbital inclination causes the observation latitude to periodically vary from approximately 75°S 
to 75°N dependent on season. The latitudes of the SAGE I and SAGE 11 measurement locations for 
February 1979 to November 1981, and October 1984 to November 1987, respectively, are displayed for 
sunset measurements in Figure 2.3-1. Due to the SAGE I orbital elements, the coverage does not repeat 
from one year to the next as it does for the SAGE II measurement location latitudes. Any gaps in the 
curves of Figure 2.3-1 are periods of time during which no satellite occultations occurred. Because of power 
problems on the spacecraft, only sunset data were taken by SAGE I after July 1979, and additional gaps 
are observed in the SAGE I curves due to this problem. A complete set of sunrise measurements for SAGE 
II, however, is available. 

For analysis purposes, the SAGE I and II data were spatially grouped into 1 1 latitude bands. The 
equatorial latitude band was 15°S to 15°N, and the midlatitude bands were 10° wide centered at 20°, 30°. 
40°, 50°, and 60° north and south. As the SAGE I or SAGE II observations passed into and subsequently 
out of one of the latitude bands above, the ozone measurements were averaged to form a zonal mean. In 
addition, the temporal duration defining a zonal mean was limited to a maximum of 1 week. Thus, a given 
latitude band contains more than one zonal mean if SAGE I or SAGE II sampled the latitude band tor a 
continuous period of greater than I week duration. These typically occur when the measurement location 
latitude changes from increasing with time to decreasing with time. 

Time series of SAGE I and II ozone zonal means at 30 km altitude are shown in Figure 2.3-2a for 
selected latitude bands. For comparison purposes, they have been time shifted to a common origin. Each 
time series was analyzed with a three-component linear regression model containing semiannual, annual, 
and 27-month period terms. The 27-month period term was chosen to model a quasi-biennial oscillation 
(QBO) in the ozone. The model fits to the data, which are intended as a visual aid, are shown as the 
continuous curves. At 30 km these time series display the expected features of seasonal ozone variability, 
namely the strong annual variation in the midlatitudes and the semiannual oscillation in the equatorial 
region. From Figure 2.3-2a, it is observed in the midlatitudes that both instruments agree in amplitude and 
phase with no apparent consistent bias. The largest differences occur at the Equator and are due to the 
fact that the QBO terms in the model are not in phase between the two data sets. 

The vertical distribution of the annual and semiannual component amplitudes between 25 km and 50 
km were computed and are shown in Figure 2.3-2b. The continuous curves on the plots represent the 
amplitudes in percentage from the model mean. Squares, either filled-in or open, mark those amplitudes 
which are statistically significant at the ±2 a confidence level. A comparison between the SAGE I and 
SAGE II annual and semiannual component amplitudes shows the similarity in the vertical distribution of 
amplitude. Both instruments agree on the altitudes at which the midlatitude annual amplitudes effectively 
vanish, i.e., at 25 km in both hemispheres and at 38 km in the Southern Hemisphere. At 25 km, a phase 
change occurs between the zonal mean time series below and above this altitude (Wang et al., 1989). 
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MONTH 

Figure 2.3-1. (a) Latitudes of the SAGE I sampling locations for sunset, (b) Latitudes of the SAGE II 
sampling locations for sunset. Because of the SAGE II orbital characteristics, the latitude coverage appears 
as a single curve even though 3 years are plotted. 


As of this writing, SAGE II has been in operation for over 4 years. The daily averages of ozone 
concentration sampled within 1 0°-w ide latitude bands symmetrically placed about the equator were modeled 
using the analysis previously described. The results for 35 km arc shown in Figure 2.3-3. The regressions 
adequately represent the seasonal and semiannual variability observed in the ozone time series. In addition, 
near the Equator, the QBO in ozone becomes clearly defined. The high vertical resolution of the SAGE II 
ozone measurements significantly improves our knowledge of vertical structure of the ozone QBO. Also 
evident in Figure 2.3-3. the semiannual oscillation appears as a significant component in the ozone vari- 
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■ SAGE I (a) ° SAGE 11 < b > 



Figure 2.3-2. (a) Zonal mean ozone number density at 30 km for SAGE I during the period February '1979 - 
November 1981 and SAGE II during the period October 1984 - November 1987 for the five latitude ba 
-roqrroc 25°S-35°S 15°S-15°N, 25°N-35°N, and 45°N-55°N. The vertical bars are the 95 /o confidence 
interval on the zonal mean. Smooth curves are the regression estimates. The data have been Time * sh'^ed 
for comparison as indicated on the abscissa, (b) Amplitudes of the seasonal and semiannual c W n ^ ts 
of the SAGE I and SAGE II zonal mean time series as derived from the regression coefficients. Units a e 

in percent relative to the model mean. 


ability, especially at 5°S. Marching poleward of the Equator, it is observed that the seasonal oscillation 
gains in amplitude with the semiannual oscillation decreasing in strength. 


2.3.2.2 Change in SAGE I and SAGE II Ozone Concentrations 

The change in stratospheric ozone concentration that occurred between the^ 33-month lifetime of 
SAGE 1 and the first 4 years of SAGE II was estimated in three ways over an altitude range from i 25 to . 50 
km, the region of the atmosphere essentially free of aerosols that could have any effect on the SAGE I 

ozone retrieval. 

The first method of estimating the ozone change, which minimizes the effects of seasonal and spatial 
sampling biases, is to select sets of profiles where SAGE I and II sampled the same month of the year 
within the same latitude band. Such a period of time is defined as an intersection and corresponds to 
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Figure 2 . 3 - 3 . Daily mean ozone number density at 35 km for SAGE II for twelve latitude bands, each 
10 wide-centered on the latitudes listed. The continuous curves represent regressions using a model 
containing components with periods of 30 months, 1 2 months, and 6 months. 


period of time in the year, of approximate duration I month, during which both SAGE I and II sampled 
the same latitude band. Note that an intersection is independent of the actual year during which measure- 
ments were made. The intersection is defined strictly in terms of the month within the year. Within the 1 1 
latitude bands considered, 105 intersections were isolated containing a total of approximately 7,300 SAGE 

I profiles and 25,000 SAGE II profiles. For each intersection, the SAGE I and II profiles were averaged 
separately. ' 

The second method of estimating the ozone change is to simply compute the zonal mean time series 
averages and form the percentage difference. This method has the possible disadvantage that seasonal 
sampling biases may be introduced. 

The third method of estimating the ozone change is to use the regression coefficients derived from 
the zonal mean time series and compute the mean value of the series over the time period of the data. This 
method’s disadvantage is the empirical nature of the model. 


Using all three methods described above, the percentage difference between the SAGE II mean and 
the SAGE I mean was computed at each altitude using the SAGE I mean as the reference. The intersection 
(method I) percentage differences and their corresponding standard errors are displayed in the upper six 
panels of Figure 2.3-4. For all latitude bands, excluding 60°N and 60°S, methods 2 and 3 gave very similar 
results to those of method I. Throughout the 1 1 latitude bands considered, except for 60°S, the ozone 
changes from SAGE I to SAGE II are within the range of -4% to +4%. The larger decreases observed in 
the Southern Hemisphere near 60°S are heavily driven by the intersection in the austral springtime and are 
likely due to vortex dynamics. 

Identical analysis (method 1) was performed using two wider latitude bands: 20°N to 50°N and 20°S 
to 50 S, areas over which there is relatively dense sampling by both instruments. The mean percentage 


224 




GLOBAL TRENDS 



-10 -5 0 5 10 

PERCENT 


Figure 2.3-4. Mean percentage difference between SAGE II and SAGE I versus altitude computed 
from the intersection of SAGE I and SAGE II. All intersections occurring between 15 N to 55 N o 
15°S to 55°S were combined into one sample to produce the results in the lower panel. The horizonta 
bars, placed at 5 km intervals, represent the standard error of the mean percentage difference. SAGE 
I is the reference in ail percentage calculations. 


differences appear in the lower panel of Figure 2.3-4. The agreement between the two hemispheres is 
apparent, and indicate for the 6-year period between 1980 (SAGE I) and 1986 (SAGE II): 

—an ozone decrease between 35 and 44 km with the maximum ozone change of - 3% ± 2 % occurring 

at 40 km, , 

—an ozone decrease of -3% ±2% at 25 km and an essentially zero ozone change at 28-33 km and 

45-48 km. 
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These values represent the changes that occur over this time period and no attempt has been made to 
separate out natural and man-induced variability. 

2.3.3 Analysis of Umkehr Data 

A trend analysis of stratospheric Umkehr profiles of ozone from 10 stations in the Northern Hemi- 
sphere is considered over the period 1977-1987. It constitutes an update of the analysis by Reinsel et al. 
(1987). The stations are Belsk, Arosa. Lisbon, Sapporo, Tateno, Kagoshima, Boulder, Edmonton, New 
Delhi, and Poona. Umkehr profile ozone data in the highest 5 Umkehr layers (4-9) are considered in the 
trend study, which covers an altitude range of approximately 19-45 km. A detailed analysis of the inter- 
dependence of the various Umkehr layers and resulting vertical resolution has been performed in the Ozone 
Trends Panel Report. It shows that ozone retrievals for Umkehr layers 4 to 8 are suitable for trend studies 
and that the effective vertical resolution for these layers varies from 15 to 1 1 km. Trends derived in layer 
9, where the retrieval is more sensitive to real ozone changes in layer 8— and has greater sensitivity to 
such changes than the layer 8 retrieval— are of limited interest. 

One of the major problems in Umkehr measurements is related to the impact of stratospheric aerosols 
from volcanic activities. This problem has become even more acute during the period following the eruption 
of the El Chichon volcano. A comprehensive theoretical examination of aerosol effects on Umkehr 
measurements was presented by Dave et al. (1979) and DeLuisi (1979) and the problem has been addressed 
by Reinsel et al. (1984) and DeLuisi et al. (I989a,b). The results of these studies indicate that accounting 
for the effects of aerosols on the Umkehr data measurements is critical in any trend analysis of these data 
The major volcanic eruption of El Chichon in April of 1982 produced significantly larger amounts of 
atmospheric aerosols than at any time in the previous 25 years. This, in turn, has produced a much more 
severe impact on Umkehr measurements from stratospheric aerosols from El Chichon than earlier periods. 
For illustration, deseasonalized monthly averages of Umkehr data in layers 4 through 9 at Arosa are plotted 
for the period 1977-1987 in Figure 2.3-5. The effects of aerosols on the Umkehr data, starting towards the 
end of 1982, are apparent. We note that low Umkehr values occur during 1982-1983 not only in the upper 
layers 7-8 but also in layers 4, 5, and 6. This is in contrast to the earlier theoretical calculations of Dave et 
al. (1979), which indicate that little error in Umkehr measurements should occur in the middle layers 4, 5, 
and 6 due to volcanic aerosols. This, in turn, suggests the possibility that real decreases in ozone occurred 
in these middle layers during 1982-1983 as a result of El Chichon and other natural events of that time 
period. Hence, large empirical adjustments of these data in layers 4, 5, and 6 for aerosol error effects during 
this period may not be appropriate. 

Because of this severe impact on the Umkehr data, considerable caution needs to be exercised when 
performing a trend analysis using Umkehr data over the period 1982-1987. In particular, the techniques of 
the previous analysis (Reinsel and Tiao, 1 987) need to be modified to account for the possibility of differences 

in the aerosols effects on Umkehr data during the El Chichon period relative to earlier effects especially 

any possible nonlinearity in the effects on Umkehr data as a function of optical thickness because of the 
much larger range of observed optical thickness values during El Chichon than in earlier periods. Because 
of this, in present trend analysis of both empirically corrected and theoretical model-corrected Umkehr 
data, two different procedures were considered to treat Umkehr data during a portion of 1982-1983. In the 
first one (i), all available Umkehr data were used. In the second one (ii), Umkehr data during the portions 
of 1982-1983 most affected by aerosols were omitted from the analysis. This resulted in the omission of, 
at most, 8 months of data at each station, generally for the period from November 1982 through June 1983 
Procedure (ii) is motivated by the fact that errors in Umkehr data appeared to follow a somewhat different 
relation to optical thickness for these extreme cases. Also in an empirical adjustment approach, these 
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Figure 2.3-5. Desea S onalizedUmkeh,dataatA,osafo,layers4.9.unadiusledioraerosols.forthepe.iod 

January 1977 through December 1987. 


extreme cases might exert a substantial influence on the overall estimation results for the aeroso effects 
and hence on estimation results for other parameters. In addition, it is rather clear that Ul " e r a * * 

this period would be the most sensitive and most difficult to accurately correct in the theoretical mo 

Ca,CU The trend analysis has been performed for each of the 10 stations usung two statistical models | which 
differ in the approaches of the aerosol effects. These models are described in Appendix 2. A. They include 

Hawaii (Reinsel et al 1984, 1987). More recently, lidar aerosol data from several locations in the north 
SKi inciuding Hampton. USA. Garmish. FRO. 

France became also available for the last several years since around 1977. It is expected that 
will be more representative of aerosols at the Umkehr station locations than the Mauna Loa data Th 
ul data have been used in ade.ailed study by DeLuisi e. al. (1989a) to obtain stra.osphene aemso^l 
thiekness derived from the lidar data for the period 1977-1987 as shown in Figure 2.3-6. DeLuis, 


227 








GLOBAL TRENDS 


OPTICAL THICKNESS (t) 



TRANSFORMED OPTICAL THICKNESS 



(1989a) also used the lidar data in a radiative transfer model to calculate theoretical model estimates of the 

d^du'rme^TH 'im ^d aerosols ’ and thus Provide theoretical corrections for the Umkehr 

data during 1978-1986. A detailed description of the correction methods are given by De Luisi et al. (1989a). 

Parameter values for both time series models (1) and (2) were estimated for each of the 10 Umkehr 
stations and each of the five Umkehr layers. 4 through 8, using both procedures (i) and (iVlndoned 
prevmusly (For comparison, a model of the form (2). which does not include any empirical aerosol 
adjustment, was also estimated for the uncorrected Umkehr data.) Therefore, two distinct sets of estimation 

Umkehr ^fu °u emp ' rical aeroso1 ad J^tment model (I) and the model (2), where 

et al 1989a) Th ^ f0r f roSo1 errors usi "g theoretical model-based calculations (DeLuisi 

al., 1989a). The results will be referred to as “empirical corrections” (model 1) and “theoretical 
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A. EMPIRICAL CORRECTIONS MODEL 



L = LISBON 


S = SAPPORO 
T = TATENO K = KAGOSHIMA 


B = BOULDER E = EDMONTON 
N = NEW DELHI P = POONA 

-rr 



7 lmk °L data trend estimates <% per year) from 10 individual Umkehr stations for layers 
4 9 tor £ pertod lnua “ 977 through December 1987. (Portions of the data, approximately from 
Noilmber 1M2 to June 1982 were deleted in the analysis.) (a) Trend estimates based on empirical 
aerosol corrections model, (b) Trend estimates based on theoretical aerosol corrections m 

.. r^j-t and f or the two procedures as “no deletions” (method i) and “with deletions 
^ethoIT ii) The individual station estimates of trend and solar cycle coefficients from both models are 
pTotted^ in Figure 2 3-7 a-b for each Umkehr layer, for the method (ii), when data were deleted during a 

portion of 1982-1983. 

Estimates of the overall trend in ozone for each layer, as well as reasonable estimates of the uncertainty 
in the overall estimates, have also been computed by using random-effects models for individual station 

p • ai M987) The resultant +2 cr confidence interval estimates of trends during 

estimates as in Reinsel et al. (198/). me resuuaru nrocedures 

1977-1987 from both the empirical and the theoretical adjustment models (1) and U) for both Proc 
are given in Table 2.3-1, together with the solar flux coefficient estimates in % per 100 F I0 , solar units. 

The results of the above analysis indicate statistically significant decreases in ozone during 1977-1987 
in Umkehr layers 7 and 8 using either method of correction for aerosol errors, and using either procedure 
,i) or (ii) for the treatment of data during a portion of 1982-1985. The delelion procedure yields slices 
negative trend estimates in layer 8 for Ihe empirical adjustment technique, and has little effect layer 
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Table 2.3.1. Overall trend estimates (in % per year) from 10 Umkehr Stations for the period 1977- 
1987 with associated standard errors (±2 a) 


Layer 

8 

7 

6 

5 

4 


Empirical Correction Model (1) 

No Deletions Deletions in 1982-83 

( -0.59± .27) (-0.51 ±.26) 

( -0.45 ± . 17) ( -0.45 ± . 16) 

( 0. 1 6 ± . 1 1 ) (-0.17zt.ll) 

( —0.08 ±.11) (— O.lOzt .11) 

(-0.08 ±.20) (-0.1 8 ±.20) 


Theoretical Correction Model (2) 

No Deletions Deletions in 1982-83 

(-0.43 ±.28) (-0.42 ±.27) 

( 0.43 ± . 19) ( -0.43 ± . 18) 

( — 0.26 ±.13) ( —0.22 ± . 1 1) 

( - 0.21 ± . 12 ) ( — 0 . 14 ± . 1 1 ) 

(- 0.21 ±. 20 ) ( - 0 . 15 ± . 20 ) 


8 

7 

6 

5 

4 


Overall Solar flux estimates (in % per 100 units of F 107 flux) from the period 1977-1987 


( 1.41 ±0.81) 
(1.25 ±0.89) 
(0.71 ±0.63) 
(0.73 ±0.61 ) 
(2.66 ±1.15) 


( 1 .48 ±0.81 ) 
(1.28 ±0.97) 
(0.81 ±0.68) 
(0.90 ±0.62) 
(2.87± 1.20) 


(1.47 ±0.94) 
(1.20 ±0.93) 
(0.38 ±0.74) 
(0.47 ±0.68) 
(2.49 ± 1.18) 


(1.61 ±0.85) 
(1.38 ±0.97) 
(0.72 ±0.66) 
(0.79±0.61 ) 
(2. 82 ± 1.15) 


for the theoretical model adjustment technique. The trend results for layer 7 are essentially the same for 
a methods. In the middle layers (4 and 5). the deletion procedure has more effect in the theoretical 
adjustment model, with less negative trend estimates occurring when data are deleted, but the reverse is 

tfU l a f mp,r,cal m « del - For the part, when the deletion procedure is used, the two correction 
methods yield fairly similar trend results overall. By contrast, the corresponding trend estimates obtained 
trom a model such as (2) but with uncorrected Umkehr data are -0.87 ±0.28, -0.71 ±0.19, -0 31±0 14 
-0.23 ±0. 14, and -0.26 ±0.22% per year, respectively, for layers 8 through 4. Hence it is seen that failure 
to incorporate corrections for aerosol errors in the Umkehr data leads to highly distorted trend results 
especially for layers 8 and 7. 

i T i“ e ll ?" d u eslimates obtained 'n this analysis for the period 1977-1987 are generally more negative in 
iq 5 ^? x. 8 than previous trend estimat es obtained using data for the period 1970-1981 (Reinsel et al. 
novo\ ^k^ 0 ' tha * these estima tes are substantially less negative than those obtained by DeLuisi et al 
1989a) which were based on a composite series of five Umkehr stations over the shorter period of 1979- 
1986 and which did not include an effect for solar flux variations. These trend results must be interpreted 
cautiously because of the shortness of the time period considered and because of the partial confounding 
between estimates of trend and solar flux effects over this relatively short period. As an illustration of the 
dependence of the trend results on the particular time period considered and on the effects of solar 
variations, a similar trend analysis for data over the period 1978-1987 was also performed. It yielded trend 
estimates in layers 4-6 that were generally positive and approximately 0.2% per year more positive than 
the trend estimates in Table 2.3-1 for those layers. 

Finally one should note that the trends reported here take into account the effects of the solar cycle 
They thus lead to an average decrease in ozone of about -4.8 ±3.1% at 40 km for the 10-year period 
e ween 1977 and 1987. When compared to the changes observed by the SAGE instruments, and within 

the uncertainty limits of SAGE and Umkehr trend results near 40 km, these two independent results are 
not inconsistent. 
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2.3.4 Analysis of Ozonesonde Data 

A trend analysis using ozonesonde data collected at 13 stations over the period June 1965 through 
December 1982 has been previously reported (Tiao et a.., .986). This study has been updated for 9 of the 
13 stations (Table 2 3-2) using all available data through December 1986. Individual ozonesonde readings 
in unt of partial pressure at various pressure levels were first integrated into ozone amounts in 15 
“fractional^ Umkehr layers, shown in Table 2.3-3. The readings were screened to meet the following 
criteria- (1) the correction factor was between 0.9 and 1.2 for the Brewer-Mast (BM) sondes and between 
0 9 and 1 15 for the ECC sondes, (2) the burst level was above the top of Umkehr layer 5B, (3) e 
corresponding daily total column ozone reading was available, and (4) readings with zero part.a pressure 
were removed In the earlier study, the correction factor in the screening procedure was set between 0.8 
Z 1 4 for the BM method and between 0.8 and 1.3 for the ECC method. The regression mode, used m 

the analysis is described in the appendix. 

Employing available data through 1986, Table 2.3-4 gives the trend estimates and their associated 
estimated 5 standard errors for each of the 15 layers of each of the 9 stations. The estimates are gi 
nercent change per year Figure 2.3-8 shows profiles of the trend estimates against the layers for ea 
station Broadly spelking over the majority of the stations, the noticeable features are a positive trend in 
^e lower layers 1A-1D, a negative trend in the middle layers from 2A to 4B, and little or no trend for the 
layers 5A and above. An overall trend profile across the nine stations can be constructed by considering 


Table 2.3.2. Ozonesonde stations, data spans 

, and measurement methods 


Qtatinn 

Latitude 

Data Span Method of Measurement 

Kagoshima 

Tateno 

Sapporo 

Payerne 

Hohenpeissenberg 

Goose 

Edmonton 

Churchill 

Resolute 

32N 

36N 

43N 

47 N 

48 N 
53N 

54N 

59N 

75N 

1/70-10/86 
3/68-1 1/86 
3/68-4/86 
9/68-12/85 
1/70-12/85 
6/69-8/80 
9/80-12/86 
10/70-8/79 
9/79-11/86 
10/73-8/79 
9/79-10/86 
3/66-11/79 
12/79-10/86 

Japan 

Japan 

Japan 

Brewer-Mast 

Brewer-Mast 

Brewer-Mast 

ECC 

Brewer-Mast 

ECC 

Brewer-Mast 

ECC 

Brewer-Mast 

ECC 


Table 2.3.3. Fractional Umkehr Layers 


Layer 

IA 

IB 

1C 

Upper Boundary 
(km) 

2.9 

5.5 

8.0 

Layer 

4A 

4B 

5 A 

Upper Boundary 
(km) 

21.3 

23.6 

25.8 


ID 

2A 

2B 

3A 

3B 

10.3 

12.5 

14.7 

16.9 

19.1 

5B 

6A 

6B 

7 


28.1 

30.5 

32.8 

35.1 
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Table 2.3-4. 


Trend estimates for nine stations (in % 
available) 


per year; data through December 1986 when 


Layer 

7 

6B 

6A 

5B 

5A 

4B 

4A 

3B 

3A 

2B 

2A 

ID 

1C 

IB 

IA 


Layer 

7 

6B 

6A 

5B 

5A 

4B 

4A 

3B 

3A 

2B 

2A 

ID 

IC 

IB 

1A 


Layer 

7 

6B 

6A 

5B 

5A 

4B 

4A 

3B 

3A 

2B 

2A 

ID 

IC 

IB 

IA 


Kagoshima 32°N 


Trend 

S.E. 

-0.0949 

0.2570 

0.1168 

0.2442 

0.4450 

0.2336 

0.1386 

0.1486 

0.2157 

0.1678 

-0.0326 

0.2033 

-0.9895 

0.1871 

-0.5869 

0.4582 

-0.6076 

0.5495 

-0.7148 

0.6638 

1.5634 

0.6002 

0.3251 

0.6275 

0.8453 

0.3820 

1.1980 

0.4418 

2.5065 

0.4815 

Payerne 47°N 

Trend 

S.E. 

-0.3927 

0.1189 

-0.2419 

0. 101 1 

0.0246 

0.0872 

-0.1226 

0.0708 

-0.2046 

0.0691 

-0.4088 

0.0680 

-0.5021 

0.0831 

-0.821 1 

0.1518 

-0.6697 

0.2427 

-0.9774 

0.2857 

-0.2042 

0.41 18 

0.3671 

0.8308 

1.6038 

0.4287 

2.1197 

0.3390 

2.4457 

0.4170 

Edmonton 54°N 

Trend 

S.E. 

3.0180 

0.3989 

2.9733 

0.3887 

2.2304 

0.3403 

1.0258 

0.2485 

0.2672 

0.2021 

-0.3388 

0.2735 

-0.2767 

0.2503 

-0.9817 

0.3416 

-0.4881 

0.5318 

-0.5528 

0.7432 

-2.1986 

1.1197 

- 1. 6091 

1.2641 

- 1.6674 

0.4212 

-0.7987 

0.5065 

-0.0882 

0.6898 


Tateno 36°N 
Trend S.E. 

-0.0288 0.2271 

0.0820 0.2188 

-0.0172 0.1501 

0.0939 0.1213 

-0.0913 0.1057 

0.0892 0.1281 

-0.3791 0.1656 

0.0672 0.1659 

0.3598 0.3561 

0.3021 0.3872 

0.4971 0.4650 

0.6034 0.4874 

0.3247 0.2446 

0.5309 0.1842 

1.2485 0.2877 


Hohenpei 

ssenberg 48°N 

Trend 

S.E. 

-0.4774 

0.1266 

-0.4549 

0.1129 

-0.2564 

0.0898 

0.0377 

0.0671 

-0.0928 

0.0731 

-0.2457 

0.0749 

-0.6372 

0.0832 

-0.5486 

0.1399 

-0.3854 

0.1817 

-0.4241 

0.2292 

0.0130 

0.3712 

2.3613 

0.3648 

2.6830 

0.2150 

2.2850 

0.1762 

2.3078 

0.2487 

Churchill 59°N 

Trend 

S.E. 

3.3806 

0.7391 

2.7236 

0.6748 

1.0688 

0.4455 

-0.1247 

0.3595 

-0.9767 

0.3092 

- 1.4080 

0.2847 

-1.9890 

0.0995 

-0.9869 

0.4493 

-1.1097 

0.5889 

0.0407 

0.3122 

-0.8426 

0.8610 

-0.3072 

1 .4824 

-1.0678 

0.6925 

- 1.3955 

0.5185 

-0.7152 

0.6638 


Sapporo 43°N 


Trend 

S.E. 

-0.2936 

0.3331 

-0.3130 

0.3102 

-0.1553 

0.2239 

0.0103 

0.1554 

0.2996 

0.1272 

0.0902 

0.1821 

0.2923 

0.2171 

0.8617 

0.3700 

0.7556 

0.5651 

1.5540 

0.6151 

2.0444 

0.6631 

3.0894 

0.6804 

-0.2249 

0.2728 

0.2234 

0.2888 

0.9364 

0.5446 

Goose Bay 53°N 

Trend 

S.E. 

1.1010 

0.3989 

1.3753 

0.3700 

1.1855 

0.2995 

0.4013 

0.2530 

0.0728 

0.2075 

-0.0806 

0.2150 

-0.6877 

0.2600 

-0.9480 

0.3735 

-0.8932 

0.5202 

- 1.5639 

0.4874 

-1.5020 

0.8146 

-3.1778 

0.9368 

-1.6966 

0.5074 

-0.9935 

0.4527 

- 1.0987 

0.4945 

Resolute 75 

°N 

Trend 

S.E. 

-0.2466 

0.3452 

-0.3685 

0.3263 

-0.3610 

0.1761 

-0.0532 

0.2658 

-0.1598 

0.0977 

-0.2328 

0.0769 

-0.6220 

0.0658 

-0.8826 

0.1834 

-0.4674 

0.1976 

-0.3612 

0.2293 

-0.5911 

0.4114 

1.7975 

0.0067 

1.4531 

0.6260 

0.6997 

0.4748 

-0.0624 

0.0204 
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Figure 2.3-8. Trend estimates derived from ozonesondes measurements at 9 stations for 15 fractional 
Umkehr layers (% per year). 


the nine estimates for each layer as a random sample from a normal population and compute the sample 
average and its estimated standard error. The results are given in Table 2.3-5. Statistically significant trend 
is observed in layers 3 and 4 (17-24 km). It corresponds to an average ozone concentration decrease ot 

0.46% per year. 

One should be extremely cautious to draw conclusions on a global scale based on the analysis of the 
ozonesonde data. Figure 2.3-8 shows that considerable differences exist in the trend profiles across the 
stations, especially in the lower (1 and 2) and the upper (6 and 7) layers. This is also reflected by the wi 
of the confidence intervals in Table 2.3-5. For example, the tropospheric trends derived for P&yerne > are 
probably too high because of an improvement in pre-launch procedure beginning in late 1983, not reflected 
by a level shift term in the statistical analysis. Otherwise, the two European m.d-lat.tude stations do show 
a rather coherent behavior with similar trends in the lower stratosphere. The situation is rather confusing 
in the troposphere: although the three Japanese stations have generally positive trends the three mid- 
latitude Canadian stations have negative trends, and large positive trends are observed at Resolute i <75 5 N) 
In addition, while statistically significant ozone decreases are observed in layers 3 and 4, consideration o 
Table 2 3-4 shows that the decrease in layer 4 is mainly due to the period 1982-1986, as no significant trend 
is observed when the data through 1982 only are considered (see Table 2.3-5b). Such a picture is rather 
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Table 2.3-5. 


Layer 

7 

6B 

6A 

5B 

5A 

4B 

4A 

3B 

3A 

2B 

2A 

ID 

1C 

IB 

IA 


Summary of averages of trend estimates (% per year) and standard errors (S.E. ±2cr) 


(a) 

Data thru 12/86 
Average 


0.6620 

0.6547 

0.4627 

0.1480 

-0.0744 

-0.2853 

-0.6434 

-0.5363 

-0.3895 

-0.2997 

-0.1356 

0.3833 

0.2504 

0.4301 

0.8312 


S.E. 

0.5041 

0.4536 

0.2890 

0.1224 

0.1292 

0.1523 

0.2052 

0.2073 

0.1967 

0.2941 

0.4549 

0.6499 

0.5135 

0.4391 

0.4646 


(b) 

Data thru 12/82 
Average 


0.7875 

0.6712 

0.3669 

0.1368 

0.0208 

0.0164 

0.1989 

-0.1615 

-0.4970 

-0.3090 

-0.0158 

0.3940 

0.1978 

0.5121 

0.6670 


S.E. 

0.6482 

0.4933 

0.4490 

0.2132 

0.1331 

0.0973 

0.5425 

0.1851 

0.1998 

0.4330 

0.6358 

1.1613 

0.4633 

0.2638 

0.4538 


consistent with the SAGE observat.ons as it might expand the ozone decrease observed between 28 and 
25 km towards lower altitude. However, in the SAGE data the maximum effect is observed at low latitudes 

™7 th Th S ° da f‘ Where ’ CXCCPt f0r Kagoshima ’ the decrease occurred at a number of mid- 
latitude stations Thus, apart from the decrease in layer 3 observed consistently over most of the stations 

it seems rather difficult to draw any definitive conclusion on a semi-global scale from this data set. 


2.3.5 Comparisons with Modal Calculations 

Adding 15 months ot new SAGE data to the analysis performed in the Ozone Trends Panel Report 
has not changed appreciably the observed change of -3 ±2% in ozone concentration near 40 km over the 
6-year period 1980-1986. The more thorough analysis of data from 10 Umkehr stations in the Northern 

3oTnd P 43 e r V h P 1977 t0 1987 reP ° rtS alS ° 3 StatiStiCa " y Significant decrease in ozone between 
30 and 43 km, the decrease near 40 km being 4.8 ± 3. 1% after allowing for seasonal and solar-cycle effects 

and correcting the data for aerosol interferences. The mean values of the observed decreases at 40 km are 

less than the mean values predicted by theoretical calculations performed in the OTP and taking into 

account both chlorine and solar influences. These are in the range of -6 to -8% over the period 1979- 

1985. However the large experimental uncertainties prevent a critical test of the theory. 

IImk i I he . SA f E i han f - °I ~ 3 ±2% at 25 km over 6 y ears when considered together with the negative 
Umkehr trend in this height range and the statistically significant, but not globally representative trend of 

nearly 0.5% per year observed by ozonesondes between 17 and 24 km, seems to suggest that the 
stratosphere below .5 km is the prime contributor to the total ozone loss at mid-latitude. Changes below 
25 km are not predicted by global models based only on gas phase processes. It is not clear whether this 
points to missing processes in the models of the lower stratosphere or that the sparse measurements are 
not representative of the global atmosphere. 
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2.4 TRENDS IN STRATOSPHERIC TEMPERATURE 

2.4.1 Introduction 

The purpose of this section is to reiterate and update the Ozone Trends Panel Report (OTP) temper- 
ature trend analysis. The objective of the OTP analysis was threefold: (I) to establish the quality of a 
variety of stratospheric temperature data sets; (2) to discuss temperature trend mechanisms, and (3) to 
determine the consistency of temperature changes with observed ozone changes. This last objective is 
particularly difficult to reach, since it requires data analyses that minimize the influences of other processes 
contributing to temperature trends. 

Global temperature data sets for the region above 30 hPa (—22 km) are only available after late 1978. 
For trend analysis, any data set covering such a relatively short period can be strongly influenced by short 
period anomalies and multi-year cycles such as the Quasi-Biennial Oscillation (QBO), the solar cycle, and 
even significant interannual variability. In particular, the OTP noted the strong thermal perturbation 
associated with the El Chichon eruption which occurred near the middle of the analysis period (1979-1986). 
Thus, rather than trying to estimate a trend, the differences between annual averages at the beginning and 
the end of the data period were calculated in the OTP. The annual averaging minimizes the QBO effect, 
and the calculation of changes minimizes the thermal perturbation associated with El Chichon, since it 
avoids using explicitly the years 1981 through 1984. The effect of the change in the solar flux over this 
period was estimated using model simulations. Finally, the OTP tried to minimize interannual variability 
associated with extratropical dynamics by focusing on temperature trends near the tropics to estimate the 
consistency of ozone and temperature changes. 

The update of the previous temperature analysis performed in the OTP report, is based on the 
temperature data set already available and updated when possible, and the recently available data from 
Rayleigh lidar observations. The possible implication for trend determination of the recent finding of a 
statistical relationship between the stratospheric temperature, the QBO, and the 1 1-year solar cycle is also 
discussed. 

2.4.2 Stratospheric Temperature Data Sets 

The stratospheric temperature data sets used in the OTP report were obtained from radiosondes, 
rocketsondes and satellite radiances. Recently lidar temperature data have become available. In this section, 
the three data sets used in OTP are briefly summarized and the new lidar data are discussed. 

2.4.2. 1 Radiosondes 

Radiosondes are balloon-borne instruments launched frequently at regular intervals. The density of 
radiosonde sites is relatively high over most continental regions. However, few measurements are available 
over the oceans, particularly in the Southern Hemisphere. Radiosondes can measure temperatures above 
10 hPa (30 km), but due primarily to the inaccuracy in the pressure measuring, the uncertainty in the 
temperature measurement above 30 hPa is large. Additional uncertainties arise from the lack of a radiosonde 
standard, and the variety of instrument types in use worldwide from different manufacturers. Also, for 
trend analysis, problems arise from the lack of documentation on instrumentation changes which have 
been introduced over the years, and application of corrections for solar short wave and long wave radiation 
errors. Although a number of authors report radiosonde temperature trends above 30 hPa, the OTP felt 
that the uncertainties of the measurements were too large above 30 hPa for trend determination. 
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2A.2.2 Rocketsondes 

Rocketsondes are rocket-borne instruments, usually thermistors, which can be carried into the 
mesosphere and lower thermosphere. Due to their larger cost and complex ground station requirements, 
rocketsonde stations are rather few and the number of sites has significantly decreased in the last few 
years. Temperature trends obtained from rocketsondes are also influenced by changes in instrument design 
in the early 1970s, nonstandard data recording procedures prior to 1969, and the estimate of aerodynamic 
and radiation corrections. The OTP used only high quality rocketsonde data subset accounting for problems, 
and eliminated rocketsonde measurements that presented temperature differences of 3 K or more with 
simultaneously launched radiosondes. 


2.4.2.3 Satellites 

Beginning with the launch of the TOVS (TIROS-N Operational Vertical Sounder) instrument aboard 
the NOAA polar orbiter series in late 1978, regular temperature soundings of the stratosphere became 
available on an operational basis (the Stratospheric Sounding Unit [SSU] and Microwave Sounding Unit 
[MSU] instruments). Satellite temperature measurements prior to this period lack sufficient accuracy or 
record length. 

Trend determination from radiance measurements from the TOVS series is complicated by calibration 
drift and differences between the performance of instruments in the series. Two approaches for the 
determination of satellite temperature trends were discussed by OTP. National Meteorological Center 
[NMC] analyses use temperature retrievals from TOVS radiances and Gelman et al. (1986, 1988) estimated 
corrections to the operational analysis temperatures for trends calculations. These corrections were derived 
from comparison with rocketsonde temperatures. Nash and Forrester (1986) also estimated temperature 
trends from TOVS radiances by intercalibrating the different TOVS instruments using the operational 
overlap of NOAA-6’s SSU and MSU. Thus they were able to estimate the change in the measured radiances 
in the subsequent SSU and MSU instruments and establish the stability of the measurements over the 
1980-1986 period. 

2.4.2.4 Lidar 

Rayleigh lidars provide density, and therefore, temperature information in the stratosphere and 
mesosphere. The Rayleigh lidar operates by detecting the Rayleigh backscattered signal from ground-based 
laser (for a description of the method see Hauchecorne and Chanin [1981] and Chanin and Hauchecorne 
[1984]). Due to the recent development of the lidar technique and the limited number of operational sites, 
data from this source were not used in the previous OTP report. To this date, two Rayleigh lidars located 
in France have been used on an operational basis: Observatoire de Haute-Provence (OHP) (44 N, 6 E) 
since 1981 and Biscarosse (44 N, I W) since 1986. The total number of profiles averaged over a nighttime 
period and defined hereafter as a sequence, exceeds 1000, with an average of two sequences a week. Since 
the lidar profile is determined by averaging several hours of data, the variability due to gravity waves, 
which appears in the instantaneous soundings, is reduced. 

The Rayleigh lidar temperature precision is height dependent and better than 1 K below 50 Km. 
However the accuracy of the temperature measurement may be affected by a number of factors: the 
presence of aerosols, the saturation of the pulse counters from low-altitude returns, and parallax effects 
between the emission and reception axis. These uncertainties can however be partly corrected. Aerosols 
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can be determined by simultaneous measurements at different wavelengths and a second channel with 
reduced efficiency can be used for the low-altitudes returns. Parallax effects can be avoided by using a co- 
axial system. To estimate the accuracy of the Rayleigh lidar, intercomparisons have been performed with 
rocketsondes. In 1986, rocketsonde and lidar measurements showed temperature differences between 1 K 
and 2 K from 30 to 70 km. This difference was smaller than the difference between the two thermistors on 
the same payload. Furthermore, intercomparison of the OHP and Biscarosse systems (located within 500 
km) over 66 simultaneous profiles showed a difference between the two systems of less than 2 K below 80 
km. 


2.4.3 Intercomparison of Stratospheric Temperature Data Sets 
2.4.3.1 Radiosonde, Rocketsonde, and Satellite Data 

The temperatures determined independently by Nash and Forrester (1986), NMC, radiosondes (Labitzke 
et al„ 1985), and rocketsondes (Angell and Korshover, 1983) were compared in the OTP report. Generally 
the analyses agreed, although there were periods of significant disagreement especially in 1982-1983. Figure 
2.4-1 (a,b,c,d) from OTP shows the results of these comparisons in the 100-30 hPa (16-22 km), 30-10 hPa 
(22-30 km), 5-1 hPa (36-48 km), and 2 hPa level (40 km). 


2.4.3.2 Comparison of Lidar Data with SSU and NMC 

A comparison was performed for the period 1981-1987 between the temperature obtained by lidar 
and the SSU and the NMC data. The data have been selected to be within 5 degrees of latitude and 10 
degrees of longitude of the OHP lidar site. It is however also important to note that both the SSU and 
NMC data used here are different from those used in the OTP report. The SSU data used in the OTP report 
were zonal means corrected for instrument changes as described in OTP Section 6.2.2. 1 (see Nash, 1988). 
The SSU data used here are the standard gridded radiances provided by the British Meteorological Office 
and have not been corrected, and the NMC data used have been adjusted as described by Gelman et al. 
(1986). 

Figure 2.4-2 shows the brightness temperature for the SSU channel 27 (approximately 1.5 hPa) and 
the lidar temperature convolved with the SSU weighting function given in Barnett and Corney (1984). 
Figure 2.4-3 shows over the same period the NMC data at 42 km (2.0 hPa) and the lidar data. The average 
lidar-SSU difference is of 4 K, with a maximum difference of 6.5 K. The average lidar-NMC difference is 
3 K at 42 km and varies from - 5 K to 7 K. 

The lidar-SSU and lidar-NMC differences are generally similar. They can be divided into four periods, 
which correspond to changes in the biases between the data sets. The first period beginning in 1981 and 
extending to early 1982 is characterized by larger values of the lidar temperatures. The second period, 
from early 1982 to mid-1983, is characterized by a close match up of the lidar SSU and NMC data. The 
third period, from mid-1983 to mid-1985, is characterized again by larger values of the lidar temperatures. 
The last period, extending from mid-1985 onward, shows smaller differences. Although the reasons for 
these differences are not readily apparent, it seems that the September 1983 satellite change from NOAA- 
7 to NOAA-8 resulted in the introduction of large bias between the lidar and SSU/NMC data sets. This 
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Figure 2.4-1 a. Seasonal temperatures with the long-term seasonal averages removed and smoothed 1 -2-1 in 
time: Angell 100-30 mb radiosonde thickness (solid dot); Berlin 100-30 mb thickness (long dash); NMC 70-30 
mb thickness (short dash); MSU channel 24 centered at approximately 90 mb (solid); and SSU channel 15X 
centered at approximately 50 mb (solid dot) for (a) Northern Hemisphere average, (b) Southern Hemisphere 
average, (c) global average, and (d) 30°S-30°N average. Berlin data were available only for the Northern 
Hemisphere. Tick marks on the abscissa correspond to the D-J-F seasons. 


effect is illustrated by the unadjusted data in OTP Figures 6.2. 2. 3 and 6.2. 2.4, where the change from 
NOAA-7 to NOAA-8 is approximately accompanied by a 4-K cooling. 

Comparisons were also performed between the lidar at 36 km and the SSU channel 26, and the lidar 
with NMC temperatures at 30 km (10 hPa), 36 km (5hPa), 48 km (1 hPa) and 55km (0.4 hPa) although this 
comparison relies on the downward extension of the lidar data by radiosondes, and is less germane to lidar- 
satellite comparisons. The best agreement between the lidar and NMC is obtained at 48 km with a maximum 
difference of 3 K. The causes of these differences will have to be studied further. 
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Fiaure 2 4-1 b. Seasonal temperatures with the long-term seasonal averages removed and smoothed 
1-2-1 in time: NMC 30-10 mb thickness (dotted); SSU channel 26X centered at approximately 20 mb (solid); 
and SSU channel 25 centered at approximately 1 5 mb (three short dashes and a long dash) for (a) Northern 
Hemisphere average, (b) Southern Hemisphere average, (c) global average, and (d) 30 S-30 N average. 
Tick marks on the abscissa correspond to the D-J-F seasons. 


2.4.4 Influence of Solar Activity 

Since the OTP report, a statistical relationship has been established between the temperature, the 1 1- 
year solar cycle, and the equatorial QBO (Labitzke, 1987; Labitzke and Van Loon, 1988). The relationship 
has first been observed in Arctic polar temperatures, when the temperatures are classified between years 
according to the west and east phases of the QBO. During the west phase, the polar temperatures are 
correlated with the solar cycle, while in the east phase, the polar temperatures are anti-correlated^The 
amplitude of this 1 1-year polar temperature cycle is approximately 10-15 K for either phase of the Q 
The relationship is illustrated in Figure 2.4-4 for the North Pole (from Labitzke and Van Loon 1988). It 
has been since extended to other latitudes and altitudes (Labitzke and Chan.n, 1988), and to the Southern 
Hemisphere. At present, no physical mechanism has been found to explain this relationship. 
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NMC (10-5 mb) 

SSU CH. 26 (Centered ~6 mb) 


Ftgure 2.4-1 c. Seasonal temperatures with the long-term seasonal averages removed and smoothed 
- -1 in time: NMC 1 0-5 mb thickness (dotted) ; and SSU channel 26 centered at approximately 6 mb (solid) 
for (a) Northern Hemisphere average, (b) Southern Hemisphere average, (c) global average, and (d) 30°S- 
30 N average. Tick marks on the abscissa correspond to the D-J-F seasons. 


The main consequence of the QBO-solar cycle relationship relevant to the present analysis is that it 
may be a source of long-term interannual variability. In addition, this relationship is variable with season, 
latitude, longitude, and height and generally, local data show a significantly larger QBO-solar cycle signal 
than zonal mean data. One should then be cautious when using local or regional data, unless the record 
length is large enough to allow the separation of trend and solar induced effects. 

2.4.5 Updating of Previously Used Data Sets 

Updates of global stratospheric temperature trends have been performed from radiosondes data by 
Angell ( 1988) and Labitzke ( 1989, private communication). Updates of the rocketsonde data have also been 
provided by Kokin et al. (1989) for the USSR stations of Heiss Island, Volgograd, Thumba, and Molod- 
ezhnaya and by the Japan Meteorological Agency for the Ryori station. 
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------ NMC (5-1 mb) 

- SSU CH. 27 (Centered 2 mb) 

SSU CH. 36X (Centered ~1.5 mb) 

Figure 2 4-1 d. Seasonal temperatures with the long-term seasonal averages removed and smoothed 
1-2-1 in time: NMC 5-1 mb thickness (dotted); SSU channel 27 centered at approximately 1.5 mb (solid 
dot) for (a) Northern Hemisphere average, (b) Southern Hemisphere average, (c) global average, and (d) 
30°S-30°N average. Tick marks on the abscissa correspond to the D-J-F seasons. 


Angell (1988) derives a global cooling of -0.62 K/decadc in the 100-50 hPa layer ovci the peiiod 
1973-1987, which mostly results from a strong cooling (-2 K in the Antarctic region). Both Labitzke 
(1989, private communication) and Angell ( 1988) conclude on a barely significant warming in the temperate 
Northern Hemisphere (0.2 K/decade at 50— 60°N from Labitzke) at the same altitude level. 


At lower latitudes (I0-30°N) at the 30 hPa pressure level, Labitzke (1989, private communication) 
reported no trend in the summer between 1961 and 1981 (excluding the periods during which volcanic 
aerosols may have warmed the lower stratosphere following the eruptions of Agung 1 1963] and El Chichon 
[1982]). Angell (1988) found a0.4-K trend for the Western Hemisphere tropics. On the other hand, a mean 
decrease in temperature of 0.4 K is reported by Labitzke ( 1989, private communication) at 50-60°N, and 
a large decrease of 2 K at the polar stations is reported by Kokin ct al. (1989) for the last 20 years. 
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Figure 2.4.2. Lidar temperature and SSU channel 27 brightness temperature from 1 981 to 1 987 (with a 
running mean of 90 days). 
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Figure 2.4-3. Lidar and NMC temperature at 42 km (with a running mean of 90 days). The periods 
corresponding to different NMC analyses are indicated at the bottom. 


At the 10-hPa (30 km) pressure level. Labit/.ke (1989. private communication) reported a -0.5 K/ 
decade trend for the latitude region between 20-5()°N over the period 1964-1987. Over the same period, 
Angell (1988) reported a trend of -0.2 K/decade based on radiosonde data. At the same altitude level, the 
updated rocket data of USSR and Japan led both to a larger trend of 1 .5 K7 decade, although with a large 
spatial variability. Using also rockelsonde data, Angell ( 1988) calculated a trend of -0.8 K/decade at 10 
hPa in the Western Hemisphere tropics (1973-1987). and I K/decade in the upper stratosphere (5-0.5 
hPa). Although the difference between the 10 hPa pressure level radiosonde and the rocketsonde-based 
trends is within a ±2 ir confidence intervals (Angell, 1988), the OTP reported serious disagreements 
between edited (see OTP Section 6.2.3) and unedited rockctsondc data. Since the Angell (1988), USSR, 
and Japan results ate ftom unedited rocketsonde data, caution should be exercised when comparing those 
data to the rocketsonde data of the OTP. Furthermore, trends at these individual stations might be mainly 
representative of local evolution. 
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SOLAR FLUX AND 30-mbar NORTH POLE 


(JAN + FEB) 12 TEMPERATURE 



1956 1960 1970 1960 1990 


Figure 2.4-4. Comparison of the 10.7 cm solar flux, in 10 22 W/(m 2 Hz) (solid line) and the 30 mb 
temperature, in degrees Centigrade, at the North Pole, averaged for January and February for (a) all 
winters, (b) winters in the West phase of the quasi-biennial oscillation (QBO), (c) winters in the East 
phase of the QBO. Asterisks denote major mid-winter stratospheric warmings (from Labitzke and 
van Loon, 1988). 

As a common conclusion, these studies give some evidence of a generally negative trend, with large 
spatial variability. The upper limit of the zonally averaged trend is between -0.6 and -0.8 K/decade 
between the 100 hPa and 10 hPa pressure levels (16-30 km) respectively, with the exception of the Antarctic 
region. These values are much lower than the local trends derived from rocketsonde data. 

2.4.6 Comparison with Trends in Ozone Distribution 

The reported SAGE 1/SAGE 11 changes between the 1979-1981 and 1984-1988 periods, in the 20-50 
degree latitude bands lead to a maximum decrease of 3% ±2% in ozone concentration at 40 km and a 
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decrease of - 5% at 50 km in the tropics (see Section 2-3-2). These values are very similar to those already 
reported in the OTP. where the SAGE data used were the zonal mean profiles for the latitude band I5°S to 
I5°N which gave an ozone concentration decrease of - 5 % peaking at 50 km. Because the predicted 
temperature change is sensitive to the ozone change above 50 km, two cases were considered: ( 1 ) no ozone 
change outside the 26-50 km altitude range were SAGE measurements are available, (2) a -5% ozone 
change above 50 km and a -2% change below 26 km. The radiative response showed a temperature 
decrease ol I K to 1 .4 K at 50 km. An additional 0.8-K cooling will occur at the same altitude due to the 
decrease in the solar flux during the declining phase of the solar cycle. Thus, the expected 50 km temperature 
was estimated to be of 1.8 to 2.2 K in the tropics. With the updated value of the tropical ozone decrease 
of 4%, the 50 km change should be 1.4 to 1.8 K. 

Changes on global and equatorial temperatures over the period from 1979-1980 to 1984-1985 are 
shown in Figures 2.4-5 and 2.4-6, reproduced from OTP Figures 6.3.2. 1 and 6. 3. 2. 2. The decrease in 
equatorial temperatures at 48 km from satellite and rocket data over the same period is about 1.75 K ± 

1 K. Global temperatures show a smaller decrease of 1.5 K ±1 K. The SAGE ozone changes and the 
observed equatorial temperature changes appear to be in satisfactory agreement within the ± 2 o confidence 
interval of the analyses. 

2.5 TRENDS IN TROPOSPHERIC GASES AND OZONE 
2.5.1 Introduction 

The Ozone Trends Panel (OTP) Report presented information on the global trends of a number of 
tropospheric trace gases that interact with stratospheric ozone, either chemically, as sources of nitrogen. 
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Figure 2.4-5. Summary of global stratospheric temperature differences ( 1 985 and 1 986 minus the average 
of 1979 and 1980). Symbols denote Angell’s radiosonde data (AR), NMC data (NMC), and satellite data 
(channel numbers, see Figure 6.2.2. 1 of the OTP Report). Error bars denote uncertainties (see text). Nash 
and Forrester (1988) estimate their satellite errors as 0.2 K rms (1 K rms for synthesized channels). 
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Figure 2.4-6. As in Figure 2.4-4, except to, tropics and induing racket ^ns^wsjatein Island (Kl; 
9°N, 168°E), Ascension Island (Al; 9°S, 14°W), and Cape Canaveral (CC, 28 N, 81 )- 

hydrogen and halogen catalysts that substantially control ozone dest. 

stratospheric temperatures. The OTP Report reviewed available data up to the end of P P 

of .his section is to update the OTP Report with data, where possible, from 1987 and 1988. Section i 0.4 
describes the trends in halocarbons and reviews the progress in the area ot trace gas ca t ration, s s ow 
in the OTP Report, the calibration uncertainties associated 

and trends observed in atmospheric nitrous oxide 

(N > 0) XI methane (CH 4 ), and carbon monoxide (CO), as these species are of particular importance in 
determining future stratospheric ozone levels. It would appear in particular that the ^rrent rate of^ ^ 

date Iva aNe on urSce and free tropospheric ozone measurements at a number of local, ons are reviewed 
Sec^ In 2 5 7 in fte igh. of recent slggestions of large-scale tropospheric increases m the Northern 
Henhsphere. As an overview. the current trends, concentrattons. and lifehmes of trace gases are snmma- 

rized in Table 2.5- 1 . 


2.5.2 Halocarbons 

A variety of anthropogenic and natural halocarbons are present in the global atmosphere. The 
abundance sources, and sinks of this group of chemicals have been recent,, d.scussed ,n Gammon e. al. 
(1986) and in the OTP Report (Ehhalt et al., 1989). 
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Table 2.5-1. 


Concentrations and global trends of tropospheric gases for 1987. Lifetimes 
where available (adapted from Ehhalt et al., 1989) 


are given, 
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The halocarbons are the most rapidly increasing species on a percentage basis in the global atmosphere 
over the last several decades. There is concern about this class of compound because they are sources of 
stratospheric halogen-free radicals, which catalytically destroy ozone (Rowland and Molina, 1975) and 
because they absorb in otherwise low opacity infrared wavelength regions of the atmospheric absorption 
spectrum and thus contribute significantly to the greenhouse eflect (Ramanathan et al., 1985). For halo- 
carbons, whose tropospheric lifetimes are much longer than global circulation times, the trends measured 
at ground level, at sites remote from halocarbon sources, can provide a reasonably accurate determination 
of global atmospheric trends (Prinn et al., 1983). The global rate of accumulation of halocarbons in the 
background atmosphere have been recently reviewed (Prinn, 1988a; Ehhalt et al., 1989). 


2.5.2.1 CCI 3 F, CCI 2 F 2 , CH 3 CCI 3) CCI 4 

The global rates of increase of these compounds in 1987 (Table 2.5-1) are similar to those given in the 
OTP Report. Sources of new observations up to the end of 1987 arc the GMCC (Geophysical Monitoring 
for Climatic Change) program with stations at Barrow, Alaska (71°N); Niwot Ridge, Colorado (40°N); 
Mauna Loa, Hawaii (20°N); Cape Matatula, Samoa (14°S), and at the South Pole (CC1,F, CC1 2 F 2 : Bodhaine 
and Rosson, 1988). An update of the observations at Cape Point, South Africa (34°S) (CCI 4 F, CCI 4 : Scheel 
et al., 1989) and Hokkaido, Japan (45°N) (CCFF, CCI 2 F 2 , CH 4 CCI 3 and CC1 4 : Makide et al., 1987) is also 
included. 

Global data, additional to that presented in the OTP Report (up to mid- 1986), are not available from 
the GAGE (Global Atmospheric Gases Experiment) program due to an ongoing calibration re-evaluation. 
Possible revisions to previously presented data (post mid-1985) (Ehhalt et al., 1989) are believed to be small 
(Prinn, 1989, private communication). The GAGE data collected at Cape Grim, Tasmania, up to the end 
of 1987 are shown in Figure 2.5-1. 

The rate of increase of CH,CCI, at Hokkaido (3.7 ±0.8 pptv per year [ ±2 ct as for all confidence 
intervals given hereafter]; Makide et al., 1987) is significantly lower than the global rate of 6.2 ±0.2 pptv 
per year deduced from the GAGE program (Prinn et al., 1987). If the low CH,CCK growth rate measured 
at Hokkaido is representative of a large area of the Northern Hemisphere, then this implies that the lifetime 
of CHjCCIj is shorter than that shown in Table 2.5-1. The observed calibration difference of about 20% 
between the two data sets cannot account for this discrepancy between the observed trends. 


As stated in the OTP Report there are considerable uncertainties associated with the absolute cali- 
bration ofCH,CCI, and CCI 4 . For example, the concentrations observed on Hokkaido for these two species 
are approximately 20% lower than observations reported from the GAGE program for similar latitudes. 
Analysis of GAGE CC1„ data suggests that the Makide CC1 4 calibration (Yokata et al., 1985) is consistent 
with the calculated industrial CC1 4 emissions and an atmospheric lifetime of about 40 years (Simmonds et 
al., 1988). These calibration uncertanties remain and need to be resolved, as they introduce considerable 
uncertainty into the deduction of atmospheric lifetimes for these species. 

2.5.2.2 Other Chlorocarbons 
CCI 2 FCCIF 2 (CFC-113) 

Observations of global background mixing ratios of CC1 2 FCCIF 2 were summarized in the OTP Report 
and found to be increasing at about 11% per year in 1986. Further observations at Cape Grim (1984-87) 
show mean concentrations and increases in 1987 of 33 pptv and 3.4 ( ±0.5) pptv per year (Figure 2.5-1). 
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GROWTH OF HALOCARBONS AND NITROUS OXIDE 
AT CAPE GRIM 



Figure 2.5.1 . Halocarbon (pptv) and nitrous oxide (ppbv) observations at Cape Grim, Tasmania, as part 
of the GAGE program (Fraser and Derek, 1989). 


Based on the observed Northern Hemispheric-Southern Hemispheric concentration ratio (Ehhalt et al., 
1989), and converting the OGC (Oregon Graduate Center) scale to the NBS scale (NBS/OGC = 1.31, 
Rasmussen, 1989, private communication), an approximate global mean mixing ratio (45 pptv) and growth 
rate (10% per year) can be derived. This change in the OGC calibration of CCKFCC1F,, if substantiated, 
will essentially remove the differences between the GAGE measurements in the Northern Hemisphere and 
those made on Hokkaido, Japan (for discussion, see OTP Report). 

CHCIF 2 (HCFC-22) 

Recent observations of global background mixing ratios of CHCIF? have been reviewed (Ehhalt et 
al., 1989) and found to be increasing at about 7% per year. Further observations at Cape Grim (1984-87) 
show a mean concentration and increase in 1987 of 91 pptv and 6.5 (±0.3) pptv per year (Figure 2.5-2). 
Based on the previously observed global distribution (Rasmussen and Khalil, 1982) an approximate global 
average mixing ratio ( 100 pptv) can be derived (Table 2.5-1). 

CH 3 CI 

Methyl chloride (CH,CI) is the most abundant halocarbon in the atmosphere, and is thought to be 
largely of natural origin. Observations of global background mixing ratios ofCH,CI have been reviewed 
(Ehhalt et al.. 1989), suggesting a global background level of 600 pptv. Further observations at Cape Grim 
( 1984-1987) show a 1987 mean level of around 570 pptv and a clear annual cycle (Figure 2.5-3), which has 
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Figure 2.5-2. CHCIF 2 (pptv) observations at Cape Grim, Tasmania, from the Oregon Graduate Center 
flask sampling program (Fraser et al., 1989). 
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Figure 2.5-3. CH 3 CI (pptv) observations at Cape Grim, Tasmania, from the Oregon Graduate Center 
flask sampling program (Fraser et al., 1989). 


been attributed to seasonal variations in the magnitude of the CH,CI-OH radical sink (Ehhalt et al., 1989). 
Cape Grim record is too short to deduce whether there exists a long-term regional trend in CH,C1. 

CHCI 3 

Global background observations of CHC1, from the OGC flask-sampling network have been published 
(Khalil et al., 1983), suggesting a global average mixing ratio of approximately 30 pptv. Further observations 
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at Cape Grim (1984-87) show a large interannual variability. There is also a suggestion of an annual cycle. 
The mean level at Cape Grim was about 10 pptv (Fraser et al., 1989; Figure 2.5-4). 

2.5.2.3 Bromocarbons 

Bromine species are relatively minor components of the stratosphere compared to chlorine species 
but are nevertheless important because bromine radicals can very efficiently catalyze ozone destruction 
(McElroy et al., 1986). 

The major sources ot stratospheric bromine are the naturally occurring organobromine species such 
as methyl bromide (CH,Br). methylene dibromide (CFFBr,) and bromoform (CHBr,), produced by oceanic 
algae (Rasmussen and Khahil, 1984; Penkett et al., 1985; Khalil and Rasmussen, 1985; Khalil et al., 1987; 
Cicerone et al., 1988; Class and Ballschmiter, 1988) and the anthropogenic organobromine species, bro- 
modifluoromethane (CBrCIF,, halon-1211) and bromotrifluoromethane (CBrF,, haIon-1301), used exclu- 
sively as fire fighting agents, and ethylenedibromide (CFFBrCfCBr), a gasoline additive. 

The concentrations ot the halons have been growing rapidly in recent years (Lai et al., 1985), and 
there could be a substantial anthropogenic contribution to the methyl bromide observed in the atmosphere 
Penkett et al. (1985) measured average concentrations of 15.4 ±1.9 pptv and 10.6 ±0.9 pptv for CH, Br in 
the Northern and Southern Hemisphere, respectively, on ship voyages down the Atlantic Ocean from 
Southampton to Antarctica in 1982 and 1983. They estimated a total source strength of approximately 10 s 
tons per year to account for the data, ot which 75% could be of anthropogenic origin (Wofsy et al., 1975). 

The significance of the growing halon concentrations and a potentially large anthropogenic component 
in the measured methyl bromide is that these are the predominant forms in which bromine is transferred 
to the stratosphere, prior to breakdown into Br atoms and BrO radicals. 

2.5.3 Nitrous Oxide 

Nitrous oxide (N : 0) is an important component of the background atmosphere, being a climatically 
significant species (Ramanathan et al., 1985), as well as a source of nitrogen oxides in the stratosphere. 



Figure 2.5-4. CHCI 3 (pptv) observations at Cape Grim, Tasmania, from the Oregon Graduate Center 
flask sampling program (Fraser et al., 1989). 
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which play a major role in regulating stratospheric ozone levels (Crutzen, 1970). Natural sources of N 2 0 
include microbiological processes on the continents and in the oceans, and anthropogenic sources arise 
from the combustion of fossil fuels (Weiss and Craig, 1976), biomass burning (Crutzen et al., 1979), an 
the use of nitrogenous fertilizers (Crutzen and Ehhalt, 1977). 


There are several long-term N-0 measurement programs on the background atmosphere at various 
locations around the globe. Data from the SIO (Scripps Institute for Oceanography) network showed that 
the global rate of increase in 1986 was 0.63 (±0.07) ppbv per year (Ehhalt et al., 1989). The Northern 
Hemispheric data from this program are shown in Figure 2.5-5 (Weiss, 1989), which show an average 
increase of 0.60 ( ±0.06) ppbv per year (Table 2.5-2). 


N 2 0 measurements have been made on air samples collected at the GMCC (Geophysical Monitoring 
for Climatic Change) sites (Barrow, Niwot Ridge, Mauna Loa, Samoa, and South Pole) since 1977 ( Bodhaine 
and Rosson, 1988). The global average mixing ratio in 1987 was 306.7 ppbv, increasing at 0.73 (±0.05) 
ppbv per year (Table 2.5-2). The globally averaged data are also shown in Figure 2.5-5. 

A long-term N : 0 data record has been collected at Cape Point, South Africa, by scientists from the 
Fraunhofer Institute for Atmospheric Environmental Research (FIAER), FRG, and CSIR, South Africa, 
using an in situ gas chromatograph (Figure 2.5-5). The 1987 average mixing ratio was 301.2 ppbv, increasing 
at 0.6 (±0.2) ppbv per year (Table 2.5-2; Scheel et al., 1988). 


A combined OGC-GAGE data analysis showed a global increase of 0.8 ±0.1 ppbv per year (Ehhalt 
et al 1989). Further analysis of the GAGE N.O data has been delayed by the GAGE calibration nre- 
evaluation. However, the preliminary GAGE N.O data from Cape Grim in 1987 (Fraser and Derek, 1989) 
have been analyzed and show an average mixing ratio of 307.8 ppbv increasing at 0.95 ( ±0.06) pp v per 
year (Table 2.5-2, Figure 2.5-5). The trends deduced from GAGE data are significantly larger than those 


Table 2.5-2. 

N,0 mixina ratios and rates of increase observed in the troposphere 





Rate of Increase 




N 2 0 (ppbv) 

ppbv/yr 


Reference 

Site 

Laboratory 

1987 

(±2<r) 

%/yr 

Global 

Average 

GMCC 

306.7 

0.73 

(±0.05) 

0.24 

Bodhaine and 
Rosson (1988) 

NH Average 

SIO 

305.1“ 

0.60 

(±0.06) 


Weiss (1989) 

Global 

Average 

SIO 

305.1“ 

0.63 

(±0.07) 


Ehhalt et al. 
(1989) 

Cape Point 
S. Africa 

FIAER- 

CSIR 

301.2 

0.62 

(±0.20) 

0.21 

Scheel et al. 
(1989) 

Global 

Average 

OGC 

308. 2 h 

0.83 

(±0.05) 

0.27 

Ehhalt et al. 
(1989) 

Cape Grim 
Tasmania 

CSIRO 

307.8 

0.95 

(±0.06) 

0.31 

Fraser and 
Derek (1989) 


a 1 986 data 
b 1985 data 
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R Hnl m . L ° ng '?I m trends of N 2° from the GMCC program (average of data from Barrow, Niwot 
*' dge \ Ma „ Una L ° a ' and Sf moa ). SIO program, the OGC-ALE program (average of data from Barrow 

pS in S Meare ®' Loa Barbados, Samoa and Cape Grim), the FIAER-CSIR program (Cape 

Tqoq*’ 1777 *'7 and CSIR0 ' GAGE data ( Ga P e Grim - Tasmania) (Bodhaine and Rosson, 1 988; Weiss 
1989, Ehhalt et al., 1989; Scheel et al., 1988; Fraser and Derek, 1989). 


from SIO, GMCC or FIAER-CSIR data (Table 2.5-2), but the significance of these differences cannot be 
fully assessed I until the GAGE re-calibration .s completed. The GAGE calibration scale is also significantly 
higher than those of GMCC and SIO, but this cannot account for the trend differences observed. The 
globally averaged trends (Table 2.5-1) are based on SIO and GMCC data. 


. k ^ he ° TP ReP ° rt da ‘ a WCre Sh ° Wn ° n ICC COre Studies of N ^° which ’ des P ite a nois y record, showed 
that N.O concentrations prior to 1800 were about 280-290 ppbv. Further data have been published (Khalil 

an Rasmussen, 1988a), that confirm this result with more precise data which show that N.O levels have 
risen by about 8% during the last century. 


2.5.4 Methane 


Methane is an important trace gas because it absorbs infrared radiation, and thus impacts on climate 
It affects tropospheric chemistry via hydroxyl radical and ozone budgets, and it impacts on stratospheric 
chemistry as a source of hydrogen and water vapor and as a sink for stratospheric chlorine. 

Recent data from three laboratories making global measurements (GMCC, CSIRO and UC (Univer 
sity of California)) are summarized in Table 2.5-3. The GMCC and CSIRO programs ^how significantly 

'77 7x7 Tr^lnT * PPbV P 7 yCar) tHan thC UC P r °8 ram <16 ppbv per yr). From the GMCC (Steele 
et al. 1989) and CSIRO (Fraser et al., 1989) programs it would appear that the global rate of CH 4 increase 

has slowed ,n recent years. Figure 2.5-6 shows the data from the CSIRO program as well as data collected 
at Cape Point, South Africa (Scheel et al., 1988), which also show a current low rate of increase (10 ppbv 
per year). The cause of the recent slowing of the CH 4 increase is not known. Clearly longer data records 
are necessary to properly define long-term CH 4 variability. 
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Table 2.5-3. Global average methane mixing ratios and rates of increase observed in the troposphere 


by various laboratories 


Laboratory 

Mixing 

Ratio 

ppbv 

Rate of Increase 
ppbv/yr 

(±2 a) %/yr 

Reference 

GMCC 

1674 a 

12.1“ 

0.7 

Steele et al., 
1987, 1989 

UC 

1680 b 

16+1 

1.0 

Blake and Rowland, 
1988 

CSIRO 

1 682 c 

12 ± 2 

0.7 

Fraser et al., 
1986b, 1989 

FIAER-CSIR 

1612 d 

10.3 + 0.4 

0.6 

Scheel et al., 
1988 


a 1 988 data, OGC calibration 
b 1 987 data, UC calibration 
c 1 987 data, NBS calibration 
d 1987 data, FIAER calibration 


The zonally averaged distribution of methane in the remote marine boundary layer is shown in Figure 
2.5-7, from 5 years of data from the GMCC sampling network, comprising 26 globally distributed sites 
(Steele et al., 1987). The seasonality in the Southern Hemisphere is very repeatable, as is the average 

latitudinal gradient. 

New data have become available on past CH 4 levels. Using infrared solar spectra recorded at the 
Jungfraujoch Station, Switzerland, in 1951-52 and 1984-87 it has been shown that the total vertical column 
abundance of CH 4 increased by 0.7(±0.2)% per year during this period (Zander et al., 1989a). Ice core 
studies of historic CH 4 levels were reviewed in the OTP Report, which showed that over the past 300 years 
CH 4 levels have more than doubled. Further studies have been published which extend the CH 4 record 
back more than 100,000 years. Stauffer et al. (1988) show that the concentration of CH 4 was about 500 
ppbv 100,000 years ago and fell to around 350 ppbv near the peak of the last ice age, 18.000 years ago, 
before rising to a typical Holocene value of about 650 ppbv. Raynaud et al. (1988) have extended the CH 4 
record back to 160,000 years before present, and have shown that CH 4 was similarly low at the end of the 
preceding ice age. These results suggest a relationship between CH 4 levels and climate that may provide 
further clues in understanding modern CH 4 variability. 


2.5.5 Carbon Monoxide 

Carbon monoxide is an important trace gas in the troposphere because it plays significant roles in 
controlling the chemistry of ozone production and hydroxyl radical destruction in the lower atmosphere. 
It directly affects the oxidizing capacity of the lower atmosphere and can thus influence the concentrations 
of climatically important gases such as methane, for which oxidation by hydroxyl radical is a major sink. 

Because of its short atmospheric residence time (2-3 months) the determination of a secular trend in 
CO mixing ratio is difficult. Increases are expected because significant CO sources are under human 
control, such as fossil fuel use, and have been growing (Cicerone, 1988). The available global CO data have 
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Figure2.5-6. Long-termtrendsof CH 4 fromtheCSIROglobalprogram(Fraseretal., 1986a 1989)and 

the FIAER-CSIR program at Cape Point, South Africa (Scheel et al., 1988). 


been reviewed (Cicerone, 1988; Ehhall et al., 1989) suggesting that CO has increased in the Northern 
Hemisphere over the past 15-30 years, but data from the Southern Hemisphere show no significant or 
uniform trends. The data have been updated in Table 2.5-4 to include observations from the OGC and 
CSIRO sampling networks, as well as from Cape Point, South Africa. The CSIRO and FIAER-CSIR data 
are shown in Fig. 2.5-8. 
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Figure 2.5-7. The global distribution, seasonality, and trend of CH 4 from the GMCC network (Steele et 
al., 1987 and unpublished data). 


All data sets show statistically significant positive CO trends in the Northern Hemisphere of about 
1% per year, except the relatively short records collected by CSIRO at Barrow and Mauna Loa, which do 
not show significant trends. The Southern Hemispheric data do not show consistent trends. The data 
collected by OGC (Khalil and Rasmussen, 1988b), over the period 1980-87, and CSIRO (1984-88) show a 
significant, positive trend at Samoa, but not at Cape Grim (CSIRO: 1978— 88). The relatively long records 
at Cape Point (1978-88) (Scheel et al., 1988) and Mawson (1980-88) do not show significant trends, whereas 
the relatively short CSIRO record at the South Pole shows a significant, positive trend. The mid-to-high 
latitude data from the Southern Hemisphere show remarkably uniform (amplitude and phase) seasonal 
variations (late summer minimum, winter maximum), which are probably in response to the hydroxyl 
radical sink (Fraser et al., 1986a, b; Ehhalt et al., 1989). 

The secular increase of the total vertical column abundance of CO over Europe has been derived 
from sets of infrared solar spectra recorded at Jungfraujoch Station, Switzerland, in 1950-51 and in 1985— 
87. The mean cumulative rate of increase of the total column abundance of CO over this 35-year period 
was found to be 0.85 (±0.2)% per year (Zander et al., 1989b). As Cicerone (1988) has pointed out, the 
detection of sustained temporal trends in atmospheric CO is very challenging. Atmospheric CO concentra- 
tions are strongly influenced by sources and sinks that vary greatly with season, latitude, and longitude. A 
successful trend-detection strategy must account for this variability and be capable of dealing with differ- 
ences between continental and marine regions and between the planetary boundary layer and the free 
troposphere. 
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Table 2.5-4. Surface carbon monoxide concentrations and trends 





Cone. 

Trend 





ppbv 

%/yr 


Site 

Lat. 

Period 

1987 

(±2 a) 

Reference 

Moscow 

56°N 

1952,53- 


1 

Dianov-Klokov & 



1970,76 



Yurganov, 1981 

Northern Hemisphere 


1970-82 


1.5 

Dvoryashina et 
al., 1984 

Jungfraujoch 

46°N 

1950,51- 


0.85( ±0.2) 

Zander et ah, 



1985-87 



1989b 

New Jersey 

40°N 

1968-77 


positive 

Graedel & 
McRae, 1980 

OGC network 1 


1980-87 



Khalil & 

Northern Hemisphere 



105 

0.9( ±0.5) 

Rasmussen, 1988b 

Southern Hemisphere 



50 

1 .9( ± 0.8) 


Cape Point h 

34°S 

1978-88 

60 

0. 1( ±0.2) 

Scheel et al., 
1988 

CSIRO network 
Barrow 

7I°N 

1984-88 

155 

- 1.5(±2.3) 

Fraser et al.. 

Mauna Loa 

20°N 

1984-88 

90 

0.2(±2.9) 

1986a, b 

Samoa 

14°S 

1984-88 

65 

3.4(±2.0) 

Fraser, 1989 

Cape Grim 

41°S 

1978-88 

60 

— 0. 1 ( ± 0.7) 


Mawson 

68°S 

1980-88 

60 

0.2( ±0.9) 


South Pole 

90°S 

1984-88 

55 

2.5(± 1.8) 


a OGC scale 
'’FIAER scale 
C NBS scale 


2.5.6 Carbon Dioxide 


Carbon dioxide (CO : ) is the major greenhouse gas in the atmosphere whose abundance is changing 
in response to human activities. CO ; plays a significant role in regulating the temperature of the stratosphere, 
which in turn influences stratospheric ozone chemistry and abundance. Precise atmospheric measurements 
ofC0 2 have been made since the late 1950s, initially in Antarctica and at Mauna Loa, Hawaii, and presently 
by an international network of more than 40 sites. The increase of CO. in the background atmosphere over 
the past 15 years has been well documented. 

Carbon dioxide data for 1987 from in situ analyzers are available from the GMCC network (Bodhaine 
and Rosson, 1988) and from Cape Grim, Tasmania (Beardsmore and Pearman, 1989). The global average 
mixing ratio in 1987 deduced from these sites is 348.1 ppmv, with a global average increase of 1.6 ppmv 
from 1986 to 1987. The long-term trends of CO : (1973-1987) from the GMCC sites are shown in Figure 2.5- 
9. On average, the global increase over this period has remained steady at about 1.5 ppmv per year, but at 
certain times and sites increases as large 3 ppmv per year and as small 0 ppmv per year have been observed. 
There are strong correlations between the rates of increase observed at these widely distributed sites. 
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Figure 2.5-8. Long-term trends of CO from the CSIRO network and from Cape Point, South Africa (Fraser 
et at., 1986a,b, 1988; Fraser and Coram, 1988; Scheel et al., 1988). 


Measurements of the global distribution of C0 2 from a network of 26 flask sampling sites (Conway et 

a,.. 1988) continued in .987 (Bodhaine and Rosson, 1988). The ^ 

this network was 348.3 ppmv. and the global average CO, increase (1986 to 1987) was 1.9 ppmv. 


2.5.7 Tropospheric Ozone 

Tropospheric ozone impacts in several important ways on the tropospheric levels of source gases 
(Levy 1971). Most notably, it is essential for the photochemical production of the hydroxyl radical, OH, 
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YEAR 

Ros^on 2 '^^). 1 ^ behavior of the c °2 trends observed at the GMCC stations (Bodhaine and 


Which determines the lifetime of many tropospheric trace gases, and thus also in part determines the ozone 
depletion potential of the hydrogen-containing source gases. In addition, tropospheric ozone makes a 
contribution ot about 10% to the total ozone column in the atmosphere. Thus, a trend in tropospheric- 
ozone has important consequences for stratospheric ozone concentrations and the detection of stratospheric 
ozone trends. I ropospheric ozone is also a significant greenhouse gas (Ramanathan et al., 1985). 

In fact, a trend in tropospheric ozone, particularly in the Northern Hemisphere, is not unexpected 
Ozone isproduced photochem.cally in the troposphere and its precursors, nitrogen oxides, CO CH 4 and 

'° ha,C increascd ' due ,0 a " ,hropoge " k emi5sions 


I9m-^k«| de ™. f ‘?j ncrc “^, in ! rOPOSpheriC OZOne has becn rcv ' ewcd "-ogan, 1985; Bojkov. 1986. 
, Penkett. 988). The available observations on current trends in ozone at or near the Earth’s surface 

— anZ A.; n p 5 ' 5 They Seem *° SUPPOrt the expected P° sitive in regions of the Northern 

Hemisphere. All seven European sites show statistically significant positive trends. However all four 

Canadian sites show negative trends, ot which only one is significant. All three Japanese sites show positive 

I h h S r hern He « iSpheriC SiteS d ° n0t ShoW statistica, 'y significant trends, 
except for the South Pole, which shows a significant negative trend. The central European trends reach as 

high as 3 \A> per year. It should be noted that these trends fluctuate in time and that some of the European 

stations have exhibited a slowing down of the positive trends and even a decrease of ozone in recent years. 

Only a few of the stations listed in Table 2.5-5 are clean air sites, making semi-continuous surface 

— '“ Sl dai ! yb " amdy , B ™ M— Loa in the Northern Hemisphere and Samoa. 

<pe Point, Cape Grim (Galbally et al., 1986) and the South Pole in the Southern Hemisphere These 
stations are sited to represent large regional air masses. Data from these stations are shown in Figure 2.5- 
o stations in the Northern Hemisphere show significant positive trends of 0.8% per year, indicating 
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Table 2.5-5. 


Surface ozone concentrations and trends deduced from ground based and balloon 


Site 


Latitude 

Period 

Average 

Mixing 

Ratio 

(PPbv) 

Trend 

%/yr 

(±2<r) _ 

75°N 

1966-85 


— 0. 1( ±0. 1) 

7PN 

1973-87 

26 bc 

0.8( ±0.4) 

59°N 

1974-85 


-0.7(± 1.3) 

55°N 

1956-84 

I9 d 

1 ,4( ±0.8) 

54°N 

1973-85 


— 0. 1 ( ± 1 .4) 

53°N 



-l.l(±1.0) 

51°N 

1952-84 

I0 d 

2.6( ± 1 .6) 

51°N 

1955-83 

I6 d 

3.l(±2.l) 

5I°N 

1972-83 

21 d 

3. 1( ±2.4) 

50°N 

1954-84 

23 d 

1 . 1( ±2.0) 

48°N 

1974-86 

31 hd 

1 . 1( ±0.3) 

48°N 

1969-86 


2.3( ±0.5) 

47°N 

1968-86 


2.4( ±0.8) 

43°N 

1969-86 


0.9( ±1.1) 

36°N 

1969-86 


1 ,3( ±0.6) 

32°N 

1969-86 


2.5(± 1.0) 

20°N 

1973-87 

IT- 

0.8( ±0.4) 

14°S 

1976-87 

14 hc 

— 0. 3( ± 0.7) 

34°S 

1982-88 

20 c 

0.3( ± 2.0) 

4I°S 

1982-86 

24 f 

0.6( ±0.7) 

90°S 

1975-87 

20 hc 

-0.5( ±0.4) 


Northern Hemisphere 

Resolute a 

Barrow 

Churchill 1 * 

Arkona 

Edmonton 1 * 

Goose 1 * 

Dresden 

Kaltennordheim 

Gr.Inselberg 

Fichtelberg 

Hohenpeissenberg 

Hohenpeissenberg 11 

Payerne 1 * 

Sapporo 1 ' 

Tateno 1 * 

Kagoshima 11 
Mauna Loa 

Southern Hemisphere 

Samoa 
Cape Point 
Cape Grim 
South Pole 


a sonde data 
b nbar 

c Oltmans and Kornhyr, 1986; Bodhaine and Rosson, 1988 
rt Feisterand Warmbt, 1987 
e Scheel et at. , 1988 

f Galbaily et al., 1986; Elsworth et aL, 1988 _ 


that tropospheric ozone increases can be observed in sites remote from industrial emissions (Oltmans and 
Komhyr, 1986; Bodhaine and Rosson, 1988). 

There are nine stations, two in Europe, three in Japan, and four in Canada which provide measure- 
ments of ozone in the free troposphere using balloon-borne ozonesondes. The data have been reviewed in 
Section 2.3-3 and are listed in Table 2.3-4. As indicated above, all European and Japanese stations exhibit 
positive trends and all Canadian stations exhibit negative trends. The Canadian data might be affected by 
an instrument change that occurred at all stations in the early 1980s (Bojkov, 1988). A previous analysis of 
these data (Tiao et al., 1986) showed some evidence for a positive trend of 0.8% per year in the lower 
troposphere. Averaging over all nine stations the current analysis shows a significant trend in the first three 
kilometers of 0.8( ±0.5)% per year. Average trends in the troposphere above this altitude were not sigm - 

icant. 
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Figure 2.5-1 0 . Surface ozone measurements at the GMCC stations (Olmans and Komhyr, 1 986- Bodhaine 
and Rosson, 1 988), at Cape Point, South Africa (Scheel et al. , 1 988), and at Cape Grim, Tasmania’fGalbally 
et al., 1986; Elsworth et al., 1988). 1 


In addition to the recent tropospheric ozone measurements, there are a few historic records of surface 
ozone measurements. The most important is a complete series of daily means collected at the Montsouris 
Observatory in Paris between the years 1876 and 1905. These data were re-analyzed by Bojkov (1986) and 
Volz and Kley (1988). The annual averages were shown to vary between 5 and 15 ppbv, about a factor of 
two lower than present ozone levels in a comparable environment. 

None of the various pieces of evidence obtained so far and reviewed above provide convincing proof 
that ozone is increasing on more than a regional scale. When combined, however, they do suggest that 
ozone has been increasing over a wide scale in the Northern Hemisphere. On the other hand, the few 
measurements existing in the Southern Hemisphere indicate more or less constant ozone levels, which is 
not unexpected due to the low levels of ozone precursors observed in the Southern Hemisphere. 

It has been suggested that typical tropospheric ozone variability means that the current ozonesonde 
network, and presumably the network of surface ozone measurement sites, are insufficient to detect a 
global trend in tropospheric ozone of about 1% per year at the 2 o confidence level, even at stations with 
records a decade in length (Prinn, 1988b). To detect an ozone trend of 1% per year on a decadal time frame, 
daily measurements of ozone are required from approximately 40 sites, chosen to represent continental 
and marine, tropical, temperate, and polar air masses in both hemispheres (Prinn, 1988b). This will require 
a doubling of the current network. Measurements are particularily lacking from Africa, South America, 
and Asia. It would appear that similar requirements may be necessary for a suitable global CO network. 
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2.6 TRENDS IN STRATOSPHERIC AEROSOLS 
2.6.1 Introduction 

Chapter 10 of the Ozone Trends Panel Report (OTP) presented a detailed description of stratospheric 
aerosols, their properties, and their possible effects on ozone and ozone observations. The specific objec- 
tives of the chapter were: ( 1 ) to define the fundamental characteristics of stratospheric aerosols, particularly 
their morphological, radiative, and chemical properties; (2) to analyze long-term aerosol data bases for 
possible trends that might influence the identification or interpretation of long-term ozone trends; (3) to 
assess the transient impacts of volcanically generated aerosols on remote sensing observations such as 
Umkehr and SBUV measurements over the last two decades, and to consider possible schemes for 
correcting aerosol errors in these baseline ozone data records; and, (4) to describe polar stratospheric 
clouds (PSCs) in terms of their frequency, spatial extent duration, and physicochemical properties, on the 
basis of observational and theoretical studies. 


The purpose of this section of the present assessment is to review the findings of the OTP Report, 
and bring the reader up to date regarding any new or extended aerosol data. Chapter 1 of this assessment 
on polar ozone described in detail the climatology of PSCs in both polar regions, their physical properties, 
and heterogeneous phase chemistry involving PSCs. Thus, this section will not include any discussion ot 
PSCs. 

2.6.2 Stratospheric Aerosols 

With the data available at the time of the OTP Report, no significant trends were detected in the 
properties of the global background stratospheric aerosol layer. Figure 1 .2-30 in Chapter I of this document 
plots the SAM II stratospheric optical depth at a wavelength of I p.m for the period November 1978 through 
February 1988 for the Arctic and Antarctic regions. Figure 10.7 in the OTP Report showed the SAM II 
data record through only 1986. Evident from these figures is the rapid impact of volcanic eruptions and 
their subsequent slow decay. Over this time period, the strongest impacts were experienced in the Arctic 
but for the case of El Chichon (17.3°, 93.2°W), which erupted in April 1982, both hemispheres were strongly 
perturbed. 

These enhanced aerosol levels from volcanic eruptions can dramatically and deleteriously impact 
remote sensors. For example, after the eruption ot El Chichon, sea surface tempeiatuie retrievals from 
AVHRR were adversely impacted by as much as a few degrees (Strong, 1984). Similarly, the retrieval of 
ozone profiles from Umkehr and SBUV measurements were affected, complicating the interpretation of 
these remote sensing data, especially when used to determine long-term changes or trends. Obviously, as 
Figure 1.2-30 shows, the El Chichon eruption caused the largest perturbation to stratospheric aerosols 
experienced over this last decade, but even the impact of smaller eruptions must be understood in order 
to remove their effects from various data sets. The problem is to determine when the aerosol loading is 
low enough not to adversely affect the interpretation of the particular sensor’s data. The last eruption 
thought to have caused a level of perturbation similar to that experienced after the eruption of El Chichon 
was the March 1963 eruption of Agung (8.3°S, 1 15.5°E). 

Routine balloon and lidar aerosol measurements of the stratosphere began in the late 1960s and early 
1970s. These long-term measurements have been obtained primarily in the Northern Hemisphere with 
some short-term data records in the Southern Hemisphere. All these measurements have shown the period 
of 1978 to 1979 to be lowest on recent record. Therefore, this period has been referred to by most 
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investigators as a background period for stratospheric aerosols. In the strictest sense, however, with a data 
record of only 20 years, no one can definitively claim this to be a true background period. 

As of this writing, mid- 1989, the SAM II record and the lidar record from the NASA Langley Research 
Center shown in Figure 2.6. 1 , as well as other similar records, show values that have gradually approached, 
albeit not yet reached, the values measured in 1978-1979. Since the El Chichon eruption, the only other 
eruption to significantly perturb the global stratospheric aerosol was the eruption of Ruiz (4.9°N, 75.4°W) 
in November 1985. The Ruiz enhancement was similar to the enhancement caused by the Mount St. Helens 
(46.2 N, 122.2 W) eruption of May 1982 and, therefore, most likely has delayed somewhat the recovery 
from the El Chichon eruption. For example, the present value at Langley for peak backscatter ratio at a 
laser wavelength ot 694 nm, is about 1.2, which represents a 20% enhancement over molecular scattering. 
In early 1979, this ratio was 1 .08 for an 8% enhancement. These are to be compared to the peak backscatter 
value in late 1982 of 24, which was caused primarily by El Chichon. Integrated backscatter over the 
stratosphere is also used to describe lidar results. It provides a measure of stratospheric total column 
loading versus time over the lidar site. Values in 1989 have fallen to 0.8 10 4 sr 1 from peak values in early 
1983 of 27 10 4 sr '. Values in early 1979 were only 0.4 10 4 sr '. The SAM II results are similar, with 
lowest Antarctic optical depths measured in October 1979 of 0.5 10 4 compared to 1.7 10 4 in October 
1988, and lowest Arctic optical depths of 1.2 10 1 in February 1979 compared to 1.7 10 4 in February 1989. 
Although the 1978-1979 aerosol minimum values have not yet been reached, there is no hard evidence that 
they will not be reached again, and therefore, there is no evidence that an upward trend in “background" 
stratospheric aerosols has occurred. 

2.6.3 Global Heterogeneous Effects 

Recently, Hofmann and Solomon (1989) suggested that stratospheric aerosols, especially during 
periods enhanced by volcanic eruptions, could possibly cause an ozone depletion through heterogeneous 
chemistry similar to that occurring on PSCs. They reached this conclusion because recent laboratory 
studies suggest that the same reactions that occur on PSCs in the cold Antarctic stratosphere, can also 



Figure 2.6-1 . Aerosol integrated backscatter coefficient over the period 1 974-1 988 as observed by the NASA 
Langley Research Center lidar. 
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i SH2££5S=£53 

the El Chichon eruption reactions will be of greatest potential importance in the middle-to-high 

calculations indicate that the Q^rtion ~> are largest and the ozonosonde 

latitudes in winter, where the reported total ozone trends sSia, deereases in NO, 

measurements show record lows after El • ........ - winter Recall from Chapter 1 that 

abundances following El Chichon. parhcularly at nuddle huttude he , e P ous ozone 

mll9 LIMS HNO, measurements but not wtlh current -non-hetcrogeneous models. 

The possibility of this 

dally after a large volcanic eruption, makes it even more urgent 
geneous reactions as they apply to ozone depletion. 


2.7 SURFACE ULTRAVIOLET RADIATION 


Predicted Trends Related To Column Ozone, 1970-1986 


Ehbui£^VfTnro P f^ r s ?rb e S °t S sca"lte:sTs n weH r l d j “\hnur^ce feldTo fttond .n 

ultraviolet irradiance at the Earth s surface. 

,,r:=:x^q 

300 nm. Backscattering of sunlight to space by clouds reduces linearly with 

for clear skies (Mo and Green, 1974; Ilyas, 1987). 

latitudes correspond to the centers otea " ^ US y , , noon and clear sky conditions. The 

52, and 53-64 degrees North. All values in Figure ^.7 I r characte rized by large variations with 

natural ultraviolet radiation environment at t e ar s s . horter t h an 330 nm results from 

month and latitude. The sharp decrease m >™ d,anc f * Huggins bands, 
absorption by ozone, and spectral structure between 310 and 330 nm arises trom 
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2s^ o’o'S 1 **£!’ r':^" cei '"T T ! hei Ear,hS: surtacecom P^ d 101 j ™ry and July alia, itudes 
• , 4b 0 , and 58.5 N. Calculations refer to clear skies, local noon, and serve as baseline cases 

appropriate to the year 1970. The top three curves refer to July and the bottom three ret to January. 


The to™ ^ rends u Panel re P° rted a percent change in column ozone for each month of the year over 
Period 1970 through 1986 in the latitude zones defined above. Table 2.7-1 lists these for the three 
atuudes represented in Figure 2.7-1 during the months of January. March, and July. Imposed on the long 

7 arC vanat, ° ns rela,ed 'o the solar cycle, the QBO, and a random interannual variability The 
calculations reported below include only the linear trend component. 

F'gure 2.7-2 presents the percent changes in noontime surface irradiance attributable to the trend in 
otal ozone between ,970 and ,986. These calculations refer to clear skies, include only the decease In 
column ozone, and neglect changes in the profile shape. All curves in Figure 2.7-2 show the expected 
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Table 2.7-1. 


Month 

January 

March 

July 


Percent Changes in Total Column Ozone, 1970—1986, Used in the Radiative Transfer 
Calculations* 


30-39°N 

(%) 

Latitude Band 
40-52°N 
(%) 

53-64°N 

(%) 

-2.2 

-2.6 

-8.3 

-3.5 

-5.6 

-4.0 

- 1.3 

-2.2 

0.0 


'Values from NASA (1988) 


increase in percent change as wavelength decreases. In addition, a given percent change in ozone has an 
increasing impact on surface irradiance as the sun nears the horizon. In general, the largest percent changes 
occur at wavelengths where the absolute irradiances are small. At wavelengths longer than 330 nm, a trend 
in total ozone has negligible influence on the surface irradiance. 

When the sun is high in the sky, as in summer, it is possible to experience a decrease in ultraviolet 
irradiance at the ground even with a decline in column ozone. This would arise from a change in the shape 
of the ozone profile, not included in Figure 2.7-2, and requires an increase in tropospheric ozone, which 
only partially offsets a reduction in the column at higher altitudes (Bruhl and Crutzen, 1989). The presence 
of aerosols and clouds leads to efficient scattering of the direct solar beam. This increases the effective 
path length taken by radiation through the lower atmosphere and consequently enhances absorption by 
tropospheric ozone (Frederick and Lubin, 1988). 


2.7.2 Observed Trends in Surface Ultraviolet Radiation 

Scotto et al. (1988) performed a trend analysis of ultraviolet irradiance data collected by a network 
of eight Robertson-Berger (RB) meters located over the United States from latitudes 30.4 to 46.8 degrees 
North. The time period covered was 1974 through 1985. The RB meter responds to broadband radiation 
convolved with an instrument response function which covers the spectral region from the ozone cutoff 
near 295 nm to approximately 340 nm. The RB data base indicates no increase in irradiance over the \~- 
year period, and decreases statistically significant to the ±2 ct confidence level or higher occurred at five 
of the stations. The discussion below assumes that the derived trend in irradiance is a real atmospheric 
phenomenon. Flowever, a thorough evaluation of the long-term stability of the RB meter network would 
be of value. 

The numerical analysis of Scotto et al. (1988) was based on the annually integrated RB meter response 
at each station. The large annual cycle evident in Figure 2.7-1 implies that the annually integrated value 
receives its greatest contribution from the summer months. This is the season in which Dobson data show 
little change in total ozone over North America between 1970 and 1988 (see Table 2.2-5). However, Scotto 
et al. (1988) also examined RB readings by individual months and still found decreases in surface irradiance. 
The output of the RB meter includes substantial contributions from solar irradiance at wavelengths where 
absorption by ozone is weak. The meter is therefore not a sensitive indicator of small changes in atmospheric 
ozone content. With the RB meter response function from Caldwell et al. (1986), the 2.2% decrease in 
ozone for 40-52 degrees North during July between 1970 and 1986 implies a 1 .5-1 .6% increase in RB meter 
response. The 5.6% decrease in ozone during March would appear as a 4.7-4. 8% increase in RB output. 
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Figure 2.7-2. Percent changes in solar spectral irradiance between 1 970 and 1 986 computed for January, 
March, and July at latitudes 34.5°, 46.0°, and 58.5°N. Calculations include only the linear trend in ozone 
for a given month between 1970 and 1986. The percentage changes in ozone appear in parentheses. 
Curves for 46.0°N and 58.5°N in March are indistinguishable. 


Given that cloud cover, aerosols, localized sources of air pollution, and likely increases in tropospheric 
ozone influence ground-based radiation measurements, it is not surprising that the RB meter network did 
not observe an increase in irradiance related to the downward trend in total ozone. The RB results are 
nonetheless very significant. In particular, they show that a trend in total ozone need not appear as an 
obvious trend in surface ultraviolet irradiance. The role of cloudiness and other factors that influence the 
transmission of the troposphere must be considered as well. The major conclusion here is that the results 
of the RB meters do not contradict the observed downward trend in ozone presented in Section 2.2. 
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2.8. OUTSTANDING ISSUES 

Although a large amount of data has been collected over many years by instruments whose charac- 
teristics are reasonably well known, few of these instruments were made for the purpose of determining 
trends or long-term changes. Taking into account the deficiencies of the present observation system and 
looking to the future, priority should be given to a system approach based on the complementarity of 
ground-based and satellite observations. Both total column and vertical distribution satellite measurements 
of ozone should be checked for instrumental drifts by comparison with a network of carefully calibrated 
ground-based instruments based on both passive— Dobson. UV, IR, and microwave spectrophotometers— 
and active-lidar systems. The implementation of such a Network for Detection of Stratospheric Changes 
requires an internationally coordinated effort with a long-term commitment for scientifically driven oper- 
ation of the stations. Special attention should be paid to the present dramatic lack of measurements outside 
of the band from 30 to 60 degrees in the Northern Hemisphere as reflected in particular in the absence of 
trend determination in the tropical regions. It is therefore of prime importance to establish new sites with 
an appropriate geographical distribution that takes into account the latitudinal and regional patterns already 
revealed by the present data analyses. Furthermore, the measurements should not be limited to strato- 
spheric ozone. Data on the vertical distribution of ozone in the troposphere are, for example, required to 
take into account a possible trend in this atmospheric region. Additional atmospheric variables such as 
temperature or aerosols have been shown in this report to be of direct influence on the ozone trend. Their 
measurements will also be complemented by those of selected stratospheric trace gases that are required 
to quantitatively relate the observed ozone trend to anthropogenic influence. This can be illustrated by the 
present difficulty in interpreting for example the observed variations of hydrogen chloride (HCI), over the 
past 15 years. 

Based on observed annual changes in CFCs and chemical partitioning arguments, it is expected that 
HCI should be increasing in the stratosphere at an annual rate of about 3.5%. Insufficient data exist, 
however, to confirm this important prediction. The longest data records have been collected from aircraft 
over the period 1978 to 1982 by Mankin and Coffey (1983) and from the Jungfraujoch station over the 
period 1976 to 1985 by Zander et al. (1987a). Both sets of observations were taken in the solar absorption 
mode using a high-resolution Fourier transform spectrometer and grating spectrometer, respectively. The 
Mankin and Coffey (1983) results give an indicated annual increase of 5% in column amount, but the 
estimated uncertainty in using the data to determine long-term change is 9% (2 a) per year when factors 
such as frequency of observations, length of data record, and long-term stability of absolute and relative 
calibration procedures are considered (Prinn, 1988a). The uncertainty estimate is probably conservative. 
Zander (1989, private communication) reports that the largest increase that can be derived from the ground- 
based results is 1.5% per year after seasonal effects are removed from his data spanning the 1976 to 1985 
time interval. Prinn (1988a) estimates the uncertainty in using these data for detection of long-term changes 
to ±2% per year (2 it). Similar Jungraujoch data covering the period 1985 to 1989 indicate a significantly 
larger trend than the one mentioned above (Zander, 1989, private communication). Tropospheric HCI 
could make some contribution to the measured total column but, since the Jungfraujoch station is at a 3.5 
km altitude, this effect should be relatively small. Both the aircraft and ground-based data sets show an 
HCI increase with time over the period of their respective data records which is greater than the individual 
measurements error bars, but they each imply different rates of increase, and their use in quantifying yearly 
changes leads to uncertainties that preclude verification of the expected annual HCI increase. A much 
broader space and time sample, combined with the investigations of sporadic potential natural sources of 
HCI (e.g., volcanic eruptions), is required to address this question. A less definitive, but useful, test of 
anthropogenic contributions to the stratospheric halogen budget is provided by contemporary measure- 
ments of the HCI/HF ratio. Since substantially all of the HF in the present-day stratosphere is of anthro- 
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pogenic origin, while part of the HC1 input is natural, the HCI/HF ratio at any time gives an indication of 
the anthropogenic input. Measurements covering the period 1977 to 1985 indicate mid-latitude values for 
this ratio of about 4 to 6 (Zander et al., 1987b), in satisfactory agreement with models and indirectly, 
therefore, supporting the predicted 3.5% annual change in HCI. 

Finally, if technology advances and improvements to old techniques will lead to the implementation 
of newer and better observing systems, these often require 10 to 15 years before homogeneous data series 
of at least a decade can be available for trend analyses. It is therefore critically essential that existing 
measurement programs be maintained and their homogeneity in time be improved, in addition, efforts must 
be expended to develop a full understanding of the differences between data series from the old and new 
instruments so that a full-length homogeneous data record can be established to facilitate trend analyses 
as the newer systems take over from the old ones to provide the mainline data record. 


APPENDIX 2.A 

2.A.1 Seasonal Trends Model (Total Ozone) 

The full statistical model for a monthly ozone series, y„ , t = 1,2,..., T, for the multiple trends case 
is given by: 

12 

~ 2 Mi A,/ + + y/Zi.t + yrZi2 : t-k + y + y/Z,4 : , + N, (A I) 

i I 

where: 

M-i = ozone mean in month i, i = 1 ... 12, 

I t;t = indicator series for the i ,h month of the year, 

= I if month t corresponds to month i of the year 
0 otherwise 

(3, = trend in Dobson units per year after 12/69 in month i of the year, 

R, = linear ramp function beginning in 12/69, 

(t - t„)/l2 if t > t„. where t„ corresponds to 12/69 

0 if t < t„ 

Z, , = solar series minus 120 flux units (minus 120 so that percentage trends will be based on ozone 

at mid solar cycle), 

Z: ; i k = QBO series lagged k months 

Z_v, = nuclear weapons series 

z 4:1 = intervention for calibration shift (some stations, when using the original published data), 

1 after the shift 

0 before the shift 

7i = coefficients to be derived for each Z term , i = I, 2, 3, 4 
N, = residual noise series 

There is a month-to-month correlation in ozone values, even after the seasonal and other effects in 
equation (Al) have been removed, i.e., in the N, series. This is modeled by allowing N, to be an autore- 
gressive series of order I : 
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N, = d> N,-, + e, 


(A2) 


where e, is an uncorrelated series (See Box and Jenkins, 1976; Tsay. 1984). 

The noise series N, is much more variable in winter than in summer. This is accounted for in the 
statistical model by fitting the above equation with individual months weighted inversely P^omoMl to 
the monthly variances of the N,; for example, all Januarys get the same weight (the UW&C model uses 
weights applied to the e, series rather than the N, series). 

Estimated changes in ozone reported in Table 2.2-1 are given as cumulative percent changes over the 
17-year period 1970-1986 inclusive, to be comparable to results reported m the OTP summary: 


Wi 


17-year change in month i in percent 


(A3) 


_ 17 x Pi 


x 100% 


Percent changes for the winter months (December through March), the summer months (May through 
August), or the whole year are the average of the percent changes in the corresponding months. 

In other tables, trends are given in percent per decade, with 17 years replaced by 10 years in equation 

(A3). 


2.A.2 Uniform Trends Model (Total Ozone) 

The uniform trends model is similar to equation (1), except that the coefficient 0 of the ramp function 
is the same for all months of the year: 


12 

V/ = 2 M'f'T + fiRt + y/Z/, 7 + y2Z 2: ,-k 

\ i 


+ + y 4Z4J + N { 


(A4) 


The percent change over the 17-year period for the uniform model is given by: 

03 = 17-year change in percent (uniform model) 
17 x p 
Ave.lp-i} 


(A5) 


Trends expressed as percent per decade replace 17 years with 10 years in equation (5). 

Monthly weighting is done as described in the previous section. 

2.A.3 Regional Models (Total Ozone) 

The individual station trend results may be combined within regions (North America, Europe. Japan) 
to give estimates of regional trends. The methodology, given in Reinsel et al. (1987), calculates the regional 
trend as a weighted average of the individual station trends. The station weights are functions of the tren 
uncertainties at each station and the spread among the station trends, with stations with less precisely 
measured trends (e.g., from shorter records or more variable data) receiving less weight. 
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The ozone trends in each of the regions are expressed as a random effects model 


where is the trend estimate of the j' h station (either year-round or by season), 6, is a random term 

accounting tor station-to-station differences, and e, represents statistical variation in the trend due to within- 
station noise. 


2.A.4 Regional/Latitudinal Models (Total Ozone) 

The regional differences in the trends can be assessed simultaneously with a latitudinal dependence 
by using the model 


w,, 8„ + 8, ( L-45°N ) + «i + (tj, + e,j (A7) 

where co„ is the trend estimate at the j* station in the i' h region, L is latitude, «, is a random term representing 
region-to-region variation in the trends, is a random term representing station-to-station variation in the 
trends, and S„ + 8, (L-45°N) provides for a linear relationship between the trend estimates and latitude (L) 
where 8„ and 8, are the intercept and slope, respectively. Also, e„ is the statistical error in the i' h station 
trend estimate within the i ,h region. 

2.A.5 Empirical Aerosol-Adjustment Model (Umkehr) 

form- 10 lhC CmpiriCal aerosol ' ad -i ustment approach, the models for the Umkehr data are of the following 


y, - p. + S, + w x, + 7 , z, + 7 , p,, + N, (^g) 

where y,. t = I, .... T denote the monthly averages of Umkehr profile ozone data at a particular layer and 

a particular station. P, - e T * - I * - t, is a transformation of the composite optical thickness series t, 
and is similar to atmospheric transmission values (see Figure 2.3-5). The value / in (A 8 ) is a time shift in 
the monthly transmission (optical thickness) series over the period; for the period of study the values of/ 
used were / = I month for the station Belsk. / = 0 for Arosa, Edmonton, and Boulder, / = | for Lisbon 

Sapporo, Tateno. and Kagoshima, and / = -6 for New Delhi and Poona. The other terms in the model 
are S,, which represents the seasonal annual and semi-annual sinusoidal components, x„ which represents 
the F„„ solar flux series, and z„ which is the hypothesized deterministic linear trend through the period 

1977-1987. Thus, this model allows for a linear relation between the Umkehr data and the composite 

transmission (optical thickness) data during the period 1977-1987 involving the El Chichon volcano to 
account tor errors in the Umkehr data due to volcanic aerosols. For the station Kagoshima a simple 
intervention level shift term, starting in January 1986, was also included in the model to account for a 
possible change in mean level associated with a change of instrument at that time. (The data for Boulder 
also suggest a possibility of a mean level shift around the early part of 1982, possibly related to a change 
to the automated Dobson No. 61 instrument at the beginning of 1982. But such an effect has not been 
considered in the analysis because there is no definite confirmation of such a possible effect). Also for the 
station Lisbon, data were not used prior to 1978 because of a level shift in the data before that time due to 
a calibration change. As there were a substantial number of missing months in the Umkehr data series for 

most stations, for convenience an autoregressive process only of order one was considered for the noise 
term N,. 
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N, = <j> N,_i + e„ (A9) 

where the e, are independent with mean 0 and constant variance err. 

2.A.6 Theoretical Aerosol-Adjustment Model (Umkehr) 

In the theoretical model approach, the monthly average Umkehr data y, were first corrected for 
aerosol errors using the error correction factors supplied by DeLuisi et al. (1989a), which resulted in the 
corrected Umkehr data series denoted as y,*. It is noted that the error correction factors were applied to 
data at the various stations with the same time shift / as the value used in model (Al). For the corrected 
data, models of a similar form to (A8), but without the empirical aerosol error adjustment factor 
were considered: 


y,* = p, + S, + (ax, + 7 ,z, + N, (A10) 

2.A.7 Trends Model Applied to Ozonesonde Data 

For each of the 15 layers at each of the 9 stations the following regression model has been used: 

y, = |jl + S, + wx, + 5u, + N, (All) 

where y, is the average of ozone for month t, m is a constant term, S, is a seasonal component consisting 
of annual and semi-annual sinusoidal terms, wx, is a trend term, 8u, is a level shift term and N, is a noise 
term. The trend term wx, is such that x, = 0, t < T„ and x, = (t-TJ/12 for t > T„, where T 0 is December, 
1969; thus w is a parameter representing annual rate of change since 1970. The level shift term 8u, , where 
u, = 0 for t < T, and u, = 1 for t > T, with T, signifying the time and 8 representing the magnitude of the 
level shift, is used for all layers at the four Canadian stations to account for possible systematic changes 
in ozone readings due to changes of the sonde instrument from BM to ECC. It is also used for layers 1A— 
ID at Payerne with T, == April 1977 for possible effects due to the change of ozonesonde release time. 
Finally, the errors N, ’s are assumed to be normally and independently distributed with zero means but 
having different variances for the 12 months of the year. 
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